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Chapter 1 

Acceptance Test Plan Overview 


1.1 GMAT Introduction 


The.information presented in this Acceptance Test Plan document shows the current status of the General Mission 
Analysis Tool (GMAT). GMAT is a software system developed by NASA Goddard Space Flight Center (GSFC) 
in collaboration with the private sector. The GMAT development team continuously performs acceptance tests in 
order to verify that the software continues to operate properly after updates are made. The GMAT Development 
team consists of NASA/ GSFC Code 583 software developers, NASA/GSFC Code 595 analysts, and contractors 
of varying professions. 

GMAT was developed to provide a development approach that maintains involvement from the private sector 
and academia, encourages collaborative funding from multiple government agencies and the private sector, and 
promotes the transfer of technology from government funded research to the private sector. 

GMAT contains many capabilities, such as integrated formation flying modeling and MATLAB compatibility. 
The propagation capabilities in GMAT allow for fully coupled dynamics modeling of multiple spacecraft, in 
any flight regime. Other capabilities in GMAT include: user definable coordinate systems, 3-D graphics in any 
coordinate system GMAT can calculate, 2-D plots, branch commands, solvers, optimizers, GMAT functions, 
planetary ephemeris sources including DE405, DE200, SLP and analytic models, script events, impulsive and 
finite maneuver models, and many more. 

GMAT runs on Windows, Mac, and Linux platforms. Both the Graphical User Interface (GUI) and the GMAT 
engine were built and tested on all of the mentioned platforms. GMAT was designed for intuitive use from both 
the GUI and with an importable script language similar to that of MATLAB. 


1.2 Testing Methodology 


Purpose 


GMAT needs to undergo a series of rigorous tests to validate the numerical implementations of its models and 
establish a set of acceptable performance times. The 595 analysts created the acceptance test plan to achieve 
this goal by comparing GMAT with flight-operational reference software packages and documenting the results. 
Results can be reproduced with the initial conditions and software setups presented in this document. 
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CHAPTER 1. ACCEPTANCE TEST PLAN OVERVIEW 


Reference software 


For this comparative study to have merit, GMAT was tested against reliable, trustworthy, and flight opera- 
tional programs, such as STK-HPOP,STK-Astrogator, Free-Flyer, Swingby, and previous GMAT Builds that 
were comparable to the aforementioned programs. To achieve accurate comparison results, each program was 
compared with equivalent, or close to equivalent, test case setups. 


Testing Categories 


The Acceptance Test Plan divides into the following testing categories: Propagation, Calculation Parameters 
[Central body(Cb) and Coordinate System(CS) dependent], Integrators, Libration Points, Stopping Conditions, 
Delta V, and Performance. 


Scripts Used 


MATLAB scripts were created to make comparisons between GMAT Builds and the reference software. The 
majority of the comparisons involved taking the difference of the data and extracting the maximum absolute 
difference observed over the propagation duration. Scripts were also created to compare performance times 
for individual GMAT test cases to the reference software. The scripts created are as followed: Compari- 
son.Tooll_Tool2JPV.m, Comparison_Tooll_Tool2_CS .m, Comparison_Tooll_Tool2_Cb.m, ComparisonJnteg.m, 
TimeComparo.m, BuildRun_Script_GMAT.m, Comparison_Tooll_Tool2_Libr.m, ComparisomStopCond, and 
STKJEtepropagate.m. 

The user of the semi-automated scripts provides input when requested, in order, to perform the script’s core 
functions. For example, a user that wants to see the position and velocity differences between STK and GMAT 
would select a few choices from a menu. Next, the script would generate the comparisons based on the report 
data available. The semi-automation scripts adhere to the naming conventions outlined in their relevant testing 
category chapter. 

Most of the scripts generate output in at least one of the following formats: ASCII, LaTex, MATLAB -mat, 
or Excel .xls files. The report files are in an ASCII space delimited format and contain the different test case 
parameters outputted after propagation. The LaTex files contain the comparison data between two programs 
and provide an easy way to include that data into a PDF document. The .mat and .xls file are two other methods 
used to save the comparison data that proved useful from the software development team. 

The details of each script and how to use them are outlined in the relevant Testing Category section and/or the 
Comparison Scripts Guide section, located in Appendix C. 


1.2.1 Propagation 


The propagation test cases account for various orbits about Earth, as well as other celestial bodies. The main 
propagation parameters to monitor for differences are the position and velocity vectors. The following script was 
generated to perform the comparisons for this category: 

ComparisomTooll _Tool2_CS .m 

See the Propagators section (Chapter 2) for more detail and comparison results. 
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1.2.2 Calculation Parameters 
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The calculation parameter test cases verify the interna] calculations used to output the various parameters pre- 
sented in the list below. This section consists of two subsections: Coordinate System(CS) and Central Body(Cb) 
dependent parameters. The following scripts were generated to perform the comparisons for this testing category: 
Comparison-Tooll_Tool2_CS.m & Comparison_Tooll_Tool2_Cb.m 


<* Coordinate Systems 
— Earth Fixed 

— Earth Mean J2000 Equator (MJ2000Eq) 

- Earth Mean J2000 Ecliptic (MJ2000Ec) 
— Earth Mean of Date Equator (MODEq) 
— Earth Mean of Date Ecliptic (MODEc) 

— Earth True of Date Equator (TODEq) 

~ Earth True of Date Ecliptic (TODEq) 

— Earth Geocentric Solar Ecliptic (GSE) 

- Earth Geocentric Solar Magnetic (GSM) 
— Mars Fixed 

- Mars MJ2000Eq 
— Mars MJ2000Ec 
— Mercury Fixed 

- Mercury MJ2000Eq 
— Mercury MJ2000Ec 

- Moon Fixed 

— Moon MJ2000Eq 

- Moon MJ2000Ec 
— Neptune Fixed 

- Neptune MJ2000Eq 

- Neptune MJ2000Ec 
— Pluto Fixed 

- Pluto MJ2000Eq 

- Pluto MJ2000Ec 

- Saturn Fixed 
Saturn MJ2000Eq 

— Saturn MJ2000Ec 
— Uranus Fixed 

- Uranus MJ2000Eq 

- Uranus MJ2000Ec 

- Venus Fixed 

— Venus MJ2000Eq 

- Venus MJ2000Ec 

® Coordinate System Parameters 

- Position (X.Y.Z) 



— Velocity(X.Y.Z) 

— Magnitude of Velocity 

- Right Ascension of Velocity 

- Specific Angular Momentum 
— Argument of Periapsis 

- Declination 

— Declination of Velocity 
— Inclination 

- Right Ascension 

- Right Ascension of Ascending Node 
» Central Body Parameters 

— Altitude 

- Beta Angle 

- C3 Energy 
— Eccentricity 
— Latitude 

— Longitude 

- Specific Angular Momentum 

- Mean Anomaly 

- Mean Motion 

- Period 

— Apo apsis Radius 
— Perigee Radius 

- Position Magnitude 
— Semi-major Axis 

- True Anomaly 

- Semilatus Rectum 
— Apoapsis Velocity 
— Periapsis Velocity 

— Greenwich Hour Angle 

- Local Sidereal Time 

See the Calculation Parameters Section (Chapter 3) for more detail and comparison results. 
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1.2.3 Integrators 

The integrator test cases isolate the differences that would occur when changing the integrators for the same 
orbit. The following script was generated to perform the test case comparisons for tins category: 
Comparison_Integ.ni 


RungaKutta(RKV) 8(9) 
DormandElMikkawyPrince(RKN) 6(8) 




1.2. TESTING METHODOLOGY 

• RungeKuttaFehlberg(RKF) 5(6) 

• PrinceDormand(PD) 4(5) 

• PrinceDormand(PD) 7(8) 

• BulirschStoer(BS) 

• AdamsBashforthMoulton(ABM) 

See the Integrators Section (Chapter 4) for more detail and comparison results. 
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1.2.4 Stopping Conditions 

The stopping condition test cases determine how effective GMAT is at stopping satellite propagation on certain 
conditions. The following script was created to perform the test case comparisons for this category: 
Comparison-StopCond.m 

The stopping conditions tested are as followed: 

• Epoch (A1 Modified Julian Date) 

• Apoapsis 

• Elapsed Days 

9 Mean Anomaly 
» Periapsis 
® Elapsed Seconds 
» True Anomaly 
o XY Plane Intersection 
® XZ Plane Intersection 

• YZ Plane Intersection 

See the Stopping Conditions Section (Chapter 5) for more detail and comparison results. 


1.2.5 Libration Point 


The libration point test cases create data about the location of several libration points. Current and future 
satellite missions use libration points as part of their mission architecture. It is important to have accurate 
data for these libration points. The following script was created to perform the test case comparisons for this 
category: 

Comparison-Tooll_Tool2_Libr.m 

See the Libration Point Section (Chapter 6) for more detail and comparison results. 
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1.2.6 Delta V 


The delta v test cases determine the effectiveness of the delta v capabilities built into GMAT. When thruster 
burns are added to the mission sequence it is important that they are added correctly. The following script was 
created to perform the test case comparisons for this category: 

Comparison-DeltaV.m 

See the Delta V Section (Chapter 7) for more detail and comparison results. 


1.2.7 Performance 


The performance fast cases generate performance time data for later comparison between GMAT and the ref- 
erence software packages. Numerical calculation accuracy is important, but the amount of computing time it 
takes for the software to run is equally as important. We extracted several test cases from previous sections and 
ran them on the reference software packages, in order to check to make sure GMAT can perform just as good or 

better. 

See the Performance Section (Chapter ??) for more detail and comparison results. 


1.2.8 Control Flow 

The control flow tests generate report data that easily allows a Matlab script to produce a table of Pass and Fail 
cases.The following script was created to generate the Pass/Fail table for this category: 

LoopTestStimmary.nl 

See the Control Flow Section (Chapter 8) for more detail and results. 



Chapter 2 

Propagation 


In order to validate the accuracy of GMAT’s propagation, the fundamental unit-level components need to he 
combined and propagated on a system level. From a software development point of view, if the program under 
development is tested in a wide range of core applications, it is more likely that problems would be found before 
each new version is released to the public. This Acceptance Test Plan tests GMAT by comparing many possible 
scenarios users of GMAT would encounter to, reference software packages. Although it is impossible to create 
all the possible scenarios each user would encounter in GMAT, this is a start to eliminate possible frustrations 
a user could experience if a component did not work correctly. 

Propagation is one of the most important aspects of GMAT. Everything from outputting parameters to perform- 
ing a thruster burns at the correct stopping condition depend on whether or not GMAT is able to propagate the 
satellite/ object for a defined time period with acceptable accuracy. 

Parts of the Initial Orbit State Conditions section are referenced from Emergent Space Technologies’ Orbit De- 
termination Toolbox (ODTBX) Spiral 1 DEMO document, 1 due to a similar objective of testing the numerical 
implementation of the program’s base functions. 


2.1 Initial Orbit State Conditions 

2.1.1 Earth Based Test Cases 


The initial orbit states for the Sun-Synchronous (SunSync), Geostationary (GEO). Molniya, International Space 
Station (ISS) and the Global Positioning Satellite (GPS) orbits were obtained from Emergent Space Technolo- 
gies’ ODTBX Spiral 1 Demo - 1 Emergent used STK-High Precision Orbit Propagator (STK-HPOP) models and 
two-Line Element (TLE) sets with an initial UTC orbit epoch of June 1st 2004, 12:00:00:00. The initial orbit 
states that were used for the test case orbits can be seen in Table 2.1, on Page 24. The perigee and apogee 
altitudes of the test case orbits can be seen in Table 2.2, on Page 24. 


The propagation duration, report output step size, and integrator step sizes were varied for the different test 
cases. For the ISS, SunSync, GPS, Molynia, and GEO cases, the propagation length and report output step 
size were chosen based on a study performed by The Aerospace Corporation to validate STK-HPOP’s. 2 The 
integrator tune steps were chosen to allow for the most accurate comparison of the test case results. These time 
steps were based on Vallado’s analysis of state vector propagation. 3 
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CHAPTER 2. PROPAGATION 


Table 2.1: Satellite Initial Conditions 


Category 

Orbit Type 

X(km) 

Y (km) 

2 (km) 

Vx(km/s) 

Vy(km/s) 

Vz(lan/s) 

LEO 

ISS 

-4453.783586 

-5038.203756 

-426.384456 

3.831888 

-2.887221 

-6.018232 

LEO 

Sun-Sync 

-2290.301063 

-6379.471940 

0 

-0.883923 

0.317338 

7.610832 

MEO 

GPS 

5525.33668 

-15871.18494 

-20998.992446 

2.750341 

2.434198 

-1.068884 

HEO 

Molniya 

-1529.894287 

-2672.877357 

-6150.115340 

8.717518 

-4.989709 

0 • 

GEO 

GEO 

36607.358256 

-20921.723703 

0 

1.525636 

2.669451 

0 



Tabic 2.2: Apogee and Perigee Altitudes for Tost Satellites 


Orbit Type 

Perigee Altitude(km) 

Apogee Altitude (km) 

ISS 

358.168 

380.387 

Sun-Sync 

400 

400 

GPS 

19757.6 

20603.8 

Molniya. 

500 

39850.5 

GEO 

35786 

35786 


The chosen parameters can be seen in Table 2.3. on Page 24. 


Table 2.3: Integra, tor, Propagator, and Output Frequency 


Orbit Type 

Integrator Step Size(s) 

Propagator Length(days) 

, -V 

Output Frequency(mins) 

ISS 

5 

1 

1 

Sun-Sync 

5 

1 

1 

GPS 

60 

2 

2 

Molniya 

5 

O 

O 

5 

GEO 

60 

7 

10 


Several test cases were created for each satellite orbit to verify GMAT’s ability to perform accurately, while 
applying various forces. The forces used for Earth-based test cases were two-body, JGM2, EGM96, and JGM3 
gravity models, third-body perturbation effects from other planets, the Jacchia-Roberts (JR) and the Mass Spec- 
trometer and Incoherent Scatter Radar Exosphere (MSISE 1990 & 2000) Atmospheric Drag Model, and Solar 
Radiation Pressure (SRP). Each of the force models were run independently within GMAT to verify their indi- 
vidual accuracy, as well as a test case that includes an atmospheric drag model, the SRP model, a non-spherical 
gravity model, and third-body perturbations. These last test cases were performed to validate the capability of 
the GMAT to accurately propagate satellite orbits while multiple force models were applied. 


The Degree and Order for Earth-based non-spherical gravity cases was set at a constant 20 by 20. 


Refer to Appendix B.l for an alternate listing of all Propagator initial orbit state conditions. 
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OTHER INITIAL STATE CONDITIONS 



X 1 j 


2.1.2 Non- Earth Based Test Cases 




GMAT is designed for accuracy in non-Earth mission scenarios. Test cases for Mars, Mercury, the Moon, Neptune, 
Pluto, Saturn, Uranus, Venus, L2 orbits, and deep space orbits were created to test various forces, individually 
and jointly, affecting a spacecraft in an orbit. Many satellite parameters, such as Cd,Cr, satellite area, and 
satellite mass, were kept the same as the Earth test cases for simplicity and consistency. The initial Keplerian 
satellite state only varied in Semi-Major Axis and gravity field degree & order for the non-Earth test cases. 


Refer to Table 2.4 for the initial Keplerian orbital elements and the Degree & Order used for the non-spherical 
gravity force cases. The integrator step size, propagation length, and output frequency for all the non-Earth 
cases are 5 seconds, 3 days, and 5 minutes, respectively. Table 2.4 excludes the deep space and L2 orbit test cases. 


Table 2.4: Non-Earth Keplerian Orbital Elements 


Orbit Type 

SMA (Ion) 

Ecc. 

Inc. (deg) 

AOP (deg) 

RAA.N (deg) 

TA (deg) 

Deg. x Ord. used 

Mars 

4603 

0.2 

45 

45 

90 

45 

20x20 

Mercury 

3640 

0.2 

45 

45 

90 

45 

4x.O 

Moon 

2500 

0.2 

45 

45 

90 

45 

20x20 

Neptune 

34999 

0.2 

45 

45 

90 

45 

4x0 

Pluto 

1795 

0.2 

45 

45 

90 

45 

4x0 

Saturn 

80000 ' 

0.2 

45 

45 

90 

45 

4x0' 

Uranus 

45000 

0.2 

45 

45 

90 

45 

4x0 

Venus 

8125 

0.2 

45 

45 

90 

45 

20x20 


SMA: Semi-Major Axis | Ecc.: Eccentricity | Inc.: Inclination j AOP: Argument of Perigee 
RAA.N: Right Ascension of Ascending Node | TA: True Anomaly 


The DeepSpace, Earth Moon L2 (EML2), and Earth Sim L2 (ESL2) cases involve propagating about libration 
points and propagating deep space orbits. Table B.14, B.15, and B.16, in Appendix B, provide the initial 
states for the DeepSpace, EML2, and ESL2 test cases, respectively. 


Refer to Appendix B.l for a listing of all Propagator initial orbit state conditions. 


2.2 Other Initial State Conditions 


In order to reduce the complications of the comparisons, certain initial orbit parameters were kept constant 
throughout all of the cases. These parameters are Cd, Cr, Spacecraft Area, each programs integrator, and 
software settings affecting the results of various force models. 

The GMAT integrator used for all the test cases was Runga Kutta 8(9), except for the STK-HPOP test cases, to 
avoid any additional differences that could occur from changing integrators. The Integrators Section compares 
the differences between the various integrators GMAT can use. 

The parameters in Table 2.5 shows the differences between GMAT and the reference programs. The ideal situation 
would be for all the programs to match perfectly, but that is not realistic due to the different approaches each 
program takes to solve each problem. 
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i&Die z.o. umv 

Parameter 

ersai * ■ 

GMAT 

STK-HPOP 

FF 

STK- Astro 

Gd 

2.2 

2.2 

2.2 

2.2 

Cr 

1.2 

1.2 

1.2 

1.2 

Area (m 2 ) 

20 

20 

20 

20 

Satellite Mass (kg) 

1000 

1000 

1000 

1000 

Integrator (excluding integrator test eases ) 

, RK8(9) 

RK7(8) 

RK8(9) 

RK8(9) 

Integrator error tolerance 

lc-013 

lc-013 

lc-013 

lc-013 

Integrator error control. Relative to 

Step 

Step? 

? 

Step 

Drag: Altitude Calculation 

Approximate 

Approximate 

Exact 

*7 

Drag: Sun. Position 

True 

True 

True 

True 

SRF: Sun Position 

True 

True 

True? 

True 

Solar Flux ( W/m 2 ) at 1 AU 

1359.38857 

1359.38857 

1358 

1359.38857 

Solid Tides 

N/A 

Disabled 

N/A 

N/A 

Ocean Tides 

N/A 

Disabled 

N/A 

N/A 

Dadlv F10.7: JR, and MSISE models 

150 

150 

N/A 

150 

Average F10.7 

150 

150 

150 

150 

Geomagnetic Index(Kp): JR and MSISE only 

3 

3 

N/A 

3 

Drag: Geomagnetic Flux Update 

Constant 

Constant 

N/A 

Constant 

Boundary Mitigation 

N/A 

Disabled 

N/A 

N/A 

Relativistic Accelerations 

N/A 

Disabled 

N/A 

N/A 

Shadow Modeling 

Dual Cone 

Dual Cone 

Dual Cone 

7 

JERS EOF format used: 

Long term C04 

Bulletin A 

Bulletin A 

Bulletin A 

Polar Motion calculation: 

Enabled 

Enabled 

Enabled 

Enabled? 

Nutation undatc interval: Earth (see) 

60 

60 

60 

50 

Planetary Ephemeris update interval (sec) 

0 

0? 

0 

0? 
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OTHER INITIAL STATE CONDITIONS 



Exceptions to Table 2.5 are as followed: 


• STK has trouble propagating the EML2 test case at an error tolerance of le-013 with a relative to step 
error control. The GMAT and STK cases were changed to both use relative to state error control 


2.2.1 Earth Orientation Parameters (EOP) data 

’’International Earth Rotation Service (IERS) Bulletins A and B provide current information on the Earth’s 
orientation in the IERS Reference System. This includes Universal Time, coordinates of the terrestrial pole, and 
celestial pole offsets. Bulletin A gives an advanced solution updated weekly by e-mail subscription or daily by 
anonymous ftp; the standard, solution is given monthly in Bulletin B and updated every week in the (IERS) C04 
solution. ” 6 

’’Bulletin A is issued by the IERS Rapid Service/Prediction Centre at the U.S. Naval Observatory (USNO), 
Washington, DC and Bulletin B”, as well as the C04 data, ”is issued by the IERS Earth Orientation Centre at 
the Paris Observatory. ” 4 

’’Bulletin A is intended for users who need accurate information before the Bulletin B finals series is available, 
i.e., those who reduce data in the very recent past ( require rapid service) or those who operate in real-time 
(require predictions). Bulletin B is intended for standard use. For scientific and long-term analyses of the 
Earth’s orientation, the long-term continuous series” A C04 (1962- present), can be used. 

”EOP (IERS) C04 is regularly recomputed to take advantage on one hand of the improvement of the various 
individual contributions and in the other hand of the refinement of the analysis procedures. To date, it is twice- 
weekly updated. ” 5 

’’The EOP (IERS) C04 is given at one-day intervals, it is free from the diumal/subdiumal terms due to the 
oceanic effects and can be interpolated linearly. The oscillations in UT and duration of the days due to zonal 
tides for periods under 35 days are present in full size in the series.” 5 

GMAT retrieves long term earth orientation IERS EOP C04 data, which includes UTC-UT1 data, from a file. 
This file includes smoothed values at 1-day intervals and data from 1962-present. 

STK and FF retrieves its EOP data from the USNO series 7 / IERS Bulletin A. 

The differences between the EOP data sets are displayed in Table 2.6. The Terrestrial Pole column refers to the 
accuracy of the pole position [x.y] and the UT1 column refers to the accuracy of the rotation angle about the 
pole UT1. 

As shown in Table 2.6, there are differences between STK and GMAT, but the Terrestrial Pole data agrees to 
within the thousandth place of a milli-arcseconds and the UT1 data agrees to within the hundredth place of 
second. 
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Table 2.6: EOP Format Accuracy 


EOP Format 

Pole Position (mas) 

UT1 (ms) 

Bulletin A obs. 1-d (1) 

0.10 

0.02 

Bulletin A. pred. 1-d (2) 

0.50 

0.14 

Bulletin A pred. 4-d (2) 

1.60 

0.52 

Bulletin A pred.. IQ-d (2) 

3.0 

1.60 

Bulletin A pred. 40-d (2) 

11.2 

7.70 

Bulletin B obs. smooth 5-d (3) 

0.15 

0.02 

Bulletin B obs. raw 5-d (3) 

0.15 

0.02 

Bulletin B pred. 5-d (3) 

1.6 

0.60 

Bulletin B pred. 10-d (3) 

3.0 

1.60 

Bulletin B pred. 30-d (3) 

10.0 

4.0 

Long-Term C04 =62-=67 

30.0 

2.0 

Long-Term C04 =68-71 

20.0 

1.5 

Long-Term C04 =72-=79 

15.0 

1.0 

Long-Term C04 =8G- : 83 

2.0 

0.4 

Long-Term C04 ’84-’95 

0.7 

0.04 

Long-Term C04 ’96-present 

0.2 

0.02 


NOTES: (1) Based on data after 1997: applies only to latest epochs in each update. 


(2) Based on data since 1995. 

(3) Based on data since 1996. 

The Terrestrial Pole and UT1 data is free from the diurnal/ subdiumal 
terms due to the oceanic effects and can be interpolated linearly. 

These terms can be added after interpolation. 
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2.2.2 Other Planetary Parameters 
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Gravitational Constant 


All programs used for comparisons utilize the DE405 Planetary Gravitational Constants listed in Table 2.7, 
except when non-spherical gravity was used. When non-spherical gravity files are used they typically call the 
gravitational constants located in the file, unless the program creates an exception. 


Table 2.7: Planetary Gravitational Constants (mu values) 


Planet 

mu (km? /s' 2 ) 

Sun 

132712440017.99 

Mercury'- 

22032.080486418 

Venus 

324858.59882646 

Earth 

398600.44150000 

Moon 

4902.8005821478 

Mars 

42828.314258067 

Jupiter 

126712767.85780 

Saturn 

37940626.061137 

Uranus 

5794549.0070719 

Neptune 

6836534.0638793 

Pluto 

981.60088770700 


Flattening Coefficient 


Whenever possible the flattening coefficients listed in Table 2.8 were used. Without the use of these values the 
planetary bodies could have wildly different shapes, which would result in large differences in parameters such 
as longitude, latitude, and altitude. 


Table 2.8: Flattening Coefficient 


Planet 

Flattening Coefficient 

Sun 

0.00000000 

Mercury 

0.00000000 

Venus 

0.00000000 

Earth 

0.00335270 

Moon 

0.00000000 

Mars 

0.00647630 

Jupiter 

0.06487439 

Saturn 

0.09796243 

Uranus 

0.02292734 

Neptune 

0.01856029 

Pluto 

0.00000000 
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Equatorial Radius 

The several celestial body equatorial radii used are listed in Table 2.9. Similar differences as the flattening 
coefficient occur when GMAT and the reference programs don’t use the same values. 

Table 2.9: Equatorial Radius 


Planet 

Radius (km) 

Sun 

695990.0 

Mercury 

2439.7 

Venus 

6051.90 

Earth 

6378.1363 

Moon 

1738.2 

Mars 

3397.0 

Jupiter 

71492.0 

Saturn 

60268.0 

Uranus 

25559.0 

Neptune 

25269.0 

Pluto 

1162.0 


Leap Seconds 

The amount of leap seconds, or A AT, has been used since 1972 in order to keep — UTC-UT1 — < 0.9sec. GMAT 
and the reference software packages use all the leap seconds up until 2004. In 2004 the amount of leap seconds 
in use were 32 seconds. 



2.3. NAMING CONVENTION 

2.3 Naming Convention 


This section describes the naming convention for propagator scripts and output reports. The naming convention 
consists of an ordered series of option strings, separated by underscores (. ). Currently, options are allowed for 
the following fields, and will be present in the file name in order: 


1. tool - The tool used to generate the test case 

2. traj - The trajectory to use. This includes initial conditions, physical parameters, and time step 

3. pmg - The point-mass gravity model to use 

4. nsg - The non-spherical gravity model to use 

5. drag - The atmospheric drag model to use 

6. other - Any other forces to include, such as SRP, secondary body gravity, etc 

The tool field should always be the first field. Future additional fields should be added to the end of the list 
of fields. If multiple other options are required, they should be added to the end of the file as required. For 
example, the file name will be tooLtraj-pmg-nsg-drag-otherNother2.report (file extensions are described later.) 
Each field has a finite list of options, as follows (future options should be added to this list): 


1. tool 

STK - Satellite Toolkit HPOP or Astrogator 
FF - FreeFlyer 

GMAT - General Mission Analysis Tool 


2. traj 


ISS 

- LEO orbit 

S unSync 

- LEO orbit 

GPS 

- MEO orbit 

GEO 

- GEO orbit 

Molniya 

- HEO orbit 

Marsl 

- eccentric low orbit 

Mercuryl 

- eccentric low orbit 

Moon 

- eccentric low orbit 

Neptimel 

- eccentric low orbit 

Pluto 1 

- eccentric low orbit 

Saturnl 

- eccentric low orbit 

Uranusl 

- eccentric low orbit 

Venus 1 

- eccentric low orbit 

DeepSpace 

deep space orbit 

EML2 

- Earth Moon L2 orbit 

ESL2 

- Earth Sim L2 orbit 


NOTE: Some test cases contain traj variations. In these cases, traj precedes the modification. For example, 
if ISS trajectory is needed with no output, then traj can be ISSnoOut. The lack of a report file is shortened 
to noOut. 


3. pmg 
Earth 
Sim 
Luna 
AllPlanets 


- Earth point mass gravity 

- Sun point mass gravity 

- Lunar point mass gravity 

- Sun, Mercury, Venus, Earth, Moon, Mars, Mercury, Jupiter, and Pluto point mass gravity included. 


NOTE: When dealing with a combination of pmg , s the first point mass is the primary body and the 
following are third body point masses. For example, LunaSunEarth would be a Lunar primary body with 
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the Earth and Sun as third body point masses. The pmg’s after the primary body are arranged based on 
the order from the sun, in order to reduce repeat filenames. 


4. nsg 

0 - no non-spherical gravity included 

JGM2 - Earth JGM2 20x20 gravity 

JGM3 - Earth JGM3 20x20 gravity 

EGM96 - Earth EGM96 20x20 gravity 
MARS50C - Maxs Mars50c 20x20 gravity 
LP165P - Moon LP165P 20x20 gravity 


5. drag 

0 - drag not included 

HP - Harris Priester 

JRXX - Jacchia-Roberts 

MSISEXX - NRL MSISE 

NOTE: XX in the drag field refers to the year. For example. JR77 would be the Jacchia-Roberts 1977 
model, and MSISE00 would be NRL MSISE 2000. Refer to Table 2.5 for the drag settings used. 


6. other 

0 - no other forces included 

SRP - Solar Radiation Pressure 

NOTE: Any of the above options may be included as an other field. Refer to Table 2.5 for the SRP settings 
used. 


2.3.1 Comparison Script Information 


The script used to perform the position and velocity comparisons needed for the Propagator section is Compar- 
ison-Tooll_Tool2JPV.m. This script takes the normalized position and velocity vector difference between two 

programs. 

Refer to Appendix C for more details about this script and others used in the Acceptance Test Plan document. 


2.4 Test Case Results 


The following results are for the Propagator section. The current GMAT Build is compared to STK and FreeFlyer 
for this section, with the maximum normalized position and velocity difference displayed in table format. 

To determine if a propagator test case comparison value was acceptable, an acceptance matrix, presented in 
Table 2.10, was created. The values in Table 2.10 were obtained from the lower position difference bounds of 
David Vallados An Analysis of State Vector Propagation Using Differing Flight Dynamics Programs, presented 
at the 2005 American Astronomical Society (AAS)/AIAA Astrodynamics Specialist Conference. These lower 
boxinds are difficult to meet in some orbits due to the wide range of orbits that are possible but they give a order 
of magnitude number to strive for. If a case has a combination of either Drag, Non-Spherical Gravity, Solar 
Radiation Pressure(SRP), or Point Mass gravity, the largest acceptable position difference is used. 

The next step beyond this acceptance matrix is to compare GMAT’s comparison data to differences seen in FF 
and STK comparisons, and peer reviews. 
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2.4. TEST CASE RESULTS 


Tab le 2.10: Acceptance Matrix 


Difference in 

Acceptable Position Difference (m) 

Non-Spherical Gravity 

< 0.001 

Point Mass Gravity 

< 0.001 

Solar Radiation Pressure 

< 0.6 

Drag 

<20 



Table 2.11: GMAT/STK GEO STK Teat Case Comparison 


Test Case 

Position Differcnce(m) 

Velocity Differcnco.(irv/s j 

AUPlanets-0-0-0 

0.0002573600464 

1.793784749e-008 

EarthLuna-0-0-0 

0.000319469177 

2.246260344e-008 

EartliSunLumvEGM9C-JR-SRP 

0.0001571119SG 

1.10'1136958e-008 

EaithSunLuna-EGM96-MSISE90-SRP 

0.0001571119996 

1.101136862e-008 

EarthSunLun;vJGM2-JR-SRP 

0.0001132816464 

7. 859884531 c-009 

E artiiS unL una- JGM2-MS ISE90- SRP 

0.0001132816525 

7.859884883C-000 

EarthSunLuna-JGM3-JR-SRP 

0.0002438654858 

1.730633303e-008 

E ar thS unL una- JGM3- M SISE90- SRP 

0.0002438654861 

1.730633311e-008 

EarthSun-0-0-0 

2.2791 68371e-005 

1 . 61446 8 7 6f>e-00 9 

Earth-0-G-Q 

2.111 395905e-005 

1.410324921e-009 

Earth-0-0-SRP 

4.571 760608e-005 

2.618694012(3-009 

Earth-O-JR-O 

2.111395905e-005 

1.410324.921e-009 

Earth-0-MSISE90-0 

2.111 395905e-005 

1.410324Q21e-009 

Earth-E GM 96-0-0 

0.0001110202624 

8.02931 0693c-00 9 

Earth-JGM2-0-0 

2.954895496e-005 

2.055988066-009 

Earth- JGM3-0-0 

1.603674008e-005 

1.08902881e-009 


Table 2.12: GMAT/STK GPS STK Test Case Comparison 



Posi tion Difference(m) 

Veloci ty Differen ee (m/s) 

AllPlancts-O-O-O 

1.857890 30 3c-005 

2.696491985e-009 

Earth! una- 0-0-0 

2.280562281e-005 

3.320021018e-009 

EarthSunLuna-EGM96-JR-SRP 

0.1959052443 

2.485425413(3-005 

EarthSunLuna-EGM96-MSISE90-SRP 

0.1959347769 

2.485881465e-005 

EarthSunLuna-JGM2-,lR,-SRP 

0.1959143599 

2.485563669e-005 

EarthSunLuna-JGM2-MSISE90-SRP 

0.1959227868 

2. 485702347c- 005 

EarthSunLunar-JGM3-JRrSRP 

0.1959136649 

2.485560304e-005 

EarthSunLuna-JGM3-MS'ISE90-SRP 

0.1959104274 

2.4855l3888e-005 

EarthSun-0-0-0 

1.81527305O-006 

2.25 74l6233e-0 10 

Earth- 0-0-0 

2. 704817671 c-00(j 

3. 741 292798o-010 

Earth-O-O-SRP 

0.0878120195 

9.617041298e-006 

Earth-O-JR-O 

2.7048L7671e-006 

3.741292798e-010 

Earth-0-MSISE90-0 

1 .54909103e-005 

2.226217304e-00.9 

Earth-EGM96-0-0 

3.230673797e-00-5 

4.2601 88402e-009 

Earth- J GM2-0-0 

3 . 20 785 9422e-005 

4.29991 2419e-009 

Earth-JGM3-0-0 

3.197574489(3-005 

4.255726699(3-00 9 
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Table 2.13: GMAT/STK ISS STK Test Case C omparison 


Test Case 

Position Difference(m) 

Vel o city Difference (m/s ) 

AllPianelS-O-O-O 

1.236533809e-005 

1.385493045e-008 

E arthLrma- 0-0-0 

2.65638111 7e-005 

3.027548132e-008 

EarthSimLuna^EGM96-JR-SRP 

278.7490457 

0.3182017608 

EaitliSmiLuna-EGM96-MSISE90-SRP 

5G.19307014 

0.06443183746 

EarihSrmLuna-JGM2-JR-SRP 

265.043753 

0.3028869043 

EarthSunLuua- J GM2-M SISE90-SRP 

56.19388165 

0.06443206028 

EarthS unLuna-JGM3-JR-SRP 

278.7489221 

0.3182016489 

EaxthSunLuna-JGM3-MSlSE90-SRI J 

56.19332541 

0.06443145135 

EarthSon-0-0-9 

1.433764674e-005 

1.600957323«-008 

Earth-0-0-0 

2.047203147e^005 

2.325251323e-008 

Earth-O-O-SRP 

0.244823116 

0.0002550285075 

Earth-0- JR-0 

251.765217 

0.2867219205 

Earfcb-0-MSI3E90-0 

54.12970949 

0.06180127344 

Eaxth-EGM96-0-0 

0.0003663016096 

4.164070917C-007 

Eaxth-JGM2-0-0 

0.0002104591217 

2.411877645C-007 

Earth-JGMTO-O 

0.0002166929687 

2.4750508(ie-007 


Tabic 2.14: GMAT/STK Moiniya STK Test Case Comparison 


Test Case 

Position Difference (m) 

Velocity Diffcrcncc(m/s) 

AIlPIanets-0-0-0 

0.0001645654561 

1 .25215 7429e-007 

EartliLunar 0-0-0 

0.0001587779228 

1.205662193e-007 

Earf,h.SuiiLnaa-EGM96- JR-S.RP 

13.53119687 

0.00.9730652052 

EarthSmiLuna-EGM96-MSISE90-SRP 

7.25773273 

0.004949260962 

E arthSuiiLuiia- JGM2- JR-SRP 

13.53346199 

0.009731885706 

EarthSunLuna-JGM2-MSISE90-SRP 

7.260932491 

0.00495105133 

EarthSimLuna- J GM.3- JR,- SRP 

13.53295584 

0.009731601686 

EarthSunLuna-JGM3-MSISE90-SRP 

7.260582648 

0.004951077623 

EarthSun-0-0-0 

0.000348717827 

2.914208443e-007 

Earfch-0-0-0 

0.0002791429809 

2.328568085e-007 

Earth-0-O-SRP 

0.571576837 

0.00047932174-59 

Earth-0- JR-0 

15.70095053 

0.01312840386 

Earfch-0-MST.SE90-0 

7.009963698 

0.005861407998 

E aith-EGM96-0-0 

0.001747777921 

1.46609522e-006 

Earth- JGM2-0-0 

0.00156516582 

1 . 313722882e-006 

Earth-JGM3-0-0 

0.001286385724 

1.0808574626-006 
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Table 2.15: GMAT/STK SunSync STK Test Case Comparison 


Test Case 

Position Difference(m) 

Velocity Difference(m /s) 

AllPlaaicts-0-0-0 

4.699222822O-005 

5.325204802o-008 

EarthLuna-0-0-0 

1.738032476e-005 

1.943588041e-008 

EarthS unLima-EGM'96-JR-SRP 

185.5519514 

0.2108884067 

EanhSimLima-EGM96-MSISE90-SRP 

45.12786362 

0.05132420504 

Eai-thSuuLunarJGM2-,TR-SRP 

185.5519345 

0.2108883121 

EarthSunLuna-JGM2-MSISE90-SRP 

45.12684175 ‘ 

0.05132310311 

E arthS iml. un a- J GM3- JR-SRP 

185.5517129 

0.2108881509 

EarthS unLuna- JGM3-MSISE90- SRP 

45.12759811 

0.0513239131 

EarthS un-0-0-0 

4.687777606O-006 

5.023Q88912c-009 

Earth-0-0-0 

4.483249915C-005 

5.068280737c-008 

Earth-O-O-SRP 

0.1578758295 

0.0001369580266 

Earth-0- JR-0 

172.0259163 

0.1950745493 

Earth-0-MSlSE90-0 

42.66516434 

0.04839009284 

Eaith-EGM96-0-0 

7.308356983e-005 

8.466691171(3-008 

Earth- .J GM2-0-0 

8 . 928348594e-005 

9.655432558e-008 

Earth- JGM3-0-0 

6 . 997911 597e-005 

7.571452902e-008 
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Table 2.16: GMAT/STK Marsl STK Test Case Comparison 


Test Case 

Position Difference(m) 

Velocity Differ ence/m/s) 

AUPlanots-0-0-0 

0.1994303898 

0.0001681334573 

Mars-0-0-0 

0.1873874394 

0.0001577557948 

Maxs-O-O-SRP 

0.7867928767 

0.0006572167727 

Maxs-MARS50C-0-0 

0.4488316637 

0.0003800813667 

Maxs-MARS50C-0-SRP 

1.214053862 

0.001021133703 


Tabic. ‘2.17: GMAT/STK Mcrcuryl STK Tost Cose Comparison 


Test Case 

Position D ifference(xu ) 

Velocity Difference(m/$) 

AUPlanets-0-0-0 

Mercury-0-0-0 

Mercury-O-Q-SRP 

0.1124163711 

0.02685481636 

58.83805392 

9.591759819^005 

2.21S352736e-005 

0.05048777188 


Table 2.18: 

GMAT/STK. Moon STK Test Case Comparison 

Test Case 

Position Diiferenee(m) 

Velocity Difference(m/s) 

AlIPlancts-0-0-0 0.0002140368472 

Lunar 0-0-0 0.03443001175 

Luna-0-0-SRP 0.03471809841 

LunarL.P 165P-0-0 0.0002701 673647 

Lima-LP 165P-Q-SRP • 0.1350043985 

1.35 5643864c-007 
2.456930533&-005 
2.41327l946e-005 
1 . 9221 74.669e-007 
9.596583571 e-005 


Table 2.19: GMAT/STK Neptunel STK Test Case Comparison 

Test Case 

Position Diifcrcucc{m) 

Velocity Diffcrcnce(m/s) 

AllPlanets-0-0-0 
Neptune- 0-0-0 
N ep tune-0-0- SRP 

0.9974264518 

1.171036908 

0.6214495234 

0.000507487521 5 
0.0005920828399 
0.0003128000507 

Table 2.20: GMAT/STK Plutol STK Test Case Comparison 

Test Case 

Position Difference(m) 

V sloci ty Difference (m/s) 

AilPlanets- 0-0-0 
Plato-0-0-0 
Pluto-0- 0-SRP 

0.8883195792 

0.2632397208 

0.7534556636 

0.0004683900657 

0.0001382385307 

0.0003932329143 



2.4, TEST CASE RESULTS 


Table 2.21: GMAT/STK Saturn! STK Test Case Comparison 


Test Case 

Position Differenee(m) 

Velocity Diifereuce(m/s) 

ALlPlanets-0-0-0 

0.14.14864331 

4.873107384e-005 

Saturn-0-0-0 

0.4497282767 

0.0001561546333 

Saturn-O-O-SRP 

0.3124746558 

0.000102600788 


Table 2.22: GMAT/STK Uranusl STK Test. Case Comparison 


Test Case 

Position Difference(m) 

Velocity Difference (m /s) 

AllPlanets-0-0-0 

0.2465709269 

7.912323063e-005 

Uranus- 0-0-0 

1.320413817 

0.0004243963023 

Uranus- 0-0-SRP 

1.156382201 

0.0002904419051 


Table 2.23: GMAT/STK Venusl STK Test Case Comparison 


lest C&se 

Position Difference (m) 

Velocity Difference (m /» } 

AllPlanets-0-0-0 

0.02561929021 

2.548815625e-005 

Vcmis-0-0-0 

0.01793453836 

1.768517227O-005 

Venus-0-0-S'RP 

0.3477682342 

0.0003453235732 

Venus-MGNP 180U-0-0 

0.007510965788 

7. 525 754724e-005 

Vein ls-M GNP 1 SOU- ( )- SRP 

0.4083837321 

0.0004049592692 
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Table 2.24: GMAT/STK DeepSpace STK Test Case Comparison 


Test Case 

Position Difference(m) 

Velocity Difference (m/s) 

AUPlancts-O-O-O 

0.02027141264 

3.996534351o-009 


Table 2.25: 

GMAT/STK EML2 STK Test Case Comparison 

Test Case 

Position Difference(m) 

Velocity Difference (m/s) 

AllPlanets-0-0-0 522765.4104 

AUPlanets-O-O-SRP 162722.736 

EarthSimLima-0-0-0 191050. 7674 

E arthSirnLuna- JGM 2- 0- 0 83170.98289 

2.863727549 

0.8844521611 

1.043593205 

0.4533442704 


Table 2.26: GMAT/STK ESL2 STK Tost Case Comparison 


Test Case 

Position Difference(m) 

Velocity Difference (m /s j 

AllPlanets-0-0-0 

29036.76192 

0.01713687194 

AllPlauets-0-0-SRP 

228360.2465 

0.1774325007 



2.4. TEST CASE RESULTS 
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Table 2.27: FF/GMAT GEO GMAT Test Case Comparison 


Test Case 

Position Difference (in) 

Velocity Difference (m/s) 

AIlPlancis-0-0-0 

0.02466296086 

2.348952399C-006 

EarthLuna-0-0-0 

0.02466037174 

2.264158243^006 

EarthSun-0-0-0 

4. 797 553846-005 

8.493947765O-007 

Earth-0-0-0 

5.214381207e-005 

6-896273096-007 

Earth-0-0-SRP 

2.899321865 

0.000124521 9489 

Earth-JGM2-0-0 

0.02515954067 

2.313999898e-006 


Tabic 2.28: 

FF/GMAT GPS GMAT Tost Case Comparison 

Test Case 

Position Difference(m) 

Velocity Difference (rn /s) 

AUPlanets-0-0-0 

EarthLmia-0-0-0 

EarthSuu-O-O-G 

Earth-0-0-0 

Earth-0-0-SRP 

Earth- JGM2-0-0 

0.001885974435 
0.001874007099 
3.4321 7 6702c-006 
5.489921041e-006 
0.5814667923 
0.01104317452 

1.0520207056-006 
9. 43361 654&-007 
8.317293859O-007 
7.9829921e-007 
6.084184505e-005 
2.140216499e-006 

Table 2.29: 

FF/GMAT ISS GMAT Test Case Comparison 

Test Case 

Position Difference ( m) 

Velocity Difference (m/s) 

AIlPlanets-0-0-0 

EarthLunarO-O-O 

EartliSun-0-0-0 

Earth-0-0-0 

Eartb-0-0-SRP 

Earth-JGM2-0-0 

5.610 904349e-006 
1.1334898246-005 
9.638402495e-006 
2.358180759e-005 
0.1095353597 
0.2076314901 

8.2240147966-007 

8.407000067e-007 

8.134703987e-007 

8.1864684346-007 

0.0001183226844 

0.0002404721693 


Table 2.30: FF/GMAT Molniya GMAT Test Case Comparison 


Test Case 

Position Difference (m) 

V elocity Difference (m/s) 

AUPlanets-0-0-0 

0.01444526657 

1.2033960446-005 

E art hL una-0-0-0 

0.01424235005 

1.19546892e-005 

EarthSim-0-0-0 

0.0003311394147 

8. 2326473 736-007 

Earth-0-0-0 

5.789754939e-005 

8.078947532e-007 

Earth-0-0-SRP 

2.735360369 

0.002258471401 

Earth- J GM2-0-0 

3.064693947 

0.002562937467 
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Table 2.31: FF/GMAT SunSync GMAT Test Case Comparison 

“Test Case Velocity Diffc.rciicc(m/s) 


AllPlanets-0-0-0 

EarfchLiuxarO-O-O 

Earth.Sun-0-0-0 

Earth-O-O-G 

Earth-O-O-SRP 

Eaxth-JGM2-0-0 


1.6244528886-005 
3.6419721526-005 
1.621 851 69b-005 
8.<)3851495o-006 
0.04776731849 
0.09458368297 


8.363014362e-007 
8. l39327674e-007 
8.4707995386-007 
8.215195c-007 
4.87371 7267C-005 
0.0001055216409 


2.4. TEST CASE RESULTS 


Table 2.32: FF/GMAT EML2 GMAT Test Case Comparison 


Test Case 

Position Difference(m) 

Velocity Difference(m/s) 

AllPlaiiets-0-0-0 

732435.6987 

4.007738969 

AllPlanets-O-O-SRP 

370563.5763 

2.010501473 

EarfchSunLuna- 0-0-0 

400326.9215 

2.185454594 

E ar thSi mLi ma- -I GM2- 0 - 0 

292445.684 

1.595241714 


Table 2.33: FF/GMAT ESL2 GMAT Test Case Comparison 


Test Case 

Position Difference (m) 

Velocity Difference (m/s) 

AllPlanets-0-0-0 
A11P1 anets-O-O-SRP 

5690806.264 

68812.33.15 

3.131093893 

5.36031143 
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2.5 FF/STK Comparison 


Table 2.34: FF/STK GEO S TK Test Case Comparison 


Test Case 

Position Difference(m) 

V elo city Difference (m/s) 

AllPlanets-0-0-0 

0.02441863704 

2.331184604e-006 

EaxthLuaa-0-0-0 

0.02435631273 

2.242954572e-006 

EarthSunLuna- J GM2-HP-SRP 

2.852626845 

0.0001210061779 

EaxthSun-0-0-0 

4.231055926e-0(J5 

8.493743904e-007 

Eaith-0-0-0 

6 . 836800 60 8e-00 5 

6.889879736-007 

Eaxth-O-O-SRP 

2.899330477 

0.0001245232993 

Earth-O-HP-O 

6.836800608O-005 

6.88987973C-007 

Earth,- J GM2-0-0 

0.02513336324 

2.312005795O-006 


Table 2.35: 

FF/STK GPS STK Test Case Comparison 

Test Case 

Position Difference(m) 

Velocity Difference (m/s) 

AllPlanets-0-0-0 

0.001903432726 

1.05396749o006 

EaxthLuna-0-0-0 

0.001896072444 

9.455905846e-007 

EaxthSunLima-JGM2-HP-SRP 0.6785588953 

6.968633054e-005 

EarthS un-0-0-0 

3.872950958e-006 

8.3163744686-007 

Earth-G-0-0 

8.l64858288e-006 

7.983002389e-007 

Eartli-0- 0-SRP 

0.6692784196 

6.82998191 e-005 

Earth-O-HP-O 

8.164858288e-006 

7.983002389e-007 

Earth-JGM2-0-0 

0.01104063243 

2.139207943O-006 


Table 2.36: 

FF/STK ISS STK Test Case 

Comparison 

Test Case 

Position Difference(m) 

Velocity Difference (m/s) 

AllPlauets-0-0-0 

1.17832665e-005 

8.2741934976-007 

EaxthLunarO-O-O 

2.023778135e-005 

8.422716288e-007 

EarthS unLuna-JGM2-HP-SRP 3.46813818S 

0.00397058468 

EarthSun-0-0-0 

2.067991 746e-005 

8.179978462e-007 

Earth-0-0-0 

4.404163183e-005 

8.186544146-007 

Earth-0-0-SRP 

0.1637932094 

0.0001561216614 

Earth-O-HP-O 

3.213782043 

0.003648358687 

Earth- J GM2-0-0 

0.2076446785 

0.000240480153 


Table 2.37: FF/STK Molniya STK Test, Case Comparison 


Test Case 

Position Difference(m) 

Velocity Difference (m/s) 

AllPlaucts-0-0-0 

0.01457313187 

1 .2154368570-005 

EarthL una-0-0-0 

0.01409431455 

1.1843763756-005 

EaxthSuaLuna-JGM2-HP-SRP 

7.784661066 

0.004770671083 

EarthS nn-0- 0-0 

0.0006798546862 

1.0036567556-006 

Eaxth-0-0-0 

0.0003365057467 

8.08351 7746e-007 

Earth-0-0-SIlP 

2.16384821 

0.001779529358 

Earth-O-HP-O 

15.27354394 

0.01277094127 

Earth- JGM.2- 0-0 

3.063128789 

0.002561623786 



. 5 . FF/STK COMPARISON 


Table 2.3S: FF/STK SunSync STK Test Case Comparison 


Test Case 

Position Difference (in) 

Velocity Difference (in/s) 

AllPi&nets-O-O-O 

3.4087722-57e~005 

8.692021282&-007 

EarthLuna-0-0-0 

1 ,971678896e-005 

8 . 20060 841 7e-00 7 

EaxtliStmLuna-JG.M2-HP-SRP 

0.7800584002 

0.0008377097648 

EarthSmi-0-0-0 

1 .5382800 96 e- 00 5 

8.472033081. e-007 

Earth-0-0-0 

3.6561230330-005 

8.215166636C-007 

Earth-O-O-SRP 

0.2052485046 

0.0001828464608 

Earth-0- HP-0 

0.6524853805 

0.0007361489016 

Earth- JGM2-0-0 

0.09463135069 

0.0001055896859 
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Table 2.39: FF/STK EML2 STK Test. Case Comparison 

Test Case Position Difierence(m) 

Velocity Difference(m/s) 

AUPlaiicts-0-0-0 209675.042 

AlIPlanets-O-O-SRP 207842.4058 

EarthSunL una-0-0-0 209277.3283 

EarthS i mLivna- .7 GM 2-0- 0 209283.0998 

1.144026162 

1.126072542 

1.141861547 

1.141942278 


Table 2.40: 

FF/STK ESL2 STK Test Case Comparison 

Test Case 

Position Diifcrcncc(m) 

Velocity DifFercncc(m/s) 

A11P Ianets-0-0-0 

5661795.443 

.3.114044494 

AUPIanets-O-O-SRP 

6652873.129 

5.182879379 



Chapter 3 


Calculation Parameters 


GMAT’s Central Body (Cb) and Coordinate System (CS) dependent parameters were tested to verify that the 
internal calculations were correct. In order to minimize the effects of other forces/elements, the two-body cases 
from the Propagators section were used, with some modification, to test both the central body and coordinate 
system parameters. The only changes to the two body cases were in the report output intervals and report 
output parameters. Data was outputted in ten minute intervals. The ISS two-body case was used for the Earth 
case and each planets respective two-body case was used for the non-Earth cases. 


3.1 Initial Orbit State Conditions 


The ISS, GEO, Marsl, Mercuryl, Moon, Neptunel, Plutol, Satuml, Uranusl and Venusl two-body case’s initial 
orbit parameters were used from the Propagation section (Chapter 2) for the test cases in this section. 

Refer to Appendix B.l Tables B.l- B.13 for a listing of all Propagator initial orbit states used for the Calculation 
Parameter test cases. 


3.2 Central Body Dependent Parameters 

3.2.1 Naming Convention 

This section describes the naming convention for central body dependent parameter scripts and output reports. 
The naming convention consists of a case sensitive ordered series of option strings, separated by underscores (_ 
). Currently, options are allowed for the following fields, and will be present in the file name: 


1. tool - The tool used to generate the test case. 

2. traj - The trajectory to use. This includes initial conditions, physical parameters, and time step. 


CbParams precedes the tool field and 2Body follows the traj field. The central body used can be determined 
based on the traj field. The final Cb file format is as followed: 

CbParams_tooLfny'_2Body.report 
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The tool field should always be the first option field. Each field has a finite list of options, as follows (future 
options should be added to this list): 

1. tool 

STK - Satellite Toolkit HPOP or Astrogator 

FF - EreeFlyer 

GMAT - General Mission Analysis Tool 

2. traj 

ISS - leo orbit 

Mar si - eccentric low orbit 

Mercuryl - eccentric low orbit 

Moon - eccentric low orbit 

Pluto 1 - eccentric low orbit 

Venusl - eccentric low orbit 

NOTE: Some test cases contain traj variations. In this case traj precedes the modification. For example, 
if an ISS trajectory is needed with a different Cd, traj could be ISSdiffCdl. 


3.2.2 Comparison Script Information 


Comparison.Tooll-Tool2_Cb.m is the script used to perform the coordinate system comparisons needed for the 
Acceptance Test Plan. Many elements of this script were extracted from the ComparisomTooll.Tool2.CS.m 
script. 

Comparison.Tooll-Tool2_Cb.m was designed to allow the user to select two programs to compare to one another. 
The comparison involves taking the difference of the variables listed in the Acceptance Tast Plan Overview 
Chapter->Testing Methodology-> Calculation Parameters section. 

Refer to Appendix C for more details on this script and others used in the Acceptance Test Plan document. 


3.2.3 Test Case Results 


The following results are for the Central Body-Calculation Parameter section. The current GMAT Build is 
compared to STK and FreeFlyer for this section. 

FF-STK comparison results presented in Tables 3.11- 3.15 are used as a way to determine if the GMAT com- 
parison values are acceptable. If GMAT comparison data is within the same order of magnitude as the FF-STK 
comparison data, that is acceptable. A more detailed acceptance metric/matrix will be developed at a later date. 
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3.2. CENTRAL BODY DEPENDENT PARAMETERS 



Table 3.1: GMAT/STK Central Body Dependent Parameter Differences (1) 


Test Case 

Altitude (m) 

Eccentricity 

M. Anomaly (deg) 

M. Motion (rad/sec) 

Period' (sec) 

GEO-2Body 

1.834476006 

1.71688e-014 

1.90007e-005 

3.01146-015 

3.59432e-009 

Hyperbolic-2Bcdy 

0.00422986 

1.98064e-013 

351.401 

0.000223215 

N/A 

ISS-2Body 

0.00543367 

5.957196014 

5.81025&-010 

4.990156015 

2.54659e-010 

Marsl-2Body 

0.0144309 

9.74698e-011 

7.847916-007 

1.625116-013 

2.29212e-006 

Mercuryl-2Body 

0.00224375 

2.36836e-011 

1.72416-007 

5.54375e-0l4 

7.660826-007 

Moou-2Body 

0.00576446 

2.74547c-011 

6.471 92C-007 

4.03859(3-014 

8.68182e-007 

Ncptunol-2Body 

0.129938 

1.316846010 

6.04911c-007 

1.00481O-013 

4.02816(3-006 

Plutol-2Body 

0.19231 

1.838236009 

2.95607e-005 

3.827926-012 

0.000141823 

Saturnl-2Body 

0.752271 

1.053626011 

6.344 9e-008 

1.11997e-014 

1.046436006 

Uranusl-2Body 

0.216614 

5.31596e-011 

6.35362e-007 

8.48005e-014 

8.31598e-006 

Venusl-2Body 

0.000632367 

1.263236011 

2.1 3253e-008 

1.535916-014 

1.51652e-007 


Table 3.2: GMAT/STK Central Body Dependent Parameter Differences (2) 


Test Case 

Semi-rnajor Axis (m) 

True Anomaly (deg) 

Semilatus Rectum (m) 

GEO-2Body 

1.17871.e-006 

1.90007e-005 

7.42148e-007 

Hypcrb o!ic-2B ody 

1.72076(3-006 

7.41807(3-012 

1. 32204O-005 

ISS-2Body 

2.08274(3-007 

5.77472&-010 

1.70985(3-007 

Marsl-2Body 

0.000741843 

1. 157476-006 

0.000794765 

Mercnryl-2Body 

0.000199977 

2. 61636e-007 

0.000175894 

Moon-2Body 

0.000103879 

9.757226-007 

7.237176-005 

Neptimel-2Body 

0.00597339 

9.2047e-007 

0.00497125 

Plntol-2Body 

0.0111278 

4.499e-005 

0.0106515 

Satnrnl-2Body 

0.00241793 

9.55595e-008 

0.00235179 

Uranu3l-2Body 

0.0100126 

9.19017e-007 

0.00930848 

Veimsl-2Body 

0.000101745 

3.07957o-008 

0.00013864 


Ta ble 3.3: GM AT/S T K Central Body Dep endent Parameter Differences (3) 


Test Case 

Apoapsis Rad. (in) 

Periapsis Rad. (in) 

Apo. Vel. ( m/sec ) 

Per. Vel. (m/sec) 

GEO-2Body 

1.25146&-006 

1.27329e-006 

8.215656-011 

1.030296010 

Hyperbolic-2B(;dy 

N/A 

4.509276-006 

N/A 

1.868736009 

ISS-2Body 

5.26597&-007 

3.38332e-007 

4.396486010 

3.472786010 

Marsl-2Body 

0.000683647 

0.000800039 

2.75849e-007 

4.863096007 

Mercury 1 - 2B o dy 

0.000280259 

0.000125894 

7.849776008 

7.333166008 

Moon-2Body 

0.000193031 

1.47281 e-005 

7.360316008 

1.262326008 

Ncptuncl-2Body 

0.0107941 

0.00538058 

2.2861c-006 

2.187036006 

Plutol-2Body 

0.0145814 

0.00998796 

2.301846006 

3.550736006 

3aturnl-2Body 

0.002828 

0.00203526 

2.946486007 

5.082626007 

Uranusl-2Body 

0.0133233 

0.00759641 

1.412396006 

1.435346006 

Venusl-2Body 

0.000149412 

0.000183741 

6.662146008 

1.500226007 


f 
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Table 3.4: GMAT/STK Centoa! Body Dependent Parameter Differences (4) 


Test Case 

CS-Energy (mi/ sec*) 

Latitude (deg) 

Longitude (deg) 

MHA (deg). 

LST (deg) 

GE(>2Body 

2.6823e-007 

3.24498e-008 

3.14812e-007 

0.00196854 

0.00196864 

Hyp erb olic-2B o dy 

1.64846e-006 

5.65537e-008 

1.35016e-007 

0,00196849 

0,00196858 

ISS-2Body 

1.74794e-006 

1.07149e-007 

1.3276e-007 

0.00196849 

0.0019686 

Marsl-2Body 

0.00149955 

5.583386-007 

9.730Gle-007 

N/A 

N/A 

Meicuryl-2Body 

0.000332532 

1.88133e-007 

2.696‘51e-007 

N/A 

N/A 

Moou-2Body 

0.000140608 

7.42402^007 

1.06364e-0f)6 

N/A 

N/A 

N/A 

McptAmol-2Body 

0.0333365 

7.39515c~007 

1.15656O-006 

N/A 

Plutol-2Body 

0.00339014 

3.20069C-005 

7.44124C-005 

N/A 

N/A 

Saturnl-2Body 

0.014334 

1. 24281 e-006 

5.18835e-007 

N/A 

N/A 

TJranus.t-2Body 

0.028651 

9.64237e-007 

2.34538e-006 

N/A 

N/A 

Vkausl-2Body 

0.000500677 

2.22686&-008 

3.82501«-008 

N/A 

N/A 


Table 3.5: GMAT/STK Central Body Dependent Parameter Differences (5) 


Test Case 

Beta Angie (deg) 

(RxV)-Mag (in? /sec) 

R-Mag (in) 

GEO-2Body 

0.000773444 

0.00180444 

1.88447e-006 

Hyp erb olic-TBody 

0.000706521 

0.025655 

1.35624e-005 

ISS-2Body 

0.00344889 

0.000749424 

4.62933e-007 

Marsl~2Body 

0.00315966 

1.23719 

0.0126011 

Mercury l-2Body 

0.00511923 

0.220796 

0.00224375 

Moon-2Body 

0.00394712 

0.000777618 

0.00576446 

Neptunel-2Body 

0.0007 41045 

35.4559 

0,0763421 

Plutol-2Body 

0.000530341 

4.01956 

0.19231 

Safcurnl-2Body 

0.00142058 

20.1241 

0.0181717 

Uraimsl-2Body 

0.000807228 

53.902 

0.101524 

Venusl-2Body 

0.00414838 

0.447268 

0.000632365 
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Table 3-6: FF/GMAT Central Body Dependent Parameter Differences (1) 


Test Case 

Altitude (m) 

Eccentricity 

M. Anomaly (deg) 

M. Motion (rad/sec) 

Period (sec) 

GEO-2Bodv 

8.20728c-006 

1.93962O-010 

6.421 53O-006 

5.69653c-016 

7.75617o-009 

ISS-2Body 

0.433346 

6.39335e-012 

1.7l4e-009 

4.2466e-0l5 

7.69433e-010 


Table 3.7: FF/GMAT Central Body Dependent Parameter Differences (2) 


Test Case 

Semi-major Axis (m) 

True Anomaly (deg) 

! Semilatus Rectum (m) 

GEO-2Body 

2.233726-006 

6.4215e-006 

2.4811e-006 

ISS-2Body 

4.09273e-007 

1.867086-009 

1.7917e-007 


Table 3.8: FF/GMAT Central Body Dependent Parameter Differences (3) 


Test Case 

Apoapsis Rad. (in) 

Pcriapsis Rad. (in) 

Apo. Vcl. (m/.s«c) 

Per. Vcl. (m/sec) 

GEQ-2Body 

2.72848e-006 

3.00497e-006 

2.39764e-008 

2.39013e-008 

ISS-2Body 

3.98359e-007 

6.893976-007 

2.918386-008 

2.88605e-008 


Tabic 3. 9: FF/GM AT Centr a l Body Dependent Paramet er Differences (4) 


Test Case 

C3-Energy (m*/sec?) 

Latitude (deg) 

Longitude (deg) 

MILA (deg) 

LST (deg) 

GEO-2Body 

0.00014728 

5.40781e-008 

0 

0.00300898 

0.00300919 

ISS-2Body 

0.000444231 

6.02545e-006 

5.81869e-007 

0.00300898 

0.00300919 


Test Case 

Beta Angle (deg) 

(RxV)-Mag- (m 2 /sec) 

R.-Mag (m) 

GEO-2Body 

45.7413 

0.00429281 

3.274186-006 

ISS-2Body 

98.3649 

0.00164437 

1.18325c-006 
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Table 3.11: F'F/S'TK Central Body Dependent Parameter Differences (1) 


Test Case 

Altitude (m) 

Eccentricity 

M. Anomaly (deg) 

M. Motion (rad/sec) 

Period (sec) 

GEO-2Body S.04721c~006 

1.9396O-010 

1.77052c-005 

2.44679c-015 

6.4756c-009 

ISS-2Body 

0.427915 

6.34e-012 

1.73014e-009 

9.18666e-015 

8.14907e-010 


Table 3.12: FF/STK Central Body Dependent Parameter Differences (2) 


Test Case 

Semi-major Axis (m) 

True Anomaly (deg) 

Scmilatus Rectum (m) 

GEO-2Body 

2.08092e-006 

1. 77044 e-005 

1.92085e-006 

ISS-2Body 

4.47471e-007 

1.67336e-009 

1.55524e-007 


Table 3.13: FF/STK Central Bod}" Dependent Parameter Differences (3) 


Test Case 

Apoapsis Rad. (m) 

Peri apsis Rad. (m) 

Apo. Vel. (m/stc) 

Per. Vel. (m/sec) 

GEO-2Body 

2. 74304c- 006 

2.03727o-006 

2.38942(3-008 

2.38796c-008 

ISS-2Body 

5.23869C-007 

6.24823c- 007 

2.87441C-008 

2.88605C-008 


Table 3.14: FF/STK Central Body Dependent Parameter Differences (4) 


Test Case 

03-Energy [v&jsecA) 

Latitude (deg) 

Longitude (deg) 

MHA (deg) 

LST (deg) 

GEO-2Body 

0.000147141 

3.06S91e-008 

0 

0.00497708 

(1.00497715 

XSS-2Body 

0.000442895 

5.92635e-006 

6. 1431e-007 

0.00497708 

0.00497716 


Tabic 3.15: FF/STK Central Body Dependent Parameter Differences (5) 


Test C&sc 

Beta Angle (deg) 

(RxV)-Mag (m z /sec) 

R-Mag (m) 

GEO-2Body 

0.000773444 

0.00355067 

4.05998e-006 

IS8-2Body 

0.00344889 

0.00178261 

1.01772c-006 


3.3 Coordinate System Dependent Parameters 
3.3.1 Naming Convention 

This section describes the naming convention for coordinate system dependent parameter scripts and output 
reports. The naming convention consists of a case sensitive ordered series of option strings, separated by under- 
scores (_ ). Currently, options are allowed for the following fields, and will be present in the file name: 


1. tool - The tool used to generate the test case 

2. traj - The trajectory to use. This includes initial conditions, physical parameters, and time step 

3. CS - The coordinate system to use. The celestial body to use is followed by the CS in the name 

CSParams precedes the tool field and 2Body precedes the CS field. The final CS file format is as followed: 
CSParams_tooLir<y _2Body_C5.report 
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The tool field should always be the first option field. Each field has a finite list of options, as follows (future 
options should be added to this list): 


1. tool 


STK 

Satellite Toolkit HPOP or Astro] 

FF 

FreeFlyer 

GMAT - 

General Mission Analysis Tool 

traj 

ISS 

- leo orbit 

S unSync 

- leo orbit 

GPS 

- meo orbit 

GEO 

- geo orbit 

Molniya 

- heo orbit 

Marsl 

- eccentric low orbit 

Mercuryl 

- eccentric low orbit 

Moon 

- eccentric low orbit 

Plutol 

- eccentric low orbit 

Venusl 

- eccentric low orbit 


NOTE: Some test cases contain traj variations. In this case traj precedes the modification. For example, if 
ISS trajectory is needed with no output, then traj can be ISSnoOut. The lack of a report file is shortened 
to noOut. 


3. CS 

EarthFixed 

EarthTODEq 

EarthMODEc 

MarsFixed 

MercuryFixed 

MoonFixed 

NeptuneFixed 

PlutoFixed 

SaturnF faced 

UranusFixed 

VenusFixed 


EarthMJ2000Eq 

EarthTODEc 

EarthGSM 

MarsMJ2000Eq 

MercuryMJ2000Eq 

MoonMJ2000Eq 

NeptuneMJ2000Eq 

PlutoMJ2000Eq 

SaturnMJ2000Eq 

UranusMJ2000Eq 

VenusMJ2000Eq 


EarthMJ2000Ec 

EarthMODEq 

EarthGSE 

MarsMJ2000Ec 

MercuryMJ2000Ec 

MoonMJ2000Ec 

NeptuneMJ2000Ec 

PlutoM J 2000Ec 

SaturnMJ2000Ec 

UranusMJ 2000Ec 

VenusM J2000Ec 


3.3.2 Comparison Script Information 

The script used to perform the Coordinate System comparisons needed for the Acceptance Test Plan is Compar- 
ison_Tooll_Tool2-CS.m. Many elements of this script were extracted from the ComparisondToolldXbol2_PV.nl 
script. 

Comparison_Tooll_Tool2_CS.m was designed to allow the user to select two programs to compare to one another. 
The comparison involves taking the difference of the variables listed in the Acceptance Test Plan Overview 
Chapter-> Testing Methodology-> Calculation Parameters section. 

Refer to Appendix C for more details of this script and others used in the Acceptance Test Plan document. 


3.3.3 Test Case Results 


The following results are for the Coordinate System-Calculation Parameter section. The current GMAT Build 
is compared to STK and FreeFlyer for this section. 






CHAPTER 3. " CALCULATION PARAMETERS 


FF-STK comparison results presented in Tables 3.26- 3.30 are used as a way to determine if the GMAT com- 
parison values are acceptable. If GMAT comparison data is within the same order of magnitude as the FF-STK 
comparison data that is acceptable. A more detailed acceptance metric/matrix will be developed at a later date. 
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3.3. COORDINATE SYSTEM DEPENDENT PARAMETERS 


Table 3.16: GMAT/STK Coordinate System Dependent Parameter Differences (Position) 


Test Case 

GEO-2Body-EarthFixcd 

GEO-2Body-EarthMJ2000Ec 

GEO-2Body-EarfchMJ2000Eq 

GEO-2Body-EarthMODEc 

GEO-2Body-EarthMODEq 

GEO-2Body-EarthMOEEc 

GEO-2Body-EartbMOEEq 

GEO-2Body-EartbTODEc 

GEO-2Body-EarthTODEq 

G EO - 2 B o dy-E art hT OEEe 

GEO-2Body-EarthTOEEq 

Hyperbolic-2Budy-EarthFixed 

Hyperbolic-2Body-EarthMJ2000Ec 

Hyp erb olic- 2B ody-EartbM J 2000Eq 

Hyp erbolio- 2B o dy- E artbMODEc 

Hyq>erbolic.-2Body-EartbMODEq 

Hypcrbolic-2Body-EartbMOEEc 

Hyp erbolio- 2B o dy-Ear tliMOEEq 

Hyp erbolic-2B o dy-EarthTODEc 

Hyperbolic- 2Body-EartiTODEq 

Hyp erb olio 2B o dy-E ar fchT OEEe 

Hyp erb olic- 2B o dy- E arthT OEEe; 

ISS-2Body-EarthFixed 

ISS-2Body-EarthGSE 

ISS-2Body-EarthGSM 

ISS-2Body-EarthMJ2000Ec 

ISS-2Body-EarthMJ2000Eq 

ISS-2Body-EarthM0DEc 

ISS-2Body-EarthM0DEq 

ISS-2Body-EartbMOEEc 

lSS-2Body-EarthMOEEq 

I S S - 2B o dy- E artliTO D Ec 

ISS-2Body-EarthTODEq 

ISS-2Body-EartbTOEEc 

ISS-2Body-EarthT0EEq 

Marsl-2Body-MarsFixed 

Marsl-2Body-MarsMJ2000Ec 

Marfil-2Body-MarsMJ2000Eq 

Mercury l-2Body-MerciiryFixed 

Mercuryl-2Body-MercuryMJ2000Ec 

Mercury l-2Body-MercuryMJ 2000Eq 

Moon-2Body-MoonFixed 

Moon-2Body-MoonMJ2000Ec. 

Mo on- 2Body-Moo nM J 2000Eq 
N eptunel- 2Body-NeptuneFixed 
Neptunel- 2Body-NeptuneMJ2000Ec 
Nepti me 1 - 2B o dy- Nepti meM J20 OOEq 
P!utol-2Body-Pluto Fixed 
Plutol-2Body-PhitoMJ2000Ec 
Plutol-2Body-PlutoMJ2000Eq 
S atu r o. 1 - 2B o dy- S atumFixed 
Satunil-2Body-SaturxiMJ2000Ec 
S aturnl - 2B o dy- S aturnM J2000Eq 
Ur anus 1- 2B o dy- U ranusFixed 
Urannsl-2Body-UranusMJ2000Ec 
Uranusl-2Body-UranusMJ2000Eq 
Venus 1 - 2Bo dy- VenusFixed 
Venus 1 - 2Body- VenusM J2000Ec 


X-Pos (m) 
0.2281413654 
2.2940980(i7e-005 
2.29409S067e-005 
0.0012728946 79 
0.001272895247 
2.295053037^005 
2.295053037e-005 
0.006357750408 
0.006357749726 
0.005487542296 
0.005487542239 
0.3488040966 
1.6822014e-005 
1. 682201 4e-005 
0.008121191058 
0.008121191058 
1.693842933O-005 
1 .6938429330-005 
0.01367428012 
0.01367420191 • 
0.03536941949 
0.03536941949 
0.01006907041 
5.3005067e-005 
4.9540460626-005 
1.033299668O-005 
1.0332996686-005 
0.000210427288 
0.0002104279702 
1.0331405066-005 
1.033140506e-005 
0.0003583738817 
0.0003583786565 
0.0009193927326 
0.0009193925052 
0.05898887696 
0.04854124924 
0.04854124947 
0.009607810114 
0.008947269691 
0.008947269691 
0.02659000859 
0.02282095541 
0.02282095546 
0.4543228424 
0.2946186363 
0.2946186363 
0.6182226091 
0.7527405821 
0.7527405819 
0.6340523651 
0.07298453147 
0.07298453147 
0.4395935848 
0.385162361 
0.3851023601 
0.003286658284 
0.002135182285 


Y-Pos (m) 
'0.04023080692 
0.00032580283 
2.433989721e-005 
0.001382489927 
0.001271371502 
2.2371750670-005 
2.433989721e-005 
0.006834058695 
0.006271528946 
0.005975813906 
0.005486127634 
0.5436684769 
0.001917796908 
1.979060471 e-006 
0.01077665365 
0.0098834862 
2.211891115C-006 
1.979060471 e-006 
0.01960279769 
0.02924879664 
0.04713787348 
0.04508509301 
0.01444909003 
3.604054655&-005 
0.003799987553 
8.047527444c-005 
9.6G7473932e-006 
0.0001799348865 
0.0001652065293 
1.099556357e-005 
9.6665644380-006 
0.0002511692401 
0.0003260520316 
0.0007893422662 
0.0006983391359 
0.05711927406 
0.06378650704 
0.06874534654 
0.01088151507 
0.01142762343 
0.01220055815 
0.01816027718 
0.02955069925 
0.03242772568 
0.474102093 
0.379445413 
0.4144292316 
1.111364997 
0.9494719067 
1.085654527 
0.6134889563 
0.09-547611444 
0.1048656477 
0.5642743054 
0.5115039512 
0.5548048448 
0.002153384223 
0.003261576126 


Z-Pos (in) 
0.02037045935 
0.0007504095265 
4.46(1-009 
8.563802112e-006 
0.000543207701 
9.603354556e-006 
4.162467864e-009 
0.0006267073331 
0.003423117541 
9.603354556e-006 
0.002404046552 
0.1095402695 
0.004851564881 
2.1245 79623e-006 
7.392372936&-006 
0.004294328392 
2.153683454c- 006 
2.066371962o-006 
0.003635126632 
0.005363341188 
2.1536834540-006 
0.01691197394 
0.002770384526 
7.67059305e-006 
0.00293166886 
0.0001261635134 
9.215000318e-006 
7.613095931e-006 
7.1365775510-005 
7.695803106e-006 
9.215000318e-006 
0.0001079586127 
0.0003840964382 
7.6960 3048c- 006 
0.0002934716576 
0.03425882358 
0.05515928102 
0.04854127501 
0.009587843124 
0.01022404834 
0.008947248972 
0.0259134155 
0.02583996015 
0.02282108797 
0.2913599029 
0.3350099975 
0.30129434 
0.7629325749 
0.8475838174 
0.7521445426 
0.06666205172 
0.08217663634 
0.0737594446 
0.5180439475 
0.433244951 
0.3789037312 
0.002565443879 
0.002492622912 
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Table 3.17: GMAT/STK Coordinate System Dependent Parameter Differences (Velocity) 


Test Case 


X-Vel (m/s) 


Y-Vel (m/s) 


Z-Vel (m/s) 


GE 0-2Body--Ear fcbPixod 

GEO-2Body-EarthMJ2000Ec 

GEO-2Body-EarthMJ2GQ0Eq 

GE0-2Body-EarthM0DEc 

GEO-2Body-EarthMODEq 

GEO-2Body-EarthMOEEc 

GE0-2Body-EartliM0EEq 

GEO-2Body-EarthTODEc 

GEO-2Body-EarthT ODEq 

GEO-2Body-EarthTOEEe 

GEO-2Body-EarthTOEEq 

Hyperbolic-2Budy-EarthPix«l 

Hyperbolic-2Body-EarthMJ2000Ec 

Hyperbolio2Body-EarthMJ2000Eq 

Hyperb olic-2B ody-EarthMODEc 

Hyperbobc-2Body-EartliMODEq 

Hypcrbolic-2Body-EarthMOEEc 

Hypcrbolio-2Bady-EarthMOEEq 

Hyperbolio-2Body-EarthTODEc 

HyperboliC'2Body-EarthTODEq 

Hyperbolio2Body-EarthTOEEc 

Hyperbo.lio2Body-EarthTOEEq 

IS S- 2B ody-E artliFixed 

ISS-2Body-EartliGSE 

ISS-2Body-EartbGSM 

ISS-2Body-EarthMJ2000Ec 

ISS- 2B ody-E ar t hM J 2000Eq 

ISS-2Body-EarthMODEe 

ISS-2Body-EartbM0DEq 

IS S- 2B o dy-E artbMOEEe 

ISS-2Body-EarthMOEEq 

LSS-2Body-EarthT0DEc 

ISS-2Body-EarthTODEq 

IS S- 2B 0 dy-EarthT 0 EEc 

IS S-2Body-EarthT OEEq 

Mars 1-2B ody-M arsF ixed 

Marsl-2Body-MarsMJ2000Ee 

Marsl-2Body-Mai\sMJ2000Eq 

Mercuryl-2Body-MercuryFixed 

Mercurjd.-2Body-MercaryMJ2000Ec 

Mercuryl-2B ody-MercuryM J2000Eq 

Moon-2Body-MooiiFixed 

Moon-2Body-MoonMJ2000Ee 

Moon-2Body-MootiMJ2000Eq 

N eptimel- 2B 0 dy-N eptuueF ixed 

Neptaael-2Body-NepfcuneMJ2000Ec 

Neptun.el-2Body-NeptiineMJ2000Eq 

Pluto l-2Body-PlutoFixed 

Plutol-2Bo dy-P lutoM J2000Ec 

Plutol-2Body-PlutoMJ2000Eq 

8atixml-2Body-SaturuFixed 

Saturnl-2Body-SatiiriiMJ2Q00Ec 

8aturnl-2Body-SaturnMJ2000Eq 

TJranusI-2Body-UraniisFixed 

Uranusl-2Body-UramwMJ2000Ec 

UranusI-2Body-UramisMJ2000Eq 

Venusl-2Body-VenusFixed 

Vmms l-2Bod v- VetiusM J2000Ec 


1 .01 4408074e-005 
1.65 8007065©- 009 
1 . 658007065e-009 
9.276047985e-008 
9.2760479S5©-008 
1 . 658895243e-009 
1.6588952436-009 
4.572472123e-007 
4.572472123C-007 
4.000988846e-007 
4.0009888466-007 
3.9610398526-005 
1.829647545e-010 
1.900701818e-010 
1.2022205456-007 
1 . 202220545e-00 7 
1.829647545O-010 
1.8296475456-010 
1.646727199e-007 
1 . 6467271996-00 7 
•5.242357659e-007 
5.242357659e-007 
1. 12218190 3e-Q05 
3.753886091e-008 
3.523581427e-008 
1 .08441 86660-OO8 
1.0844186666-008 
2.404640931e-007 
2.404640931e-007 
1. 08423825-5e-008 
1.0842382556-008 
4.0732328626-007 
4.073232862e-007 
1.051030152c-006 
1.0510301526-006 
4.119460101e-G05 
3.9224898e-005 
3.9224898e-005 
8.297573739e-006 
7.406872116e-006 
7.4068721 1 6e-006 
2.01 5.969569e-005 
1.554381 4016-005 
1. 554381 401e-Q05 
0.002311382259 
0.0001425149572 
0.0001425149572 
0.0003613682175 
0.0003872987759 
0.0003872987759 
0.0001164046752 
2. 48890774 7c-005 
2.488907747e-005 
0.0001044264559 
0.0001229694768 
0.00012296947G8 
3.2984535836-006 
2.213821571e-0 06 


1.788716823o006 
2.383426789e-008 
1.708189146e-009 
1.009415607e-007 
9.28549633&-008 
1.566302643e-009 
1.708S5528e-009 
5.0514091526-007 
4.636103862O-007 
4.3591429776-007 
4.002142819e-007 
4.825957234e-005 
2.841105129e-008 
7.416289804e-011 
3.212767836-007 
2.9465979666-007 
8.48210390 Sc-011 
7.283063042C-01.1 
2.454381143e-007 
3.61354946e-007 
1.403959171e-006 
1.284048596-006 
1.427654617©-005 
0.989109796-008 
4.26608171e-006 
9.146905455O-008 
1.1376122276-008 
2.0532398e-007 
1.883493361e-007 
1.291988 738e-008 
1.137490102e-008 
2. S78270955e-007 
3.764446532e-007 
8.982550259e-007 
7.945391012c-007 
5.0401936856-005 
5.6123893096-005 
5.5061041e-005 
1.056043l31e-005 
1.0845782056-005 
1.0469359826-005 
1.201479294e-005 
2.227603246c-005 
2.2340361516-005 
0.00232459858 
0.0002082204933 
0.0002067147653 
0.0004667307785 
0.0005522331561 
0.0005410307473 
0.0001207579849 
3.488300981O-005 
3.502235213e-005 
0.0001481343457 
0.000170892658 
0.0001702024259 
1.926899973e-006 
3. 1 0787484Ce-006 


2 . 60 3041 31 8c- 00 7 
5.4726889696-008 
0 

6.031009026e-01 0 
3.9615844036-008 
6.8270389346-010 
3.035338586-013 
4.6019854596-008 
2.603101342O-007 
6.827038934e-010 
1.75307769e-007 
7.70518005e-006 
7.1950001516-008 
7.5051076466-011 
2.7339241986-010 
1.2804208766-007 
5.306866058o-011 
7.416289804c-011 
6.046185774e-008 
2.50673704e-007 
5.306866058e-011 
5.6251658896-007 
3.0566100766-006 
8.1583628746-009 
3.338248789e-006 
1.441 839981 c-007 
1.0165868156-008 
8.086198378e-009 
8.1784801156-008 
8.1923356996-009 
1.0165368556-008 
1.225990426-007 
4.359921 313e-007 
8.192335699O-009 
3.3386 89325o-007 
2.244996011e-005 
3.7734461116-005 
3.923428094e-005 
6.283098936-006 
6.91.1486272e-006 
7.379157396e-00C 
1.474075317e-005 
1.46090765O-005 
1.5543840026-005 
0.0001281027884 
0.0001364492723 
0.0001464608648 
0.0004682860251 
0.000344448214 
0.0003868007483 
2.3541264056-005 
2.32026131 66-005 
2.510042485e-005 
0.000160607621 
0.0001139574768 
0.000120316952 
2.241729913©- 006 
2.148337952e-006 
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Table 3,18: GMAT/STK Coordinate System Dependent Parameter Differences (Specific Angular Momentum) 


Test Case 


X-(H) ( m 2 /sec ) Y-(H) hrr/sec) Z-(H) {m?/stc) 


GEO-2Body-EarthFixcd 

10.77395299 

2.103745961 

434.3113329 

GEO-2Body-EarthMJ2000Ec 

1.13710075a- 005 

2.306616807 

1.000254997 

GEO-2Body-EarthMJ2000Eq 

0 

0 

0.001542503014 

GEO-2Body-Ead;b.MODE<: 

1.670208945 

0.001840817276 

0.001789885573 

GEO-2Body-EaithMODEq 

1.670203311 

0.001656591829 

0.002357410267 

GEO-2Body-EarthMOEEc 

2.159488468e-005 

Q. 000 8076312952 

0.001688022166 

GEO-2Body-EaxtliMOEEq ' 

1. 27982535 7e-005 

2.844632424e-021 

0.001542503014 

GE 0-2Body-EarthT ODEc 

8.5322038 

1.939137292 

0.8412171155 

GE 0-2Body-EarthT ODEq 

8.532204845 

9.575566487 

0.004365574569 

GEO-2Body-EarthTOEEc 

7.285380082 

0.001295120455 

0.001702574082 

GEO~2Body-Eai:tliTOEEq 

7.285369859 

1.252006488 

0-001396983862 

HyperboUe-2Body-EarthFbced 

12999.11505 

5061.392323 

11584.82954 

Hypeibolic-2Body-EartliMJ2000Ec 

0.0002509729892 

1.80931238 

0.7227790775 

Hyperbolic-2Body-EarthMJ2000Eq 

0.000197488248 

0.01781154424 

0.01760781743 

HypeiboBc-2Body-EartbMODEc 

1.183702615 

0.007530616131 

0.02339947969 

Hyperbolic-2Body-EarthMODEq 

1.183709749 

0.01674925443 

0.01817534212 

Hyp crb olic- 2B ody-EarthM OEEc 

0.00026047443 

0.008927599993 

0.0233121682 

Hypcrl)olic-2Body-EaxthMOEEq 

0.000217340812 

0.01784064807 

0.01766602509 

Hyperbolic-2Body-Earth.TODEe 

2.066649429 

1.515742042 ' 

0.6022310117 

Hyperbolie-2Body-EaithTODEq 

2.066892929 

5.087713362 

5.08732046 

Hypeibolic-2Body-EarthTOEEe 

5.182947816 

0.00926957 

0.02366141416 

Hyperbolic-2Body-EarthTOEEq 

5.182958375 

0.6900227163 

0.6991467671 

ISS-2Body-EarthFixed 

89.91428695 

67.10288835 

31.77834151 

ISS-2Body-EartliGSE 

0.2475171641 

0.0290456228 

0.0004511093721 

ISS-2Body-Ea.rthGSM 

3.261639904 

22.63645729 

29.37958561 

ISS-2Body-EartLMJ2000Ec 

0.0007530616131 

0.7919024938 

0.2579763532 

ISS-2Body-EarthMJ2000Eq 

0.0007603375707 

0.0005566107575 

0.000301952241 

ISS-2Body-EarthMODEc 

0.4284956958 

0.9411305655 

0.001266016625 

ISS-2Body-EarthMODEq 

0.4284956958 

0.8629867807 

0.3752284101 

ISS-2Body-EarthMOEEc 

0.0007530616131 

0.0002983142622 

0.0004874891602 

ISS-2Body-EarthMOEEq 

0.0007003375707 

0.0005566107575 

0.000301952241 

ISS-2Body-EarthTOD Ee 

0.7472881407 

2.260314432 

0.2173910616 

ISS-2Body-EarthTODEq 

0.7472772268 

3.630950232 

2.533914085 

ISS-2Body-EarthTOEEc 

1.879539923 

4.109922884 

0.0004874891602 

lSS-2Body-EarthTOEEq 

1.879554475 

4.081139195 

1.909433195 

Marsl-2Body-MarsFixed 

23.98886318 

22.21717 

10.38143819 

Marsl-2Body-MarsMJ2000Ec 

0.9952709661 

0.5993915693 

0.5886722647 

Marsl-2Body-MarsMJ200GEq 

0.9952691471 

0.2261380806 

0.7740909496 

Mercuryl-2Body-MercuryFixed 

0.327698217 

0.4808271115 

0.2502738425 

Mercury 1 -2Body-Mercur yM J 2000Ec 

0.1878624971 

0.2510159902 

0.2059750841 

Mercury l-2Body-Mercm-yMJ2000Eq 

0.1878634066 

0.1111429479 

0.2888264135 

Moon-2Body-MoonFixed 

5.71125554 

5.283694918 

0.1482353582 

Moon-2Body-MoonMJ2000Ec 

0.001037733455 

0.04335993253 

0.01902935765 

Moon-2Body-MoonMJ2000Eq 

0.00.1038642949 

0.0002576208312 

0.0003556124284 

Nept\mel-2Body-NeptiineFixed 

42849.13634 

43033.72977 

96195.52572 

Nep fcuneI-2Body-N eptuneM J20(K)Ec 

29.59038829 

77.10719365 

79.3626532 

Nepfcvmel-2Body-Neptnae.MJ2000Eq 

29.59038829 

84.00197472 

57.43693328 

Pluto l-2Bo dy-PlutoFbced 

12.66916559 

15.03451797 

18.39295624 

Plutol-2Body-PlutoMJ2000Ec 

2.923440004 

0.8473472235 

3.204913128 

Pluto 1- 2Bo dy-P lutoM J2000Eq 

2.923445891 

0.8302175236 

3.147571306 

Saturnl-2Body-SafcuruFixed 

7341.511315 

7984.99-5318 

884.71 7701 9 

Satunil-2Body-SatumMJ2000Ec 

13.2243149 

25.11031926 

27.86074765 

Saturnl-2Body-Sa.turnMJ2000Eq 

13.22408207 

16.91115992 

24.03883263 

Uranus l-2Body-UramisFixed 

2955.545322 

3574.169823 

2353-711105 

Uranusl-2Body-UranusMJ2000Ec 

5.186127964 

39.40385068 

90.66710481 ■ 

Uraniisl-2Body-UranusMJ2000Eq 

5.186419003 

61.77691968 

77.83877663 

Venus l-2Bo dy- VenusF ixed 

0.7762544101 

0.3817622201 

0.2173546818 

Venusl-2Body-VenusMJ2000Ec 

0.2975648385 

0.6312984624 

0.4800604074 
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Table 3.19: GMAT/STK Coordinate Sj^stem Dependent Parameter Differences (Velocity Vector-based) 


Test Case 

GEO-2Body-EaxthFixcd 

GEO-2Body-EarthMJ2000Ee 

GEO-2Body-Eartb.MJ2()OOEq 

GEO-2Body-EarthMODEc 

GE 0- 2Bo dy-E arthM O DEq 

GE 0- 2B o dy-Ear thMO EEc 

GEO-2Body-EartliMOEEq 

GEO- 2Body-EarthTODEc 

GEO-2Body-EarthTODEq 

GEO-2Body-EarthTOEEc 

GEO-2Body-EarthTOEEq 

Hyp erb uhc- 2Bo dy-E ar thFixed 

Hyperbolic-2Body-Earth.MJ2000Ec 

Hyperbohc-2Body-EarthM32000Eq 

Hyperbolic-2Body-EartbJVIODEc 

Hyperbolic- 2Body-EarthMODEq 

Hypcrbolic-2Body-EartbMOEEc 

Hypcpbolic-2Body-EartidVfOEEq 

Hyperbolic- 2Body-EartkTODEe 

Hyperbolic-2Body-Earth.TODEq 

Hyp erbolic- 2Body-EartbT0EEc 

Hyperbolic-2Body-Eaa:tbT0EEq 

ISS-2Body-EarthFixed 

ISS-2Bcdy-EaxtliGSE 

IS S - 2B o dy-E arth GSM 

IS S- 2B o dy- E arthM J2000Ec 

IS S- 2B o dy-E arthM J 2000Eq 

ISS-2Body-EarthMODEc 

I.SS-2Body-EarthMODEq 

ISS-2Body-EarthMOEEc 

IS S- 2B o dy-E arthMOEEq 

IS S- 2B o dy-E artkT 0 D Ec 

I S S - 2B o dy- E arthT 0 D Eq 

ISS-2Body-EaithTOEEc 

ISS-2Body-EarthTOEEq 

M arsl-2Body-MarsFixed 

Marsl-2Body-MarsMJ2000Ec 

M ar s 1 - 2B o dy-MaxsM. J 2000Eq 

Mercury l-2Body-MercuryF.bc.ed 

Mercury l-2Body-MercuryMJ2000Ec 

Mercury l-2Body-MercuryMJ2000Eq 

M o on-2Bo dy-Mo onF ixed 

Moon-2Body-MoonMJ2000Ec 

Mo on-2Body-Mo onM J 2000Eq 

N ep tunel-2Body-NeptuneFixed 

Neptunei-2Body-NeptuneMJ2000Ec 

Neptunel-2Body-NeptuneMJ2000Eq 

Plutol-2Body-PIutoFixed 

P Into 1- 2B o dy-PlutoM ,J2000Ec 

Plutol-2Body-PlutoMJ2000Eq 

S at urn 1 -2Body- S at urnFixed 

S aturnl - 2B o dy- S aturnM J2000Ec 

Satuml-2Body-SatumMJ2000Eq 

Ur am is 1 -2B o dy- Ur anu sFixed 

UrannKl-2Body-UranusM,T2000Ec 

Uranusl-2Bcdy-UraimsMJ2000£q 

Venusl-2Body-VenusFixed 

Venusl.-2Body-VenusMJ2000Ec 


Mag-Vel (m/s) 
3.2785831 87o-006 
1.07913678&-010 
1.07913678e-010 
1.070254996e-010 
1 .0 70254996e-01 0 
1.07913678e-010 
1.07.91.3678e-010 
1,070254996(3010 
1.070254996O-010 
1.07913678e-010 
1.07913678e-010 
4.71472088e-005 
1 . 6 1 6484 724e-0 1 0 
1.616484724e-010 
1.616484724e-010 
1.616484724e-010 
1.616484724C-010 
1.61 6484724e-0 1 0 
1.616484724e-010 
1.616484724e-010 
1.616484724e-010 
1.6164847246-010 
1.0002674246-006 
4.574118861e-010 
3.1252334006-008 
4. 6362 91 351 c-010 
4.636291351e-010 
4.52970994e-010 
4..627409567»-010 
4.636291351e-010 
4.627409567e-010 
4.627409567e-010 
4.6274095(57e-010 
4.627409567c-010 
4.627409567(3010 
9.5928069936-006 
8.899687209e-006 
8.89967 6995(5-006 
1.59941 127e-006 
1.597606047e-006 
1 .597 60G047e-006 
3.394075021e-006 
3.384815983c-006 
3.384815983e-006 
0.001546170623 
3.205422061(5-005 
3 . 205422 949e-005 
8.2178436176-005 
8 . 345 2 92 61 3e-00 5 
8.3452926136-005 
8.7160199066-006 
5.184041 64c-006 
5 . 1 84041 64e-006 
3.448082886e-005 
2.627271023(5-005 
2.6272710236-005 
5.125109226e-007 
5.1260951046-007 


Right Asc. of Vel. (deg) Dec. of Vel. (deg) 


1.162789617 
2.742908123(5-010 
3.426237072(5-011 
1.8879120496-009 
1.733170052e-009 
3.G28031209e-011 
3.427658157(5-011 
9.472358897e-009 
8.702897958^009 
8.130868423e-009 
7.460869256(5-009 
2.061432909e-007 
2.155502443(5-010 
1.421085472e-012 
1.854090215e-009 
1.910734682e-009 
1. 591615728O-012 
1.44950 71 Slo-012 
3.401112281e-009 
5.0393C8034e-009 
8.096151305(5-009 
8.594184919e-009 
1.289334008(5-007 
5.30000932e-010 
2.998996962e-008 
8.2936857386-010 
1.526245796e-010 
1.852903608e-009 
2.479552563e-009 
1.305977548e-010 
1.526245796(3-010 
3.180247177(5-009 
4.23S856377e-009 
8.094446002C-009 
1.066266009o-008 
9.947378707(5-007 
1.432667602e-006 
1.215923845e-006 
3.112265077(3-007 
3.378602003e-007 
2.85895851e-007 
1.115449194e-006 
1.24027315c-006 
1.048296383e-006 
1.309179918(5-005 
1.140866715(5-006 
9. 779481 331 e-007 
4.784315172e-005 
5. G52545477e-005 
4.844967677e-005 
4.920609911e-007 
1.242314625O-007 
1.061025898(3-007 
7. 15 7807744 e- 006 
1 .159199428(5-006 
9.843604971e-007 
3.693500616&-008 
3.7 34896836e- 008 


0.178531639 

1.110862513e-009 

0 

1.108624303e-011 
7.381565394e-010 
1.261213356e-011 
■5. 656303434(3-01 5 
9. 347012053e-01 0 
4. 8507531 33C-009 
1.261213356e-011 
3.266776493e-009 
5.501368605e-008 
7.097700205(5-010 
9.201528428e-013 
1.578293052e-012 
7.375620115e-010 
2. 41584530 2c-013 
■9.166001291(3-013 
5.87760951(1-010 
1.719470788&-009 
2.415845302(5-013 
3.265999027e-009 
3 . 8 71 IS 74 59e-00 8 
6.073808123(5-011 
3.987214825e-008 
1. 106641 889C-009 
7.80993048e-011 
6.026290578o011 
7.349498787e-01() 
6.103206829e-011 
7.80993048e-011 
9.443787974e-010 
4.168732914(3009 
6.103206829O-011 
3.25S224002c-009 
3.948958787e-007 
6.136789708e-007 
6.201466789e-007 
1.2746111016-007 
1.397476694(3007 
1.42502526-007 
5.325181203e-007 
5.264534106(3007 
5.364797317e-007 
1.351642366e-005 
4. 847768089(3007 
4.998895804e-007 
2.956015004e-005 
2.377500014e-005 
2.5127144856-005 
9. 227543885 e-008 
5.306250017C-008 
5.502680978e-008 
7.584381549e-007 
5.016063795e-007 
5.082584709e-007 
1.783747905e-008 
1.708282227e-008 
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Table 3.20: GMAT/STK Coordinate System Dependent Parameter Differences (Angle-based) 

Test Case Arg. of Per. (deg) Decl. (deg) Inc. (deg) PA (deg) R A AN (deg) 

GEO-2Body-Eartb.Fi>:cd 

3.357c-005 

2.76SC-008 

0.2044 

3.148c-007 

9.529c-005 

GEO-2Body-EaitbMJ2000Ec 

1.349e-005 

l.llle-009 

l.llle-009 

2.715e-010 

0 

GEO-2Body-EarthMJ2000Eq 

1.309e-005 

0 

le-01.4 

3.513e-011 

0 

GE 0-2Bod y-E arthM ODEc 

1.381e-005 

1.148e-011 

3.233e-013 

1.8866-009 

1 . 8646-009 

GEO-2Body-EarthMODEq 

1.334e-005 

7.381e-010 

7.548e-010 

1.734ft-009 

8.657e-010 

GEO-2Bodj'-EarthMOEEc 

1.349e-005 

1.291e-011 

3.126e-013 

3.401e-011 

0 

GEG-2Body-Ear(;liMOEEq 

1.343e-005 

5-C56e-015 

0 

3.513(5-011 

0 

GEO-2Body-EarfchTODEc 

1.356e-005 

9.2Sle-010 

9.348e-0.10 

9.56-009 

9.4826-009 

GEO-2Bodv-EarthT ODEq 

1.095c-005 

4.652c-009 

3.537O-009 

8.705(5-009 

1.046C-0Q5 

GEO-2Body-EarthTOEEc 

1.3Ge-O05 

1.291e-011 

3.233e-013 

8.129e-009 

8.107(5-009 

GEO-2Body-EaxthTOEEq 

3.3.1 4e-005 

3.267e-009 

1.764e-009 

7.462e-0Q9 

6.978e-005 

Hyperbolic-2Bod,y-EarthF ixed 

3.504e-007 

1.754e-008 

2.357e-007 

1.356-007 

1.5286-007 

Hyp erb olic-2B odv-EarthM J2000E-C 

7.617e-012 

1.096e-009 

l.llle-009 

2.172(5-010 

1.4(5-013 

Hyperboljc-2Body-EarthMJ2000Eq 

7.844e-0l2 

8.207e-013 

4.334e-013 

1.251(5-012 

2.56-013 

Hyperbolic-2Body-EartliMODEc 

1.18e-011 

1.561e-0i2 

3.411e-013 

1.854e-009 

1.858e-009 

Hj'perbolic-2Body-EartbMODEq 

1.052e-009 

7.379e-010 

5.89Se-013 

2.405e-009 

9.7o7e-010 

HypcrboEc-2Body-EarthMOEEc 

7.518C-012 

2.984c-013 

3.411C-013 

1.393(5-012 

1.4c-013 

Hypcrfco]ic-2Body-EartbMOEEq 

7.674C-012 

8.313e-013 

1. 421o-01 3 

1.251O-012 

0 

Hy perb olic-2B ody -EarthT 0 DEc 

9.422e-012 

9.216e-010 

9.348e-010 

3.396O-009 

3.2356-009 

Hypexbolic-2Body-EarthTODEq 

1.824e-009 

3.924e-009 

4.129e-009 

4.887e-009 

1.696e-009 

Hyperbolic-2Boriy-EarthTOEEc 

7.702e-012 

2.984e-013 

3.375e-013 

8.096e-009 

8.107(5-009 

Hyperbolic-2Body-Ea.rthTOEEq 

4.5dle-009 

3.265e-009 

5.531.6-010 

1.062e-008 

4.29e-009 

ISS-2Body-EartliFixed 

3.20 6e- 008 

3.275e-00S 

3.965e-008 

1.32S&-007 

1.38e-007 

ISS-2Body-EartliGSE 

7.G14e-010 

G.504e-011 

5.684e-013 

5.681e-010 

4.463e-010 

IS S- 2 B o dy-E arth GS M 

2.427e-006 

3.92e-008 

4.043e-00S 

2.862e-008 

7.1326-009 

ISS-2Body-EarthMJ2000Ec 

2.188c-009 

1.109C-009 

4.623e-010 

8.22c-010 

1.29o-009 

ISS-2Body-EarthMJ2000Eq 

5.556e-010 

7.966e-011 

5.471 e-013 

1.525e-010 

5.4e-013 

ISS-2Body-EarthMODEc 

5.451e-010 

6.455e-011 

1.592e-012 

1.853e-009 

1.852e-009 

ISS-2Body-EarthMODEq 

1.206e-009 

7.378e-010 

5.283e-010 

2.496-009 

1.291(5-009 

ISS-2Body-EartliMOEEc 

5.49e-010 

6.525e-011 

4.121(^013 

1.2696-010 

7.1056-013 

ISS-2Body-EartiMOEEq 

5.5256-010 

7.966e-011 

5.471e-013 

1.525e-010 

5.542e-013 

ISS-2Bod.y-EarthTO.DEc 

1.664e-009 

9.461e-010 

3.891e-0l0 

3.1786-009 

4.239&-009 

ISS-2Body-EarthTODEq 

3.32e-009 

4.163e-009 

3.569o-009 

4.2486-00.9 

4.395e-009 

ISS-2Body-EarthTOEEc 

5.487c-010 

6. 5260-011 

4.05c-013 

8.095c-009 

8.095C-009 

ISS-2Body-EaiihTOEEq 

2.92C-009 

3.264C-009 

2.689c- 009 

1.063c-008 

5.96c-009 

M ars 1- 2 B o dy- Mar sFixed 

2.247e-006 

5.221e-007 

5.014e-008 

9.731e-007 

2.604e-007 

Maisl-2Body-MarsMJ2000Ec 

4.431e-008 

8.568(^007 

2.152e-009 

1.285(5-006 

1.562e-009 

M ais 1-2B odv-M arsM J2000Eq 

4.543e-008 

7.205e-007 

2.049e-009 

1.377(5-006 

1.332(5-009 

MerOTiyl-2Body-Mercurj'Fixed 

7. 2C3O-009 

1.881O-007 

1.462cs-009 

2.697(5-007 

5.0666-009 

Mercurjd.-2Body-Mercury.MJ2000.Ec 

8.609e-009 

1.99e-007 

1.0046-009 

3.065(5-007 

2.2376-009 

Mercuryl-2Body-MercuryMJ2000Eq 

7.477e-009 

1.735&-007 

1.6646-009 

3.307&-007 

1.026e-009 

Moon-2Body-MoonFixed 

2.966e-008 

7.424e-007 

2.672e-009 

1.064e-006 

1.4186-007 

Moon-2Body-MoonMJ2000Ec 

4.765c-009 

7.387C-007 

4.114C-010 

1.173C-006 

8.903C-010 

Moon-2Body-MoonMJ2000Eq 

6.126e-009 

6.356(5-007 

1.155e-011 

1.246e-00C 

6.068(5-012 

Neptunel-2Body-NeptuneFixed 

4.002e-005 

5.911e-007 

1.341e-005 

1.157e-006 

7.627e-006 

N ep tune 1- 2 Body-N ep fcuneM J 200 OEc 

2.001e-008 

6.812e-007 

1.011e-008 

1.058e-006 

.1 ,0426-008 

Neptunel-2Body-NeptimeMJ2000Eq 

2.142e-008 

6.033(5-007 

6.696-009 

1.147(5-006 

1.42(5-008 

Plutoi-2Body-PlutoFixed 

5.512e-006 

3.201e-005 

7.484 e-007 

7.441&-005 

8.1176-007 

Plutol-2Body-PlutoMJ2000Ec 

5.338e-007 

3.396-005 

4.962e-008 

5,l05e-005 

3.472e-008 

Plutoi-2Body-PlutoMJ2000Eq 

5.077e-007 

2.963e-005 

3.89Se-008 

5.8786-005 

5.172e-00S 

SaturD.l-2Body-SaturnFix-.ed 

1.508e-007 

5.809e-008 

5.751e-008 

5.1886-007 

4.571 e-007 

Satuml-2Body-SatuiuMJ2000Ec 

3.259C-009 

7.328C-008 

6.576c-010 

1.179&-007 

9.333c-010 

8 aturn 1-2 B o dy-SatumM J20 OOEq 

4.55e-009 

6.548e-Q08 

5.3836-010 

1.273e-007 

8.028e-010 

TJranusl-2Body-UranusFbced 

4.529e-006 

8.4e-007 

3.03e-007 

2.3456-006 

4.615e-007 

Uranusl-2Body-UrarmsMJ2000Ec 

1.129e-008 

ti.9e-007 

8.766(5-009 

1.112(5-006 

5.248e-009 

Uranusl-2Body-UramisMJ2000Eq 

1.392e-008 

6.011e-007 

6.433e-009 

1.201(5-006 

1 e-008 

Venusl-2Body~VenusFixed 

1.012e-009 

2.227e-008 

1.379e-010 

3.825e-008 

1.158e-009 

Venusl-2Body-VeuusMJ2000Ec 

2.027e-009 

2.162e-008 

4.20.9(5-010 

3.81 66-008 

s 7nQo.ni n 



58 


CHAPTER ^". "'"'CALCULATION PARAMETERS 
Table 3.21: FF/GMAT Coordinate System Dependent Parameter Differences (Position) 


Test Case 

X-Pos (m) 

Y-Pos (m) 

Z-Pos (m) 

ISS-2Body-EarthMJ2000Eq 

2.54757424e-005 

2.516884479c-005 

2.321598913e-005 


Table 3.22: FF/GMAT Coordinate System Dependent Parameter Differences (Velocity) 

Test Case X-Vel (m/a) Y-Vel (m/s) 

Z-Vel (m/s) 

ISS-2Body-EarthMJ2000Eq 4.993130354e-007 5.00476105&-007 

4.961786537e-007 


Table 3.23: FF/GMAT Coordinate System Dependent Parameter Differences (Specific Angular Momentum) 
Test Case X-(H) (ru'V'sec) Y-(H) (: irP/sec, ) Z-(H) ( m'-jsec ) 

lSS-2Body-EarthM,T2000Eq 0.0005493347999 0.001349690137 0.001567968866 


T able 3.24: FF/GMAT Coordinate System Dependent Parameter Differences (Velocity Vector-base d) 
Test Case Mag-Vel (m/s) Right Asc. of Vel. (deg) Dec. of Vel. (deg) 

ISS-2Body-EarthMJ2000Eq. ■ 4860Q67904el0G7 ■ : "29612928006 ' 67.36750845 


Table 3.25: FF/GMAT Coordinate System Dependent Parameter Differences (Angle-based) 

Test Case Arg. of Per. (deg) Dccl. (deg) Inc. (deg) RA (deg) RAAN (deg) 

ISS-2Bodv-EarthMJ2000Eq 1.681e-009 5.8S4e-010 4.684e-010 7.628e-010 2.418&-010 
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COORDINATE SYSTEM DEPENDENT PARAMETERS 



Table 3.26: FF /STK Coordinate System Dependent Parameter Differenc es (Position) 
Test Case X-Pos (mj Y-Pos (m) Z-Pos (m) 

.ISS-2Body-EajthMJ200QEq 3.580873909m005 3.470836418e-005 3.243G98945e-005 


Tab le 3.2 7: F.F /S TK Coordi nate System Dependent Parameter Differences (Velocity) 
Test Case X-Vel (m/s) Y-Vel (m/s) Z-Vel (m/s) 

ISS-2Body-Earth.MJ2000Eq 5.092002375e-007 5.105600387e-007 5.025699S56e-007 


Table 3.2 8: FF/STK Coordinate System Dependent Parameter Differences (Specific Angular Momentum) 
"Test Case X-(H) (in? /sec) Y-(H) (mf/se c) Z-(H) {nP/se^ 

ISS-2Body-EarthMJ2000Eq 0.000571102672? 0.001346052159 0.00146246748 


Table 3.29: FF/S TK C oordinate System Dependent Parameter Differences (Velocity Vector-based) 

Test Case _ __ Mag- Vel (m/a ) Right Asc. of Vel. (deg) Dec. of Vel. (deg) 

IS S-2 B o dy-’E ar fc'h M J 2 OOOEq 4.857803049e-007 296.2928006 67.36750845' 


Table 3.30: FF/STK Coordinate System Dependent Parameter Differences (Angle-based) 

Test Case Arg. of Per, (deg) Peel. "(dcgT Inc, (deg) RA (dog) R A AN (dog) 

ISS-2Body-EarthMJ2000Eq 1.343e-009 C.518e-010 4.685e-010 8.637e-010 2.4146-010’ 
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CHAPTER 3. "'CALCULATION PARAMETERS 



Chapter 4 





Integrators 


GMAT’s integrators were tested on a system level in order to verify that all the integrators were working correctly. 
In order to minimize the effects of other forces/elements, the two-body cases from the Propagators section were 
used, with some modification, for testing the integrators. The report output interval and integrators were the 
only parameter changed when using the Propagator two body test cases. Data was outputted in ten minute 
intervals. 


4.1 Initial Orbit Conditions 


The ISS and GEO two-body case’s initial orbit parameters were used from the Propagation section (Chapter 2) 
for the test cases in this section. 

Refer to Appendix B.l Tables B.l- B.13 for a listing of all Propagation initial orbit states used for the Integrator 
test cases. 


4.2 Naming Convention 

This section describes the naming convention for integrator scripts and output reports. The naming convention 
consists of an ordered series of option strings, separated by underscores (_ ). Currently, options are allowed for 
the following fields, and will be present in the file name in order: 

1. tool - The tool used to generate the test case 

2. traj - The trajectory to use. This includes initial conditions, physical parameters, and time step 

3. integ - The integrator to use 


The word Integrator precedes the tool field and 2Body follows the integ field. The final integrator file format is 
as followed: 

Integrator_fooLtraj'Jnteg_2Body.report 

The tool field should always be the first option field. Future additional fields should be added to the end of the 
list of fields. Each field has a finite list of options, as follows (future options should be added to this list): 


6.1 
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CHAPTER 4. INTEGRATORS 


STK 

1. tool FF 

GMAT 


- Satellite Toolkit HPOP or Astrogator 

- FreeFlyer 

- General Mission Analysis Tool 


. ISS - leo orbit 

t-rO,J Q£Q _ g e0 or bit 

NOTE: Some test cases contain traj variations. In this case traj precedes the modification. For example, 
if an ISS trajectory is needed with a different Cd, traj could be ISSdiffCdl. 


RKV89 

RKN68 

RKF56 

3. integ PD45 
PD78 
BS 
ABM 


- RungaKutta 8(9) 

- DormandElMikkawyPrinee 6(8) 

- RungeKuttaFehlberg 5(6) 

- PrinceDormand 4(5) 

- PrinceDormand 7(8) 

- BulirschStoer 

- AdamsBashforthMoulton 


4.3 Comparison Script Information 

ComparisonJnteg.m is the script used to perform the integrator comparisons needed for the Integrators section 
of the Acceptance Test Plan. This script was designed 'to allow the user to select a GMAT Build or the exact' 
analytic solution to compare to one another. 

Refer to Appendix C for more details about this script and others used in the Acceptance Test Plan document. 


4.4 Test Case Results 


The following results are for the Integrator section. The GMAT Integrator results are being compared to an 
exact analytical two-body solution. We’d like the comparison data to be as close to zero as possible. A detailed 
acceptance metric/matrix will be developed at a later date. 


Table 4.1: Exact /GMAT GEO Integrator Test Case Comparison 
Test Case Position Differeuce(m) Velocity Differencefm/s) 

ABM-2Body 0.0001127712184 8.17S70?643e-009 

BS-2Body 6.290317204e-005 4.5351896966-009 

PD45-2Bodv 3.297321562e-005 2.3209081 48e-009 

PI)7S-2Body 1.82237034 6e-005 1.253046654e-009 

RKF56-2Body 3.459697578c-005 2.564071918cv009 

RKN68-2Body 1.793126562c-005 I.29779348Qe-009 


RKV89-2Body 0.0001127712184 8,178707643^009 



4.4. TEST CASE RESULTS 


Table 4.2: Exact/GMAT ISS Integrator Tfest Case Comparison 
Test Case Position Difference(m) Velocity Difference(m/s) 


ABM-2Body 

BS-2Body 

PD45-2Body 

PD7S-2Body 

RKF56-2Body 

RKN68-2Body 

RKVS9-2Body 


1 .015221 1 1 23o-005 
1.691079091e-005 
3.080434654e-005 
2.772621462(8-005 
3.336029534&-005 
3.6170251256-006 
1 .052211123(3-005 


1.172665432e-008 

1.926879767e-008 

3.521816658e-008 

3.15224572e-008 

3.80U9779e-008 

3.875411333e-009 

1.1 72065432(3-008 
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Chapter 5 


Stopping Conditions 


GMAT’s Stopping Conditions were tested on a system level in order to determine if it stops correctly on user 
selected conditions. Refer to Table 5.1 for a list of the stopping conditions tested and the units used for the 
output of data. 


Table 5.1: Stopping Conditions 


Stopping Condition 

Stopping Value 

Unit Used 

Epoch (A1 Modified Julian Date) 

23158.042037434974 

Days 

Apo apsis 

180 

TA in degrees 

Elapsed Days 

6 Hours 

Days 

Elapsed Days (Multiple Satellites) 

5 Hours 

Days 

Mean Anomaly 

45, 90, & 180 

MA in degrees 

Periapsis 

0 

TA in degrees 

Elapsed Seconds 

3600 

Seconds 

True Anomaly 

45, 90, & ISO 

TA in degrees 

XY Plane Intersection 

z=o 

Km 

XZ Plane Intersection 

Y— 0 

Km 

YZ Plaue Intersection 

x=o 

Km 
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CHAPTER 5. STOPPING CONDITIONS 


5.1 Initial Orbit Conditions 


For a listing of the initial conditions used to produce the data for the Stopping Conditions Section, refer to 
Table 5.2 for the Earth based non-hyperbolic point mass test cases, Table 5.3 for the Earth based non-hyperbolic 
all forces test cases, and Table 5.4 for the Earth based hyperbolic point mass test cases. 


Table 5.2: Initial Orbit Parameters (EarthPM k EarthM J 2000EqPM) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System (CS’j 

Earth Mean J2000 Equator 

X 

-8043.9600382.977915 (km) 

Y 

-1564.9950345568864 (km) 

Z 

3750.9601677510364 (km) 

VX 

0.99861303787927636 (km) 

VY 

-6.8834168529193462 (km) 

vz 

-0.46566090709653452 (km) 

Mass (No Fuel) 

850 (kg) 

Cd 

2.2 

Cr 

1.8 

Drag Area 

15 (m 2 ) 

Drag Model 

None 

NSG Model 

None 

.SRP Area 

1 (m 2 ) 

SRP 

Off 

Integrator Type 

RuugaKutta 8(9) 

Integrator Init. StepSize 

60 (sec) 

Integrator Accuracy 

le-13 

Integrator Max. StepSize 

2700 (sec) 

Report Precision 

16 significant figures 

Report StepSize 

Only initial state 

Report CS/Cb 

Same as initial state CS 
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5.1. INITIAL ORBIT CONDITIONS 

Table 5.3: Initial Orbit Parameters (EajthAil fc EartliMJ2Q00EqA]l) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System (CS) 

Earth Mean J2000 Equator 

X 

-8043.9600382977915 (km) 

Y 

-1564.9950345568864 (km) 

Z 

3750.9601677510364 (km) 

vx 

0.99861303787927636 (km) 

VY 

- 6 . 8834168529193462 (km) 

vz 

-0.46566090709653452 (km) 

Mass (No Fuel) 

850 (kg) 

Cd 

2.2 

Cr 

1.8 

Drag Area 

15 (w 2 ) 

Drag Model 

Jacehia Roberts 

Drag F107/F107A 

150/150 

Drag Kp 

3 

NSG Model 

JGM2 

Degree x Order 

4x4 

SRP Area 

1 ( m 2 ) 

SRP 

On 

Integrator Type 

RungaKutta 8(9) 

Integrator Init. StepSize 

60 (sec) 

Integrator Accuracy 

le-13 

Integrator Max. StepSize 

2700 (sec) 

Report Precision 

16 significant figures 

Report StepSize 

Only initial state 

Report CS/Ob 

Same as initial state GS 



CHAPTER 5. STOPPING CONDITIONS 
Table 5.4: Initial Orbit:- Parameters (EarthPMhyper & EaxthMJ2000EqPMhyper) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System (CS) 

Earth Mean J2000 Equator 

X 

12371.791482634855 (km) 

Y 

5050.7627227610719 (km) 

Z 

5050.762722761071 (km) 

vx 

-7. 9859921 51 2608487 (km) 

VY 

2.44520073255755 (Ion) 

VZ 

2.4452007325575495 (km) 

Mass (No Fuel) 

850 (kg) 

Cd 

2.2 

Cr 

1.8 

Drag Area 

15 (m 2 ) 

Drag Model 

None 

MSG Model 

None 

SRP Area 

1 (m 2 ) 

SEP 

Off 

Integrator Type 

RungaKutta 8(9) 

Integrator Init. StepSizc 

60 (sec) 

Integrator Accuracy 

lo-13 

Integrator Max. StepSize 

2700 (sec) 

Report Precision 

16 significant figures 

Report StepSize 

Only initial state 

Report CS/Cb 

Same as initial state CS 


Table 5.5: Initial Orbit Parameters (MoonPM) 

Initial State Parameter 

Parameter Value (unit) 

Coordinate System (CS) 

Earth Mean J2000 Equator 

X 

-1486.792117191545200 (inn.) 

Y 

0.0 (km) 

Z 

1486.792117191543000 (km) 

VX 

-0.142927729144255 (kin) 

VY 

-1 .631407624437537 (km) 

VZ 

0.142927729144255 (km) 

Mass (No Fuel) 

850 (kg) 

Cd 

2.2 

Cr 

1.8 

Drag Area 

15 (m 2 ) 

Drag Model 

None 

NSG Model 

None 

SPiP Area 

1 (to 2 ) 

SRP 

Off 

Integrator Type 

RungaKutta 8(9) 

Integrator Init. StepSize 

60 (sec) 

Integrator Accuracy 

le-13 

Integrator Max. StepSize 

2700 (sec) 

Report Precision 

1 6 significant figures 

Report StepSize 

Only initial state 

Report CS/Cb 

Same as initial state CS 




.1. INITIAL ORBIT CONDITIONS 


Table 5.6: Initial Orbit Parameters (MarsPM) 


Initial State Parameter Parameter Value (unit) 


Coordinate System (OS) 

Earth Mean J2000 Equator 

X 

-2737.481646173082000 (kin) 

Y 

0.0 (km) 

Z 

2737.4816461 78082000 (km) 

vx 

-0.311321695052649 (km) 

VY 

-3.553492313930950 (km) 

vz 

0.311321695052650 (km) 

Mass (No Fuel) 

850 (kg) 

Cd 

2.2 

Cr 

1.8 

Drag Area 

15 (to 2 ) 

Drag Model 

None 

NSG Model 

None 

SEP Area 

1 (m 2 ) 

SEP 

Off 

Integrator Type 

R.ungaKutta 8(9) 

Integrator Init. StepSize 

60 (sec) 

Integrator Accuracy 

le-13 

Integrator Max. StepSize 

2700 (sec) 

Ecport. Precision 

16 significant figures 

Report StepSize 

Only initial state 

Report CS/Cb 

Same as initial state C3 
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5.2 Naming Convention 


CHAPTER 5. STOPPING CONDITIONS 


This section describes the naming convention for stopping condition scripts and output reports generated for 
use in GMAT’s Acceptance Test Plan. The naming convention consists of an ordered series of option strings, 
separated by underscores (_ ). Currently, options are allowed for the following fields, and will be present in the 
file name in order: 


1. tool - The tool used to generate the test cases. 

2. Cb - The Central Body used for the stopping condition, including the force model setup. 

3. stopCond - The stopping condition used for the test case. 

The word StopCond precedes the tool field. The final stopping condition file format is as followed: 
StopCond- tool -Cb. stop Cond.y&poxi 


The tool field should always he the first option field. Future additional fields should be added to the end of the 
list of fields. Each field has a finite list of options, as follows (future options should be added to this list): 


1. tool 


STK 

FF 

GMAT 


- Satellite Toolkit HPOP or Astrogator 

- FreeFlyer 

- General Mission Analysis Tool 


EarthPM 

EarthMJ2000EqPM 

EarthAll 

2. Cb EarthMJ2000EqAll 
Earth^MultiSatsPM 
MarsPM 
MoonPM 


- Earth central body with point mass force model 

- Same as EarthPM and plane intersection calculations based on EarthMJ2000Eq 

- Earth central body with all force model types turned on 

- Same as EarthAll and plane intersection calculations based on EarthMJ2000Eq 

- Same as EarthPM and the test case involves # satellites 

- Mars central body with point mass force model 

- Moon central body with point mass force model 


3. stopCond 


AlMod Julian 

Apoapsis 

Days 

MA### 

Periapsis 

Seconds 

TA### 

XY plane 
XZplane 
YZplane 


- Al Modified Julian Date 

- Apoapsis of orbit based on TA 

- Elapsed Days 

- ### degree Mean Anomaly 

- Periapsis of orbit based on TA 

- Elapsed Seconds 

- ### degree True Anomaly 

- XY Plane Intersection 

- XZ Plane Intersection 

- YZ Plane Intersection 


5.3 Comparison Script Information 


Comparison-StopCond.m is the script used to perform the comparisons needed for the Stopping Condition chap- 
ter of the Acceptance Test Plan. This script was designed to allow the user to select a GMAT Build and/or the 
exact solution to compare to one another. 


Refer to Appendix C for more details on this script and others used in the Acceptance Test Plan document. 



5.4. TEST CASE RESULTS 

5.4 Test Case Results 
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The following results are for the Stopping Condition section. The GMAT Stopping Condition results are being 
compared to the exact desired stopping condition. We’d like the comparison data to be as close to zero as 
possible. A detailed acceptance metric/matrix will be developed at a later date. 

Refer to Table 5.1 to determine the units used in the Difference(s) column of the Stopping Condition comparison 
results. 


Table 5.7: Exact/GMAT EarthPM StopCond Test Case Comparison 


Stopping Condition Difference (s) 


AlMod.Tulian 

3.637978807e-012 

Apoapsis 

1. 20741 8308e-006 

Days 

3.6379788(J7e-012 

EA180 

1 . 33485 1 106e-006 

EA45 

1 . 1 2245S571e-009 

EA90 

5.233573575c-010 

MA18CI 

1 ,468336194e-006 

MA90 

5.793253877e-009 

Peri apsis 

0 

Secs 

2.095480 34e-007 

TA180 

1.20741 8308e-006 

TAQO 

2.69331224e-009 


Tabic 5.8: Exact/GMAT Earth All StopCond Test Case Comparison 


Stopping Condition Difference(s) 


AlModJuliau 

7.275957614e-012 

Apoapsis 

0 

Days 

7.27601 3125C-012 

EA180 

0 

EA45 

1.3701 33873e-008 

EA90 

1.320385934e-008 

MAI 80 

0 

MA90 

1.408977823e-008 

Peri apsis 

0 

Secs 

4.190951oS6e-007 

TA180 

0 

TA90 

8.0765971 69e-010 
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Table 5.9: Exact /GM AT EartbPMhyper StopCond Test Case Comparison 


Stopping Condition 

Difference(s) 

A 1 Mo cl Julian. 

3.637978807e-012 

Days 

0 

HA45 

5.94486238e-009 

HA90 

6.062620628O-009 

MA45 

4.0322731 82e-009 

MA90 

7.014335779e-010 

Peri apsi s 

0 

Secs 

2.09548034e-007 

TA45 

1.530786164e-008 

TA90 

3.041570551e-008 


Table 5.10: Bxact/CMAT MoonPM StopCond Tost Cas e Comparison 


Stopping Condition Difference^) 


AlModJidian. 

3.6379788076-0.12 

Apoapsis 

1.2074183086-006 

Days 

3.637978807c-012 

EA180 

0 

EA90 

3.1936338026-009 

MA180 

1.774535207e-006 

MA90 

8.9537905986-009 

Periapsis 

0 

Secs 

2.095480346-007 

TA180 

1.2074183086-006 

TA90 

9.5943448746-009 


Table 5.11: Bxact/GMAT MarsPM StopCond Test Ca.s e Comparison 


Stopping Condition Difference^) 


AlModJidian 

3.637978807e-012 

Apoapsis 

0 

Days 

0 

EA180 

0 

EA90 

2.2204986296-008 

MAI 80 

0 

MA90 

1 .276124806e-008 

Periapsis 

0 

Secs 

2.095480346-007 

TA180 

0 

TA90 

1.5057395336-008 


Table 5.12: Exa.ct/ 'GMAT EarthMJ2000EqPM SuipCu.ad Tes t Case Comparison 
Stopping Condition Difference (s) 

XYplane l.!)67496619e-007 

XZplauc 1. 187821 727c-00C 

YZplane 4.2192789e-007 



. 4 . TEST CASE RESULTS 


Table 5.13: Exact/GMAT EartliMJ2000EqAll StopCond Test Case Comparison 


Stopping Condition 

Difference(s) 

XYplane 

8.57494129e-008 

XZplane 

1.575498128e-006 

YZplane 

2.401038834e-008 


Table 5.14: Exact/GMAT EarthMJ2000EqPMliyper StopCond Test Case Comparison 


Stopping Condition 

Diiference(s) 

XYplane 

9.138221522e-007 

XZplane 

9.138204291O-007 

YZplane 

3. 63783191 4e-006 


Table 5.15: Exact/GMAT MoonMJ2000EqPM StopCond Test Case Comparison 


Stopping Condition 

Difference (s) 

XYplane 

1.752778189e-007 

XZplane 

0 

YZplane 

1.752923708e-007 


Table 5.1b: Exact/GMAT Mai:sMJ2000EgPM StopCond Test Case Comparison 


Stopping Condition 

Difference^) 

XYplane 

2.3841 85791c-007 

XZplane 

0 

YZplane 

2.2351 74179e-007 


Table 5.17: Exact/GMAT Earth3MultiSatsPM StopCond Test Case Comparison 


Stopping Condition 

Difference(s) 

Days 

7.276e-012 7.276e-012 7.276e-012 
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Chapter 6 


Libration Points 


The libration point tests are designed to verify the location of various libration points. In each software tool 
used the location of the libration point is defined by the tool and then converted to a MJ2000Eq representation 
for comparison purposes. 


6.1 Initial Orbit Conditions 


Table 6.1: Initial Orbit Parameters (Sun-Earth (SE) Libration Points) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System fCS) 

Sun-Earth L#* MJ2000Eq 

X 

0 (k m) 

Y 

0 (km) 

Z 

0 (km) 

VX 

0 (km) 

VY 

0 (km) 

VZ 

0 (km) 

Mass (No Fuel) 

850 (kg) 

Cd 

2.2 

Or 

1.8 

Drag Area 

15 (m 2 ) 

Drag Model 

None 

NSG Model 

None 

SRP Area 

1 (m 2 ) 

SRP 

Off 

Integrator Type 

N/A 

Integrator Init. StepSize 

N/A 

Integrator Accuracy 

N/A 

Integrator Max. StepSize 

N/A 

Report Precision 

16 significant figures 

Report StepSize 

Only initial state 

Report CS/Cb 

Earth MJ2000Eq 


NOTES: (*) All five libration points are defined in the test script. 
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CHAPTER 6. LIBRATION POINTS 


Table 6.2: Initial Orbit Parameters (Sim-Earth-Mooii(SEM) Libration Points) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System (CS) 
X 

Sun-Earth-Moon L#* MJ2000Eq 

0 (km) 

Y 

0 (km) 

Z 

0 (1cm) 

vx 

0 (1cm) 

VY 

0 (1cm) 

VZ 

0 (kra) 

Mass (No Fuel) 

850 (kg) 

Cd 

2.2 

Cr 

1.8 

Drag Area 

15 (m 2 ) 

Drag Model 

None 

NSG Model 

None 

SllP Area 

1 (m 2 ) 

SRP 

Off 

Integrator Type 
Integrator Init. StcpSizc 

N/A 

N/A 

Integr ator A ccuracy 

N/A 

Integrator Max. StepSize 

N/A 

Report Precision 

16 significant figures 

Report StepSize 

Only initial state 

Report CS/Cb 

Earth MJ2000E<j 


NOTES: (*) All five libration points are defined in tile test script. 


6.2 Naming Convention 


This section describes the naming convention for libration point scripts and output reports generated for use in 
GMAT’s Acceptance Test Plan. The naming convention consists of an ordered series of option strings, separated 
by underscores (_ ). Currently, options are allowed for the following fields, and will he present in the file name 

in order: 

1. tool - The tool used to generate the trajectory. 

2. libType - The type of libration point used for the test case. 

3. UbPoint - The libration point used for the test case. 


The word LibrationTest precedes the tool field. The final stopping condition file format is as followed: 
Libr&tio-nTest-tooLlibTypeJibPoint.report 


The tool field should always be the first option field. Future additional fields should be added to the end of the 
list of fields. Each field has a finite fist of options, as follows (future options should be added to this list): 

STK - Satellite Toolkit HPOP or Astrogator 

FF - FreeFlyer 

00 GMAT - General Mission Analysis Tool 

OD - Orbital Determination Toolbox 

„ ... _ SEML - Sun Earth Moon 

2. hbType S£L _ SuQ Earth 

3. UbPoint # - The libration Point number 
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COMPARISON SCRIPT INFORMATION 



6.3 Comparison Script Information 



Comparison_Tooll_Tool2Xibr.m is the script used to perform the comparisons needed for the Libration Points 
chapter of the Acceptance Test Plan. This script was designed to allow the user to select two tools from a list 
for comparison. The Tools available are presented in this chapter’s Naming Convention section. 


Refer to Appendix C for more details on this script and others used in the Acceptance Test Plan document. 


6.4 Test Case Results 


The following results are for the Libration Points section. The current GMAT Build is compared to STK for the 
Libration Points section. We’d like the comparison data to be as close to zero as possible. A detailed acceptance 
metric/matrix will be developed at a later date. 


Table 6.3: GMA.T/STK SEM LibrationTest Test Case Comparison 


Libration Point 

Position Difference(m) 

Velocity Difference(m/s) 

1 

2.997694537e-006 

1.082467449e-011 

2 

7.217749953e-006 

-2.8969S8205e-011 

3 

8.94069671 6e-005 

1.548983164e-009 

4 

-0.0006556510925 

18.76805413 

5 

0.0005066394806 

-1.37023761 


Tab ic 6.4: GMAT /STK SE LibrationTest Test Case Comparison 


Libration Point 

Position Difference (m) 

Velocity Difference (m/s ) 

1 

2.997694537e-006 

1.082467449e-011 

2 

7.2l7749953e-006 

-2.896988205e-011 

3 

8.94069671 6e-005 

1.5489831 84 c-009 

4 

-0.0006556510925 

18.76805413 

5 

0.0005066394806 

-1.37023761 
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LIBRAT10N POINTS 



Chapter 7 



Delta V 


Once an initial state was created for these Delta V test cases, a set amount of Delta V was applied. For the 
impulsive bum test cases delta V values of 0.1, 0.1, and 0.1 were applied in the X,Y, and Z axes respectively for 
either the Cartesian or VNB axes. Each test case generate a report that contains the Cartesian elements of the 
state before and after the impulsive burn. 


7.1 Initial Orbit Conditions 

Initial conditions for the impulsive burn delta V test cases axe presented in Tables 7.1- 7.4. 
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Table 7.1: Initial Orbit. Parameters (Earth Impulsive Burns) 


Initial State Parameter 

Parameter Value (unit) 

Start & Stop Time 

01 Jan 2000 11:09:28.000 (UTCG) 

Central Body 

Earth 

Coordinate System 

Earth Mean J2000 Equator 

X 

7378.0 (km) 

Y 

0.0 (km) 

Z 

0.0 (km) 

vx 

0.0 (km) 

VY 

0.1973811193846027 (km) 

V'Z 

0.1973811193846018 (km) 

Mass (No Fuel) 

850 (kg) 

Cd 

2.2 

Cr 

1.8 

Drag 

15 (m 2 ) 

Drag Model 

None 

PMG Bodies 

Only Central Body 

NSG Model 

None 

SRP Area 

1 (Ml 2 ) 

SRP 

Off 

Integrator Type 

RungaKutta 8(9) 

Integrator Init. 

60 (sec) 

Integrator Max. StepSize 

2700 (see) 

Integrator Accuracy 

le-13 

Report Precision 

16 significant figures 

Report CS/Cb 

Same as initial state CS 


Table 7.2: Initial Orbit Parameters (Mars Impulsive Burns) 


Initial State Parameter 

Parameter Value (unit) 

Start t, Stop Time 

01 Jan 2000 11:59:28.000 (UTCG) 

Central Body 

Mars 

Coordinate System 

Mars Mean J2000 Equator 

X 

4500.0 (km) 

Y 

0.0 (km) 

Z 

0.0 (km) 

VX 

0.0 (km) 

VY 

2.1814448386859766 (km) 

VZ 

2.1814448386859713 (km) 

Mass (No Fuel) 

850 (kg) 

Cd 

2.2 

Cr 

1.8 

Drag 

15 (m 2 ) 

Drag Model 

None 

PMG Bodies 

Only Central Body 

NSG Model 

None 

SRP Area 

1 (m 2 ) 

SRP 

Off 

Integrator Type 

RungaKutta 8(9) 

Integrator Init.. 

60 (see) 

Integrator Max. StepSize 

2700 (sec) 

Integrator Accuracy 

le-13 

Report Precision 

16 significant figures 

Report CS/Cb 

Same as initial state CS 
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Tab ie 7.3: Initial Orbit- Parameters (Moon Impulsive Bums) 


Initial State Parameter 

Parameter Value (unit) 

Start k Stop Time 

01 Jan 2000 11:59:28.000 (UTCG) 

Central Body 

Moon 

Coordinate System. 

Moon Mean J2000 Equator 

X 

2050.0 (km) 

Y 

0.0 (km) 

Z 

0.0 (1cm) 

vx 

0.0 (lan) 

VY 

1.093528701 (km) 

V'Z 

1.093528701 (km) 

Mass (No Pud) 

850 (kg) 

Cd 

2.2 

Cr 

1.8 

Drag 

15 (m 2 ) 

Drag Model 

None 

PMG Bodies 

Only Central Body 

NSG Model 

None 

SRP Area 

1 (in 2 ) 

SRP 

Off 

Integrator Type 

RungaKutta 8(9) 

Integrator Init. 

60 (sec) 

Integrator Max. StepSize 

2700 (sec) 

Integrator Accuracy 

le-13 

Report Precision 

16 significant figures 

Report CS/Cb 

Same as initial state OS 


Table 7.4: Initial Orbit Parameters (Sun Impulsive Burns) 

Initial State Parameter 

Parameter Value (emit) 

Start k Stop Time 

01 Jan 2000 11:59:28.000 (UTCG) 

Central Body 

Sun 

Coordinate System 

Sun Mean J2000 Equator 

X 

1000000.0 (km) 

Y 

0.0 (km) 

Z 

0.0 (km) 

VX 

0.0 (km) 

VY 

257.597010870 (km) 

vz 

257.597010870 (kin) 

Mass (No Fuel) 

850. (kg) 

Cd 

2.2 

Cr 

1.8 

Drag 

15 (m 2 ) 

Drag Model 

None 

PMG Bodies 

Only Central Body 

NSG Model 

None 

SRP Area 

1 (in 2 ) 

SRP 

Off 

Integrator Type 

RungaKutta 8(9) 

Integrator Init. 

60 (sec.) 

integrator Max. StepSize 

2700 (sec) 

Integrator Accuracy 

le-13 

Report Precision 

16 significant figures 

Report CS/Cb 

Same as initial state CS 
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7.2 Naming Convention 


This section describes the naming convention for Delta V scripts and output reports generated for use in GMAT’s 
Acceptance Test Plan. The naming convention consists of an ordered series of option strings, separated by 
underscores (. ). Currently, options are allowed for the following fields, and will be present in the file name in 

order: 


1. tool - The tool used to generate the trajectory. 

2. deltaVbody - The type of Delta V applied and the central body the test case uses. 

3. axes - Axes type used for burn maneuver. 

The word DeltaV precedes the tool field. The final stopping condition file format is as followed: 
D elt a' V_ too L delta Vbody.axes. report 


The tool field should always be the first option field. Future additional fields should be added to the end of the 
list of fields. Each field has a finite lust of options, as follows (future options should be added to this list): 


STK - Satellite Toolkit HPOP or Astrogator 

FF - FreeFlyer 

GMAT - General Mission Analysis Tool 

OD - Orbital Determination Toolbox 


2. deltaVbody 


IEarth - Earth centered impulsive burn 
IMars - Mars centered impulsive burn 
IMoon - Moon centered impulsive burn 
ISun - Sun centered impulsive burn 


VNB - Velocity, Velocity Normal, Velocity BiNormal axes 

Cartesian - Typical X,Y,Z,Vx,Vy,Vz axes 


7.3 Comparison Script Information 


ComparisonJDeltaV.m is the script used to perform the comparisons needed for the Delta V chapter of the 
Acceptance Test Plan. This script was designed to allow the user to select a GMAT Build and/or the exact 
solution to compare to one another. 


Refer to Appendix C for more details on this script and others used in the Acceptance Test Plan document. 


7.4 Test Case Results 


The following results are for the Delta V section. The GMAT Delta V results are being compared to the exact 
desired Delta V values. We’d like the comparison data to be as close to zero as possible. A detailed acceptance 
metric/matrix will be developed at a later date. 
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S3 


Table 7.5: Exact./GMAT lEarth DeltaV Test Case Comparison 


Test Case 

X-Pos (m) 

Y-Pos (m) 

Z-Pos (m) 

X-Vel (m/s) 

Y-Vel (m/s) 

Z-Vel (ni/s) 

Cartesian 

0 

0 

0 

0 

8.881784197e-013 

8.881784197e-013 

VNB 

0 

0 

0 

0 

2 . 6 64 5 3 5 2 o 9 e- 0 1 2 

8.881784197e-013 


Table 7.6: Exacfc/GMAT IMar.s DeltaV Test Case Comparison 


Test Case 

X-Pos (m) 

Y-Pos (m) 

Z-Pos (rn) 

X-Vel (m/s) 

Y-Vel (m/s) 

Z-Vel (rn/s) 

Cartesian 

0 

0 

0 

1. 387778781O-012 

8.8817841970-013 

0 

VNB 

0 

0 

0 

1.387778781e-012 

3.1 086244 69e-01 2 

8.8817841 97e-013 


Table 7.7: Exaet/GMAT IMoon DdtaV Tost Case Comparison 


Test Case 

X-Pos (m) 

Y-Pos (in) 

Z-Pos (m) 

X-Vcl (m/s) 

Y-Vel (m/s) Z-Vcl (m/s) 

Cartesian 

0 

0 

0 

2.775557562e-014 

0 0 

VNB 

0 

0 

0 

2.775557562e-014 

0 0 


Table 7.8: Bxact/GMAT ISun DeltaV Test Case Comparison 

Test Case” x5oT(m) Y-Pas^)" Z^Pos"(ru7 “ l-W^m/s) ^Y-Vd {xnjsT Z-Vel'(m/s) 

Cartesian 0~ 0 ~ ‘ O” 1.38777S7Sle-012 0 “ O' 

VNB 0 0 0 0 0 0 
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Chapter 8 


ControlFlow 


The Control Flow tests were designed to verify that the control flow commands (If, While, and For) function as 
expected. There are scripts that test the control flow commands by themselves, nested, and using different user 
defined parameters, such as arrays, numbers, variables, and spacecraft parameters. Each test script was designed 
to store flags that contained details of the command execution for each test case and reported to a text file. 

Due to the layout of the report, a Matlab script can easily create a pass and fail table based on the values of the 
flag variables. In the report output the main columns to pay attention to are the ranOK and ansFlag columns. 
The ansFlag variable tells us if there was an incorrect control flow execution (ansFlag— 99), correct control flow 
execution (ansFlag=l), or the control flow didn’t get executed (ansFlag=l). The ranOK variable tells us if each 
test case inside the test script ran correctly (ranOK = 1) or not (ranOK = 0). The only scripts that do not 
generate the ranOK column are the IfLoopTest##_##, IfLoopTest##, and IflfLoopTest##„## .m, because 
the ansFlag column is sufficient. 


8.1 Test Case Results 

The results in Table 8.1 display the Pass and Fail outcome of each Control Flow test script. 
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Table 8.1: Loop Test Case Results 


TestName 

Pass/Fail Failed /Total Tests 

For 

Pass 

0/15 

ForFor 

Pass 

0/9 

ForIf41-14 

Pass 

0/20 

Foiif 42-24 

Pass 

0/21 

ForIf43-34 

PaSS 

0/20 

ForIf51-15 

Pass 

0/16 

ForIf52-25 

Pass 

0/20 

ForIf53-35 

Pass 

0/20 

ForWhile42 

Pass 

0/8 

Ifll 

Pass 

0/9 

Ifl2-21 

Pass 

0/18 

If22 

Pass 

0/9 

K32-23 

Pass 

0/18 

If33 

Pass 

0/9 

If42-24 

Pass 

0/16 

I£44 

Pass 

0/8 

I:f52-25 

Pass 

0/16 

If55 

Pass 

0/8 

IfFor 

Pass 

0/9 

mf41-14 

Pass 

0/16 

If[f42-24 

Pass 

0/16 

IfIf43-34 

Pass 

0/20 

IfEfol-15 

Pass 

0/16 

I:fIfo2-25 

Pass 

0/16 

Iftf53-35 

Pass ■ 

0/16 

IfWhile 

Pass . 

0/8 

WMle42-24 

Pass 

0/16 

Wiiile43-34 

Pass 

0/16 

While52-25 

Pass 

0/16 

Wliile53-35 

Pass 

0/16 

WkileFor 

Pass 

0/9 

WmloIf41-14 

Pass 

0/16 

WliileIf42-24 

Pass 

0/16 

WMleIf43-34 

Pass 

0/16 

WMleJf51-15 

Pass 

0/17 

WMleIf52-25 

Pass 

0/16 

Whilelf53-3r> 

Pass 

0/16 

WluleTarget 

Pass 

0/1 

Wail e While42- 24 

Pass 

0/16 







Appendix A 
Acronyms 


Cb 

Cd 

COTS 

C:r 

CS 

EGM 

EOP 

FDAB 

FF 

GEO 

GMAT 

GPS 

HEO 

HP OP 

IAU 

IERS 

ISS 

JGM 

JR 

GMAT 

LEO 

LOD 

MEO 

mas 

MSISE 

NSG 

PIP 

PMG 

SRP 

STK 

TAI 

TLE 

USNO 

UT 

UTC 


- Central Body 

- Drag coefficient 

- Commercial Off Th e She lf software 

- Coefficient, j-eflexivity^A^ — " " ' ~ 

- Coordinate ^ — 

- Earth Gravity Model 

- Earth Orientation Parameters 

- Flight Dynamics Analysis Branch 

- Free Flyer 

- Geosynchronous Orbit 

- General Mission Analysis Tool 

- Global Positioning Satellite 

- Highly Elliptical Orbit 

- High Precision Orbit Propagator 

- International Astronomical Union 

- International Earth Rotation and Reference Systems Service 

- International Space Station 

- Joint Gravity Model 

- Jacchia Roberts 

- General Mission Analysis Tool 

- Low Earth Orbit 

- Length of Day 

- Medium Earth Orbit 

- milliarcseconds 

- Mass Spectrometer and Incoherent Scatter Radar Exosphere 

- Non-Spherical Gravity 

- Professional Intern Program 

- Point Mass Gravity 

- Solar Radiation Pressure 

- Satellite Tool Kit 

- International Atomic Time 

- Two Line Element 

- U.S. Naval Observatory 

- Universal time 

- Coordinated Universal 'Time 
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ACRONYMS 



Appendix B 



Initial Conditions 


The text and tables in this appendix allow someone to have an easier time duplicating the test case setups de- 
scribed in the previous chapters. The same information was presented in the previous chapters hut this appendix 
has the information arranged different for ease of duplication reasons. 


B.l Propagator Test Cases 


Table B.l: Initial Orbit Parameters (ISS) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System 

Earth Mean .J2000 Equator 

Start Time 

1 Jim 2004 12:00:00.000 (UTCG) 

Stop Time 

2 Jim 2004 12:00:00.000 (UTCG) 

X 

-4453.783586 (kin) 

Y 

-5038.203756 (km) 

Z 

-426.384456 (km) 

VX 

3.831 88S (km) 

VY 

-2.887221 (km) 

vz 

-6.018232 (km) 

Mass (No Fuel) 

1000 (kg) (some exceptions) 

Cd 

2.2 

Cr 

1.2 

Drag Area 

20 (m 2 ) 

Drag Model 

Case dependent 

NSG Model 

Case dependent 

SRP Area 

20 (to 2 ) 

SR.P 

Case dependent 

Integration Method 

Tool/Program Dependent 

Integrator StepSize control 

Fixed 

Integrator Iuit. StepSize 

5 (sec) 

Integrator Accuracy 

lc-13 

Integrator Max. StepSize 

•5 (sec) 

Report Precision 

16 significant figures 

Report StepSize 

60 (sec) 

Report GS/'Cb 

Same as initial state GS 
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Table B.2: Initial C 

Initial State Parameter 
Coordinate System 
Start Time 
Stop Time 
X 

y 

z 

vx 

VY 

vz 

Mass (No Pad) 

Cd 

Cr 

Drag Area 
Drag Model 
NSG Model 
SRP Area 
SRP 

Integration Method 
Integrator StepSize control 
Integrator Init. StepSize 
Integrator Accuracy 
Integrator Max. StepSize 
Report Precision 
Report StepSize 
Report CS/Cb 
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B.l. PROPAGATOR TEST CASES 
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Table B.3: Initial Orbit Parameters (GPS) 

Initial State Parameter Parameter Value (unit) 

Coordinate System Earth Mean J2000 Equator 

Start Time 1 Jmi 2004 12:00:00.000 (UTCG) 

Stop Time 


X 

Y 

Z 

vx 

VY 

vz 

Mass (No Fuel) 

Cd 

Cr 

Drag Area 
Drag Model 
NSG Model 
SRP Area 
3RP 

Integration Method 
Integrator StepSize control 
Integrator Init. StepSize 
Integrator Accuracy 
Integrator Max. StepSize 
Report Precision 
Report StepSize 
Report CS/Cb 


3 Jim 2004 12:00:00.000 (UTCG) 
5525.33668 (km) 
-15871.18494 (km) 
-20998.992446 (km) 
2.750341 (km) 

2.434198 (km) 

-1.068884 (km) 

1000 (kg) (some exceptions) 

2.2 

1.2 

20 (m 2 ) 

Case dependent 
Case dependent 
20 («x 2 ) 

Case dependent 
Tool/Program Dependent 
Fixed 
60 (sec) 
le-13 
60 (sec) 

16 significant figures 
120 (sec) 

Same as initial state CS 
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Table B.4: Initial Orbit Parameters (Molniya) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System 

Earth Mean J2000 Equator 

Start Time 

1 Jim 2004 12:00:00.000 (UTCG) 

Stop Time 

4 Jim 2004 12:00:00.000 (UTCG) 

X 

-1529.894287 (km) 

Y 

-2672.877357 (km) 

Z 

-6150.115340 (km) 

vx 

8.717518 (luii) 

VY 

-4.989709 (km) 

VZ 

0.0 (km) 

Mass (No Pael) 

1000 (kg) (some exceptions) 

Cd 

2.2 

Cr 

1.2 

Drag Area 

20 (m 2 ) 

Drag Model 

Case dependent 

NSC Model 

Case dependent 

SRP Area 

20 (to 2 ) 

SRP 

Case dependent 

Integration Method 

Tool/Program Dependent 

Integrator StepSize control 

Fixed 

Integrator Init. StepSize 

5 (sec) 

Integrator Accuracy 

le-13 

Integrator Max. StepSize 

5 (see) 

Report Precision 

16 significant figures 

Report StepSize 

300 (see) 

Report CS/Cb 

Same as initial state CS 
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Table E.5: Initial Orbit Parameters (GEO) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System 

Earth Mean J2000 Equator 

Start- Time 

1 Jun 2004 12:00:00.000 (UTCG) 

Stop Time 

8 Jun 2004 12:00:00.000 (UTCG) 

X 

36607.358256 (km) 

Y 

-20921.723703 (km) 

Z 

0.0 (km) 

vx 

1.525636 (km) 

VY 

2.669451 (km) 

vz 

0.0 (km) 

Mass (No Fuel) 

1000 (kg) (some exceptions) 

Cd 

Cr 

2.2 

1.2 

Drag Area 

20 hi 2 ) 

Drag Model 

Case dependent 

NSG Model 

Case dependent 

3RP Area 

20 (m 2 ) 

3RP 

Case dependent 

Integration Method 

Tool/Program Dependent 

Integrator StepSize control 

Fixed 

Integrator Init. StepSize 

60 (sec) 

Integrator Accuracy 

le-13 

Integrator Max. StepSize 

60 (see) 

Report Precision 

16 significant figures 

Report StepSize 

600 (sec) 

Report CS/Cb 

Same as initial state CS 
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Table B.6: Initial Orbit Parameters (Mars) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System 

Mars Mean J2000 Equator 

Start Time 

1 Jun 2004 12:00:00.000 (UTCG) 

Stop Time 

4 Jim 2004 12:00:00.000 (UTCG) 

X 

-2737.481646173082000 (km) 

Y 

0.0 (Ion) 

Z 

2737.481 646173082000 (km) 

VX 

-0.311321 695062643 (km) 

VY 

-3.553492313930950 (km) 

VZ 

0.311321695052650 (km) 

Mass (No Fad) 

1000 (kg) 

Cd 

2.2 

Cr 

1.2 

Drag Area 

20 (ml 2 ) 

Drag Model 

Case dependent 

NSG Model 

Case dependent 

SRP Area 

20 (m 2 ) 

SRP 

Case dependent 

Integration Method 

Tool/Program Dependent 

Integrator StepSize control 

Fixed 

Integrator Init. StepSize 

5 (sec) 

Integrator Accuracy 

le-13 

Integrator Max. StepSize 

5 (see) 

Report Precision 

16 significant figures 

Report StepSize 

300 (sec) 

Report CS/Cb 

Same as initial state CS 
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Table B.7: Initial Orbit Parameters (Mercury) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System 

Mercury Mem J2000 Equator 

Start- Time 

1 Jun 2004 12:00:00.000 (UTCG) 

Stop Time 

4 Jun 2004 12:00:00.000 (UTCG) 

X 

-2164.769322630887000 (km) 

Y 

0.0 (Ion) 

Z 

2164.769322630886100 (km) 

VX 

-0.251096955137200 (km) 

VY 

-2.866074270797602 (km) 

vz 

0.251096955137201 (km) 

Mass (No Fuel) 

1000 (kg) 

Cd 

2.2 

Cr 

1.2 

Drag Area 

20 (m?) 

Drag Model 

Case dependent 

NSG Model 

Case dependent 

SEP Area- 

20 (ni 2 ) 

SEP 

Case dependent 

Integration Method 

Tool/Program Dependent 

Integrator StepSizc control 

Fixed 

Integrator Init. StepSize 

5 (sec) 

Integrator Accuracy 

le-13 

Integrator Max. StepSize 

5 (sec) 

Report Precision 

16 significant figures 

Report StepSize 

300 (sec) 

Report CS/Cb 

Same as initial state CS 
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Table B.8: Initial Orbit Parameters (Moon) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System 

Moon Mean J2000 Equator 

Start Time 

1 .Jun 2004 12:00:00.000 (UTCG) 

Stop Time 

4 Jun 2004 12:00:00.000 (UTCG) 

X 

-I486. 7921 17191 545200 (km) 

Y 

0.0 (km) 

Z 

1486.79211 71 91543000 (km) 

VX 

-0.142927729144255 (km) 

VY 

-1.631407624437537 (km) 

vz 

0.142027729144255 (km) 

Mass (No Pad) 

1000 (kg) 

Cd 

2.2 

Cr 

1.2 

Drag Area 

20 (m?) 

Drag Model 

Case dependent 

NSG Model 

Case dependent 

SRP Area 

20 (m 2 ) 

SRP 

Case dependent 

Integration Method 

Tool /Program Dependent 

Integrator StepSiac control 

Fixed 

Integrator Init. StepSize 

5 (sec) 

Integrator Accuracy 

le-13 

Integrator Max. StepSize 

5 (see) 

Report Precision 

16 significant figures 

Report StepSize 

300 (sec) 

Report CS/Cb 

Same as initial state CS 
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Table B.9: Initial Orbit Parameters (Neptune) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System 

Neptune Mean J200Q Equator 

Start Time 

1 Jun 2004 12:00:00.000 (UTCG) 

Stop Time 

4 Jun 2004 12:00:00.000 (UTCG) 

X 

-20815.089640681723000 (km) 

Y 

0.0 (km) 

Z 

20815.089640681723000 (km) 

VX 

- 1 .426423063858300 (km) 

VY 

-16.281497481173282 (km) 

vz 

1.426423063858303 (km) 

Mass (No Riel) 

1000 (kg) 

Cd 

2.2 

Cr 

1.2 

Drag; Area 

20 (m 2 ) 

Drag Model 

Case dependent 

NSG Model 

Case dependent 

SEP Area 

20 (m 2 ) 

3RP 

Case dependent 

Integration Method 

Tool/Program Dependent 

Integrator StepSize control 

Fixed 

Integrator Init. StepSize 

' 5 (sec) 

Integrator Accuracy 

le-13 

Integrator Max. StepSize 

5 (sec) 

Report Precision 

16 significant figures 

Report StepSize 

300 (sec) 

Report CS/Gb 

Same as initial state CS 
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Table B.10: Initial Orbit Parameters (Pluto) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System 

Pluto Mean J2000 Equator 

Start Time 

1 Jun 2004 12:00:00.000 (UTCG) 

Stop Time 

4 Jun 2004 12:00:00.000 (UTCG) 

X 

-1067.516740143530600 (kra) 

Y 

0.0 (lan) 

Z 

1067.51674014352,9700 (km) 

VX 

-0.075474392S86505 (km) 

VY 

-0.861480838897026 (km) 

VZ 

0.075474392886505 (km) 

Mass (No Fuel) 

1000 (kg) 

Cd 

2.2 

Cr 

1.2 

Drag Area 

20 (m‘ 2 ) 

Drag Model 

Case dependent 

NSG Model 

Case dependent 

SRP Area 

20 (m 2 ) 

SRP 

Case dependent 

Integration Method 

Tool/ Program Dependent 

Integrator StepSize control 

Fixed 

Integrator Init. StepSize 

5 (sec) 

Integrator Accuracy 

le-13 

Integrator Max. StepSize 

5 (sec) 

Report Precision 

16 significant figures 

Report StepSize 

300 (sec) 

Report CS/Cb 

Same as initial state CS 
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Table B.ll: Initial Orbit Parameters (Saturn) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate. System 

Saturn Mean J2000 Equator 

Start Time 

1 Jun 2004 12:00:00.000 (UTCG) 

Stop Time 

4 Jun 2004 12:00:00.000 (UTCG) 

X 

-47577.347750129338000 (km) 

Y 

0.0 (km) 

Z 

47577.347750129360000 (km) 

vx 

-2.222602848522210 (km) 

VY 

-25.369834288049386 (km) 

VZ 

2.222652848522210 (km) 

Mass (No Riel) 

1000 (kg) 

Cd 

2.2 

Cr 

1.2 

Drag Area 

20 (m?) 

Drag Model 

Case dependent 

NSG Model 

Case dependent 

SEP Area 

20 (to 2 ) 

SRP 

Case dependent 

Integration Method 

Tool /Program Dependent 

Integrator StepSize control 

Fixed 

Integrator Init. StepSize 

5 (see) 

Integrator Accuracy 

le-13 

Integrator Max. StepSize 

5 (see) 

Report Precision 

16 significant figures 

Report StepSize 

300 (sec) 

Report CS/Cb 

Same as initial state CS 
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Table B.12: Initial Orbit Parameters (Uranus) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System 

Uranus Mean J2000 Equator 

Start Time 

1 Juu 2004 12:00:00.000 (UTCG) 

Stop Time 

4 Jun 2004 12:00:00.000 (UTCG) 

X 

-26762.258109447845000 (km) 

Y 

0.0 (km) 

Z 

26762.258109447823000 (km) 

VX 

-1.158158360792704 (km) 

VY 

-1 3.21 9466S691 35891 (km) 

vz 

1.158158360702704 (km) 

Mass (No Fuel) 

1000 (kg) 

Cd 

2.2 

Cr 

1.2 

Drag Area 

20 (m 2 ) 

Drag Model 

Case dependent 

NSG Model 

Case dependent 

SRP Area 

20 (to 2 ) 

SR,P 

Case dependent 

Integration Method 

Tool/Program Dependent 

Integrator StepSsze control 

Fixed 

Integrator Init. StepSize 

5 (sec) 

Integrator Accuracy 

le-13 

Integrator M^ax. StepSize 

5 (sec) 

Report Precision 

1 6 significant figures 

Report StepSize 

300 (sec) 

Report CS/Cb 

Same as initial state CS 
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Table B.13: Initial Orbit Parameters (Venus) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System 

Venus Mean J2000 Equator 

Start Time 

1 Jan 2004 12:00:00.000 (UTCG) 

Stop Time 

4 Jun 2004 12:00:00.000 (UTCG) 

X 

-4832.074380872521000 (km) 

Y 

0.0 (km) 

Z 

4832.07438087251 7400 (km) 

vx 

-0.645356787452373 (km) 

VY 

-7.366240195908405 (kin) 

VZ 

0.645356787452373 (km) 

Mass (No Fuel) 

1000 (kg) 

Cd 

2.2 

Cr 

1.2 

Drag Area 

20 (m 2 ) 

Drag Model 

Case dependent 

NSG Model 

Case dependent 

SEP Area 

20 (m 2 ) 

SEP 

Case dependent 

Integration Method 

Tool/Program Dependent 

Integrator StepSize control 

Fixed 

Integrator Init. StepSize 

5 (sec) 

Integrator Accuracy 

le-13 

Integrator Max. StepSize 

5 (sec) 

Report Precision 

16 significant figures 

Report StepSize 

300 (sec) 

Report CS/Cb 

Same as initial state CS 
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Table B.14: Initial Orbit Parameters (DeepS'pace) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System 

Sun Mean J2000 Ecliptic 

Start Time 

01 jail 2000 12:00:00.000 (UTCG) 

Stop Time 

01 Jan 2001 12:00:00.000 (UTCG) 

X 

300434.12.094803076000000 (km) 

Y 

143707423.481292670000000 (km) 

Z 

2198384.040184043300000 (km) 

vx 

-29.715920923036403 (km) 

VY 

6.056690472247896 (km) 

vz 

0.123271169290614 (km) 

Mass (No Fuel) 

1000 (kg) 

Cd 

TV 

2.2 

L.'T 

Drag Area 

X..J* 

20 (m 2 ) 

Drag Model 

Case dependent 

NSC Model 

Case dependent 

SEP Area 

20 (m 2 ) 

SEP 

Case dependent 

Integration Method 

Tool/Program Dependent 

Integrator StepSize control 

Variable 

Integrator Init. StepSize 

30000 (sec) 

Integrator Accuracy 

le-13 

Integrator Max. StepSize 

30000 (sec) 

Report Precision 

16 significant figures 

Report StepSize 

86400 (sec) 

Report CS/Cb 

Same as initial state CS 
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Table B.15: Initial Orbit Parameters (EML2) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System 

Earth Mean J2000 Equator 

Start Time 

23 Jan 2010 00:00:03.999 (UTCG) 

Stop Time 

6 Feb 2010 00:00:03.999 (UTCG) 

X 

406326.22661300009 (km) 

Y 

1 77458.38761599999 (km) 

Z 

145838.58078999998 (km) 

VX 

-0.517274673822 (km) 

VY 

0.774650366561 (Ion) 

VZ 

0.331416602654 (km) 

Mass (No Fuel) 

1000 (kg) 

Cd 

2.2 

Gr 

1.2 

Drag Area 

20 (m 2 ) 

Drag Model 

Case dependent 

NSC Model 

Case dependent 

SRP Area 

20 (m 2 ) 

SRP 

Case dependent 

Integration Method 

Tool /Program Dependent 

Integrator StepSize control 

Variable 

Integrator Init. StepSize 

1200 (sec) 

Integrator Accuracy 

le-13 

Integrator Max, StepSize 

1200 (sec) 

Report Precision 

16 significant figures 

Report StepSize 

2400 (sec) 

Report CS/Cb 

Same as initial state CS 
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Table B.16: Initial Orbit Parameters (ESL2) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System 

Earth Moan J2000 Equator 

Start Tune 

5 Feb 2006 17:05:48.772 (UTCG) 

Stop Time 

4 Aug 2006 17:05:48.772 (UTCG) 

X 

1010800.968074723 (fan) 

Y 

-910963.5377102628 (fan) 

Z 

-295145.6311353027 (fan) 

vx 

0.2642852647102676 (fan) 

VY 

0.286744175490658 (km) 

VZ 

0.07338744995264675 (Ion) 

Mass (No Fuel") 

1000 (kg) 

Cd 

2.2 

Cr 

1.2 

Drag Area 

20 (m 2 ) 

Drag Model 

Case dependent 

NSG Model 

Case dependent 

SRP Area 

20 (to 2 ) 

SRP 

Case dependent 

Integration Method 

Tool /Program Dependent 

Integrator StepSize control 

Variable 

Integrator Init. StepSize 

15000 (sec) 

Integrator Accuracy 

le-13 

Integrator Max. StepSize 

15000 (sec) 

Report Precision 

1 6 significant figures 

Report StepSize 

43200 (sec) 

Report CS/Cb 

Same as initial state CS 
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The Cb and. CS initial orbit state conditions for the ISS, GEO. Marsl, Mercuryl, Moon, Neptunel, Pluto 1, 
Saturnl, Uranusl and Venusl test cases are the same as the vales presented in Tables B.l, B.5, B.6- B.13. The 
only exception is that the Report StepSize is 600 seconds. 


Table B.l 7: Initial Orbit Parameters (Hyperbolic) 


Initial State Parameter 

Parameter Value (unit) 

Coordinate System 

Earth Mean J2000 Equator 

Start Time 

01 Jim 2004 12:00:00.000 (UTCG) 

Stop Time 

02 Jun 2004 12:00:00.000 (UTCG) 

X 

12371.701482634855 (km) 

Y 

5050.7027227610710 (km) 

Z 

5050.762722761071 (km) 

vx 

-7.9859921512608487 (km) 

VY 

2.44520073255755 (km) 

vz 

2.4452007325575495 (km) 

Mass (No Fuel) 

1000 (kg) 

Cd 

2.2 

Cr 

1.2 

Drag Area 

20 (m 2 ) 

Drag Mode] 

None 

NSG Model 

None 

SRP Area 

20 (m?) 

SRP 

None 

Integration Method 

Tool /Program Dependent 

Integrator StepSize control 

Fixed 

Integrator Init. StepSize 

5 (sec) 

Integrator Accuracy 

le-13 

Integrator Max. StepSize 

5 (sec) 

Report Precision 

16 significant figures 

Report StepSize 

600 (see) 

Report CS/Ob 

Case Dependent 
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B.3 Integrator Test Cases 


The integrator initial orbit state conditions for the ISS and GEO test cases are the same as the vales presented 
in Tables B.l and B.5. The only exception is that each test case uses one of the following Integration Methods: 
ABM, BS, PD45, PD78, RKF56, RKN68, and RKV89. 


B.4 Stopping Condition Test Cases 


Refer to Tables 5.2- 5.6 in Chapter 5 for initial orbit state values for the Stopping Condition test cases. 


B.5 Libration Point Test Cases 


Refer to Tables 6.1- 6.2 in Chapter 6 for initial orbit state values for the Stopping Condition test cases. 


B .6 DeltaV Test Cases 


Refer to Tables 7.1- 7.4 in Chapter 7 for initial orbit state values for the Performance test cases. 


B.7 Performance 


Refer to Tables ??- ?? in Chapter ?? for initial orbit state values for the Performance test cases. 



Appendix C 

Comparison Scripts Guide 


Using specific naming conventions, outlined in this acceptance test plan, and a folder architecture, highlighted 
below for the test cases, several semi-automated scripts were generated to compare all of the GMAT test case 
results with other tools. Most of these scripts have the ability to also compare results of older versions of GMAT. 


C.l Folder Architecture 

The folder architecture for the files needed for the comparison scripts is presented below: 

® GMAT_RegSetup/ 

• output / AcceptTest /CompareResults / 

- [Tooll]_[Tool2] 

• input/ AcceptTest/ 

• output /AcceptTest / [Tool] .reports / 

• output / AcceptTest /Good-reports / 

e FF 

• STK 

® GMATDocuments/ 

• AcceptTest 

C.2 Install Instructions 

1. Copy the GMAT-RegSetup and GMATDocuments folders to the same location on your hard drive. 

2. Check to make sure the folders listed in Section C.l 

3. Make sure the GMAT executable folder has a Matlab folder with the latest GMAT commands and keywords 

4. Open Matlab 

5. Set the path in Matlab to include the GMAT matlab folder 
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6 . Open GMAT and start the Matlab server 

7. Make sure Excel is closed 

8. Set current directory to the main 

GMAT uesults/ GMAT -scripts folder and then type BuildCompare_GMATteam in the command window 

or 

9 . Open one of the following files and run (F5 in Windows): 

Click ok to change the current Matlab directory if prompted 

| a) BuildRim-Script-GMAT.m 

(b) Comparison-TooIl.Tool2_PV.m 

(c) Comparison_Tooll-Tool2_CS.m 

(d) Comparison_Tooll-Tool2_Cb.m 
(c) Comparison.Tooll.Tool2-Libr.m 

(f) ComparisonJnteg.m 

(g) Comparison_DeltaV.m 
|h) Comparison-StopCond.m 

(i) TimeComparo.m 

(j) LoopTestSununaxy.m 

C.3 Warnings/Script Hints 

• The following scripts were not designed for the user to hit one button, run multiple calculations, and output 
data without user interaction. User interaction is necessary in all of these scripts. 

« The [ ] notation indicates multiple words can be used. For example, [Tool] means replace the bracketed 
expression with words such as FF,STK, and GMAT. 

• As of November 2005, the Excel data created by the comparison scripts will be saved in one file name 
[Tooll]_[Tool2]_Resiilts-[DD-MMM-YYYY].xls in their respective [Tooll]_[Tool2] folder in the CompareRe- 
sults folder. 

• Be careful adding MATLAB .m files that contain the text GMAT in their file name. These scripts use the 
text GMAT as an indicator to know if the file is a GMAT compatible script. 

• The adherence to the naming conventions are very strict in these scripts. Make sure when adding reports, 
scripts, or folders to use case sensitive filenames that agree with the naming conventions at all times. 

« Several output formats are overwritten when running each script. The Excel documents are saved with 
the current date as part of the filename. If the filename exists it will be replaced. The same overwriting 
process occurs with the Matlab .mat files. 

• In order to compare old GMAT Builds to one another, a new folder must be created with the date of 
the GMAT Build. For example, the [Month] [Day] build performed well so after running the scripts that 
generate all the comparison data for the [Month] [Day] Build, a new folder must be created. The folder can 
be named YYMMDDGMATjreports with the appropriate GMAT Build Year Month and Date replacing 
YYMMDD. Now simply copy the contents of GMAT_reports into YYMMDDGMATjreports and the data 
can be used to compare future build of GMAT to one another. 

• All the .mat files, except for the TimeComparo files, are formatted in a similar way. They contain the 
following variables mat.Toolll, mat.Tool21, mat-header, maxDiffs, and diffMat_Tooll-Tool2 or norm- 
Mat_Tooll_Tool2. 

The mat-header variable contains what parameters each column represents for the other variables The 
mat_Toolll and mat_Too!21 variables contain the report file data minus any headers. Tooll is the alpha- 
betical first tool selected and Tool2 is the other tool selected 
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• Common mistakes made with these scripts: 

- Adding report files and not following the naming convention (i.e. FFJSSJEarthJLO.O was misnamed 
as FFJSS_Earth_0_0). 

- Outputting other tools and not following the proper ordering of parameters for comparison 

- Report files were outputted in the wrong time interval increment 

- Report files were outputted without enough numerical precision. GMAT outputs data at a fixed width 
of 12 numerical characters (default). Other programs should he the same or better. 


C.4 BuildRun_Script_GMAT.m script 


Win compatible / Mac & Linux ? 


C.4.1 Purpose 


This script was designed to send multiple GMAT scripts from Matlah to GMAT to be built and ran. In early 
versions of GMAT, there was no capability to run multiple scripts, but since current versions of GMAT contain 
this capability this script is not as vital in the Acceptance test plan. 

Its secondary purpose is to record the individual time elapsed to run each test case and output the results to a 
.mat file for later use. 


C.4. 2 Inputs 

® All .m files located in the GMAT-RegSetup/input/AcceptTest/ folder that follow the GMAT comparison 
naming convention described in Sections 2.3, 3.2.1, 3.3.1, and 4.2 


C.4.3 Outputs 

• When choosing to build multiple scripts and the all command is used a mat file is created, which keeps a 
log of the time it takes to run each case. The script also displays the case name and the time it takes for 
GMAT to run the file in Matlab’s command window. 

» Any other output data is dependent on the GMAT script. 


C.4. 4 Script Algorithm 

C.5 Comparison_Tooll_Tool2_PV.m 


Win/Mac/Linux compatible 


C.5.1 Purpose 


This script is used to perform the position and velocity comparisons needed for the Propagator Section ( 2) of 
the Acceptance Test Plan. This script takes the normalized position and velocity vector differences between the 
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two selected programs. 


C.5.2 Inputs 

• Folder to search for files: [RootFolder]/[Too\]jcepoTts. 

(Refer to Section 2.3 for the naming convention of these report files.) 

« The report files must be formatted the same way. The first column is time (Mod. Julian Date), second- 
fourth columns are the position vector components ( x,y,z ), and fifth-seventh columns are the velocity 
vector components ( Vx , Vy, Vz). The data must be separated by spaces. 

• Currently the location of the report’s first row of numerical data is coded into the script. FreeFlyer is the 
fourth row, STK is the seventh row, GMAT is the first row, OD toolbox is the first row, and any other 
tools added will automatically search the first row. 


C.5.3 Outputs 

• Comparison data is displayed in MATLAB’s command window for all test cases. 

• Excel documents with comparison data and pass/fail information. 

• MATLAB .mat files with comparison data. 

® Latex documents with comparison data. 


C.5.4 Script Algorithm 

a Display welcome message 

« Display menu for Tooll options based on [Tool] .reports folder 

• Wait for user to choose tool from menu 

Implement error system for incorrect choice 

a Display menu for Tool2 options based on [Tool] .reports folder 

® Wait for user to choose tool from menu 

Implement error system for incorrect choice 

® Alphabetize Tooll and Tool2 for naming purposes and rename if necessary 

• Store the row location of the first instance of numerical data for both Tools’ report files 

• Initialize variables and folder locations 

• Display a menu of Tooll *. report files 

Generate error report if no *. report files are located in [Tooll]_report folder 

• Wait for user to choose report comparison option 

Implement error system for incorrect choice 

Open Excel Connection if compare all choice selected 

• Display filename, position difference, and velocity difference header 

» Begin Loop. Loop once for single comparison and several times for comparing all files. 

— Check the Tool2 folder for the same report 
Display error message if no report found 
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-- Continue if match found or exit loop if no match found 
— Read both output files and save the data to different matrices 
— Check to see if the row sizes are the same in both matrices 
Display error if row sizes do not match 
— Take difference of both Tools report data 
~ Normalize the results based on position and velocity 
— Determine the maximum normalized position and velocity difference 
- Store propagation duration of the test cases in a variable 
— Add acceptable differences values for Excel output 
— Save comparison data to .mat file 

* If compare all reports chosen, format data for output to Latex 

* Use BasicLatexTable script to save data to LaTex file 

* Save comparison results, acceptance errors, and duration to Excel 
— Close Excel connection if open. 

• End Loop. Allow user to rerun script 


C.6 Comparison_Tooll_Tool2_CS.m 


Win/Mac/Linux compatible 


C.6.1 Purpose 


This script is used to perform the coordinate system dependent comparisons needed for the Calculation Param- 
eters Section ( 3) of the Acceptance Test Plan. The comparison involves taking the maximum absolute value of 
the differences of the variables listed in the Inputs section of this help guide. 


C.6. 2 Inputs 

• Folder to search for files: [RootFoicier]/ [Tool]_reports. 

(Refer to Section 3.3.1 for the naming convention of these report files.) 

• The report files must be formatted the same way. Time, [X, Y, andZ] Position(km), [X, Y, andZ] Ve- 
locity (km/ sec), Mag. of Velocity(km/sec), Right Ascension of Velocity (deg), [X,Y,andZ] RxV-Specific 
Angular Momentum(fcm 2 /sec), Arg. of Perigee (deg), Declination(deg), Declination of Velocity (deg), Incli- 
nation(deg), Right Ascension (deg), Right Ascension of Ascending Node(deg) 

• Currently the location of the reports first row of numerical data is coded into the script. FreeFlyer is the 
fourth row, STK is the seventh row, GMAT is the second row, OD toolbox is the first row, and any other 
tools added will automatically search the first row. 


C.6. 3 Outputs 

• Excel documents with comparison data 

• MATLAB .mat files with comparison data 


* Latex documents with comparison data 
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C.6.4 Script Algorithm 


• Display welcome message 

• Display menu for Tooll options based on [Tool]_reports folder 

• Wait for user to choose tool from menu 

Implement error system for incorrect choice 

« Display menu for Tool2 options based on [Tool] .reports folder 

a Wait for user to choose tool from menu 

Implement error system for incorrect choice 

• Alphabetize Tooll and Tool2 for naming purposes and rename if necessary 

• Store the row location of the first instance of numerical data for both Tools’ report files 

« Initialize variables and folder locations 

« Display a menu of Tooll *. report files 

Generate error report if no *. report files are located in [Toollj_report folder 

« Wait for user to choose report comparison option 
Implement error system for incorrect choice 
Open Excel Connection if compare all choice selected 

o Display filename, position difference, and velocity difference header 

8 Begin Loop. Loop once for single comparison and several times for comparing all files. 

- Check the Tool2 folder for the same report 

Display error message if no report found 
— Continue if match found or exit loop if no match found 

- Read both output files and save the data to different matrices 
— Check to see if the row sizes are the same in both matrices 

Display error if row sizes do not match 
— Take difference of both Tools report data 

-- Determine the maximum difference for each coordinate system dependent parameter 
— Store propagation duration of the test cases in a variable 

- Save comparison data to .mat file 

* If compare all reports chosen, format data for output to Latex 

* Use BasicLatexTable script to save data to LaTex file 

* Save comparison results, acceptance errors, and duration to Excel 
— Close Excel connection if open. 

• End Loop. Allow user to rerun script 


C.T Comparison_Tooll_Tool2_Cb.m 


Win/Mac/Linux compatible 



C.7. COMPARIS ON. TO QL1.T0 0L2.CB.M 113 

C.7.1 Purpose 

This script is used to perforin the central body dependent comparisons needed for the Calculation Parameters 
Section ( 3) of the Acceptance Test Plan. The comparison involves talcing the maximum absolute value of the 
differences of the variables listed in the Inputs section of this help guide. 


C.7.2 Inputs 

• Folder to search for files: [RootF older]/ [Tool] jreports. 

(Refer to Section 3.2.1 for the naming convention of these report files.) 

• The report files must be formatted the same way. Time, Altitude (km), Beta Angle (deg), C3JBnergy 
(km 2 / sec 2 ), Eccentricity, Latitude (deg), Longitude (deg), (RxV)-Mag (km 2 /sec), Mean Anomaly (cleg), 
Mean Motion ( rad/sec ), Period (sec), Apoapsis Radius (km), Perigee Radius (km), R_Mag (km), Semi- 
major Axis (km), True Anomaly (deg), Semilatus Rectum(km), Apoapsis Velocity (km/ sec), Periapsis 
Velocity (km/ sec), Greenwich Hour Angle(deg), Local Sidereal Time 

• Due to the inability of FF and STK tools to output all the parameters listed in the previous bullet, excep- 
tions for FF and STK are built into the code. 

STK: Semilatus Rectum, Apoapsis Velocity, Perigee Velocity, Greenwich Hour Angle,- and Local Sidereal 
Time can not be outputted in the same report file. Out of the aforementioned parameters Greenwich 
Hour Angle is the only parameter that can be outputted in a separate file. All the other parameters are 
calculated in MATLAB based on the results STK could generate. 

FF: (RxV) JVIag and R_Mag could not be outputted easily. Instead the [XYZ] components were outputted 
and the script computes the magnitude of the vectors. Apoapsis Velocity, Periapsis Velocity, and Local 
Sidereal Time are all created in this script based on available parameters. 

Modifications to the code will need to be made to add a new tool that cannot output all of the central 
body parameters that GMAT does. 


• Currently the location of the reports first row of numerical data is coded into the script. FreeFlyer is the 
fourth row, STK is the seventh row, GMAT is the second row, OD toolbox is the first row, and any other 
tools added will automatically search the first row. 


C.7.3 Outputs 

® Excel documents with comparison data 

• MATLAB .mat files with comparison data 

• Latex documents with comparison data 


C.7.4 Script Algorithm 

® Display welcome message 

• Display menu for Tooll options based on [Tool] .reports folder 

• Wait for user to choose tool from menu 

Implement error system for incorrect choice 
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« Display menu for Tool2 options based on [Tool] .reports folder 

« Wait for user to choose tool from menu 

Implement error system for incorrect choice 

• Alphabetize Tooll and Tool2 for naming purposes and rename if necessary 
a Store the row location of the first instance of numerical data for both Tools’ report files 
« Initialize variables and folder locations 

a Display a menu of Tooll *. report files 

Generate error report if no *. report files are located in [Tooll]_report folder 

a Wait for user to choose report comparison option 
Implement error system for incorrect choice 
Open Excel Connection if compare all choice selected 

a Display filename, position difference, and velocity difference header 
a Begin Loop. Loop once for single comparison and several times for comparing all files. 

- Check the Tool2 folder for the same report 

Display error message if no report found 
— Continue if match found or exit loop if no match found 

- Read both output files and save the data to different matrices 
-- Check to see if the row sizes are the same in both matrices 

Display error if row sizes do not match 
— Take difference of both Tools report data 

— ■ Code in exceptions for STK, FF, and any other tools that don’t output all the desired GMAT central 
body dependent parameters 

— Determine the maximum difference for each central body dependent parameter 
— Store propagation duration of the test cases in a variable 
— Save comparison data to .mat file 

* If compare all reports chosen, format data for output to Latex 

* Use BasicLatexTable script to save data to LaTex file 

* Save comparison results, acceptance errors, and duration to Excel 
-- Close Excel connection if open. 

a End Loop. Allow user to rerun script 


C.8 Comparison_Tooll_Tool2_Libr.m 


Win/Mac/Linux compatible 


C.8.1 Purpose 


This script is used to perform the position and velocity comparisons needed for the Libration Points Section ( 6) 
of the Acceptance Test Plan. This script takes the normalized position and velocity vector differences between 
the two selected programs. 
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C.8.2 Inputs 

• Folder to search for files: [Root Folder]/[Tool]_reports. 

(Refer to Section 6.2 for the naming convention of these report files.) 

• The report files must be formatted the same way. The first column is time (Mod. Julian Date), second- 
fourth columns are the position vector components ( x,y,z ), and fifth-seventh columns are the velocity 
vector components (Vx,Vy,Vz). The data must be separated by spaces. 

• Currently the location of the report’s first row of numerical data is coded into the script. FreeFlyer is the 
fourth row, STK is the seventh row, GMAT is the first row, OD toolbox is the first row, and any other 
tools added will automatically search the first row. 

C.8.3 Outputs 

• Comparison data is displayed in MATLAB’s command window for all test cases. 

• Excel documents with comparison data and pass/fail information. 

• MATLAB .mat files with comparison data. 

• Latex documents with comparison data. 


C.8.4 Script Algorithm 

[INSERT script Algorithm] 


C.9 Comparison_Integ.m 


Win/Mac/Linux compatible 


C.9.1 Purpose 


This script is used to perform the integrator comparisons needed for the Integrator Section ( 4) of the Accep- 
tance Test Plan. The comparison involves taking the difference of the position and velocity vector and then 
normalizing these two vectors to get the position and velocity difference. This script behaves similar to the 
Comparison_Tooll_Tool2JPV.m script but the components being varied are the integrators for two body test 
cases. 


C.9. 2 Inputs 

* Folder to search for files: [RootF older] /[Tool]- reports. 

® Naming convention: Integrator-[Tool]-[Trajectory]_[/nte(;raforTy_pe]_2Body.report 

• The report files must be formatted the same way. The first column is time, second-fourth columns are 
the position vector (x,y,z), and fifth-seventh columns are the velocity vector (x,y,z). The data must be 
separated by spaces. 

® Currently the location of the reports’ first row of numerical data is coded into the script. FreeFlyer is the 
fourth row, STK is the seventh row, GMAT is the second row, OD toolbox is the first row, and any other 
tools added will automatically search the first row. 
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C.9.3 Outputs 

« Excel documents with comparison data into the following folder: 

[RootF older] /Com.pareR,esvHts/[Tooll\-[Tool2] 

• MATLAB .mat files with comparison data into the following folder: 

[RootFol der] / CompareResults/ [T ooll ]. [Tool2] 

« Latex documents with comparison data into the following folder: 

[RootF ol der]/LatexJDocs 


C.9.4 Script Algorithm 

« Display welcome message 

« Display menu for Tooll options based on [Tool] .reports folder (Can only be Exact or GMAT folders) 

a Wait for user to choose tool from menu 

Implement error system for incorrect choice 

® Display menu for Tool2 options based on [Tool] .reports folder (Can only be Exact or GMAT folders) 

« Wait for user to choose tool from menu 

Implement error system for incorrect choice 

« Alphabetize Tooll and Tool2 for naming purposes and rename if necessary 
s Store the row location of the first instance of numerical data for both Tools’ report files 
• Initialize variables and folder locations 

« Display a menu of Tooll *. report files 

Generate error report if no *. report files are located in [Tooll]_report folder 

u Wait for user to choose report comparison option 
Implement error system for incorrect choice 
Open Excel Connection if compare all choice selected 

® Display filename, position difference, and velocity difference header 
® Begin Loop. Loop once for single comparison and several times for comparing all files. 

— Check the Tool2 folder for the same report 
Display error message if no report found 
— Continue if match found or exit loop if no match found 
— Read both output files and save the data to different matrices 
— Check to see if the row sizes are the same in both matrices 
Display error if row sizes do not match 

- Take difference of both Tools report data 

- Normalize the results based on position and velocity 

— Determine the maximum normalized position and velocity difference 
— Store propagation duration of the test cases in a variable 
— Save comparison data to .mat file 

* If compare all reports chosen, format data for output to Latex 

* Use BasicLatexTable script to save data to LaTex file 

* Save comparison results, acceptance errors, and duration to Excel 
— Close Excel connection if open. 

« End Loop. Allow user to rerun script 
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C.10 

Comparison JDeltaV.m 

Win/Mac 

:/Linux compatible 

C.10.1 

Purpose 

C.I0.2 

Inputs 

C.10. 3 

Outputs 

C.10.4 

Script Algorithm 

C.ll 

Comparison_StopCond.m 

Win/Mac/Linux compatible 

C.ll.l 

Purpose 

C.ll. 2 

Inputs 

C.ll. 3 

Outputs 

C.11.4 

Script Algorithm 

C.12 

STKJRepropagate.m 

Win compatible 


C.12.1 Purpose 


The STK_Repropagate script was designed to reduce the time it took to generate STK report files, after mod- 
ifications to the STK scenario were made, and obtain more accurate STK run times. Through STK’s connect 
module Matlab connect with STK and propagates satellites, generates reports, and outputs run times. 


C.12.2 Inputs 

b STK scenario folders that follow the GMAT Acceptance Test Plan naming convention, in the following 
folder: [KtfofFolderJ/TVuthFiles/STK 
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C.12.3 Outputs 

• STK report file saved into [RootF older]/ STK-reports. 

« Matlab .mat file with the time taken to propagate each satellite. 


C.12.4 Script Algorithm 

C.13 TimeComparo.m 


Win/Mac/Linux compatible 


C.13.1 Purpose 


When running the BuildRunJ3cript_GMAT.nl script there is an all option after selecting the build & rim multiple 
cases choice. By using this all option, the GMAT performance times for all the test cases are saved to a .mat 
file. This script uses those saved performance times and, based on a pre-selected amount of test cases, creates a 
new excel file that contains GMAT, FF," and STK performance times for those pre-selected cases. 


C.13. 2 Inputs 

* Template file containing pre-selected cases: [RootF older] /NonGMATrunTimes.xls 

a Folder to search for files: [RoofFolderJ/CompareResults 

a Naming convention: 

[Date]/ _Time2RunAll . mat 


C.13.3 Outputs 

® Excel document with comparison data 


C.13.4 Script Algorithm 



Appendix D 


STK Setup 


The STK setups were very crucial in determining a preliminary standard to compare GMAT to. In the initial 
stages of the Acceptance Test Plan STK scenarios were obtained from Emergent Space, as mentioned in Section 
2.1. These Earth based cases were created in STK-HPOP and modified to provide a setup that was as equivalent 
to GMAT as was possible. Non-Earth test cases could not he created with STK-HPOP, so STK-Astrogator was 
used. In order to use STK-Astrogator several Astrogator elements needed to be created. Section D.4 details all 
the elements needed for the non-Earth STK test cases. 


D.l Support Files Needed 


All alterations to STK scenarios were performed with STK 6.1 on the desktop machine described in the Perfor- 
mance Section (Section ??). 

In order to duplicate data generated by STK for this Acceptance Test Plan, the same files presented in Table 
D.l are needed. 


Table D.l: STK support files 


Filenainc(s) 

Filc( ! s) and/or Foldcr(’s) Location 

Information in Section 

EOP.dat 

[STK Root Directory] \DynaniieEarthData 

2.2.1 

EOP.dat.all 

[STK Root D irectory] \ D ynamicEart hD at a 

2.2.1 

All planetary cb files 

[STK Root Directory]\STKData\Centra!Bodies 

2.2.2 

Non-sphcrical grav (nsgj files 

[STK Root Dircctoryj\STKData\GcntralBodics 

2.3 

Report Styles 

[STK individual user folder] \Config\Stylcs 

D.5 

Non-Earth Astrogator Elements 

[STK Root Directory] \8TKData\ Astrogator 

D.4 

Leap Seconds File 

[STK Root D ireeto ry]\STKD at a\ Astro 

2.2.2 


D.2 STK modules used 


Refer to Appendix B for a easy initial orbit state format to use as input into STK. Table D.2 displays the STK 
module to use in order to duplicate the results seen in theis Acceptance Test Plan document. 
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Table D.2: STK modules used 


Test Cases 

Test Group (s) 

STK-type 

ISS. GEO, SunSync, Molniya, GPS, Hyperbolic 
DeepSpaee, EML2, ESL2. All non-Earth cases 
ISS, GEO 

Propagator, Cb, CS 
Propagator 
Integrator 

STK- HP OP 
STK-Astrogator 
STK-ZBody 


D.3 Scenario Setup 


The following are guidelines we followed when creating STK Scenarios, excluding HPOP or Astrogator specific 

guidelines: 

ci Initial Scenario Epoch is the same as the satellite epoch 

« Make sure all support files (Table D.l) are present and match the information in this document. 

® *.cb planetary files need to contain the values in Tables 2.5, 2.7, and 2.9 

« Remember the Modified Julian Date used in GMAT and STK are different. 

(i.e. UTC-GMATJVfodJulian = UTC-STK-ModJulian + 29999.5) 

D.3.1 HPOP 

a Planetary information should be taken from the JPL DE file 
® Satellite properties must be consistent with table the values in Table 2.5 


[INSERT Screen captures of HPOP] 


D.3. 2 STK 

* All Astrogator support files (Table D.l) are essential in reproducing the values in this document. 

« Make sure the maximum propagation is greater than the test case propagation duration or turn off the 
feature in STK. 


[INSERT Screen captures of Astrogator] 
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D.4 Astrogator Elements 


For the non-Earth cases a large number of Astrogator elements were created in order to compare STK to GMAT’ 
results. 


D.4.1 Calculation Objects: Cartesian Elements 


Vx_EaxthMJ2000Ec 

VxJEaxthMODEe 

Vx_EarthMODEq 

Vx_E ar thM 0 EEc 

Vx_EarthMOEEq 

VxJEarthTODEc 

Vx_E arthT 0 D Eq 

Vx_E arthT 0 EEc 

VxJEarthTOEEq 

VyJEarthMJ2000Ec 

Vy JE arthM 0 DEu 

VyJEarthMODEq 

Vy JE arthM 0 EEc 

Vy_E arthM 0 EE q 

Vy_EarthTODEc 

V y JE arthTO D Eq 

VyJEarthTOEEc 

VyJEartbTOEEq 

Vz_EarthMJ2000Ec 

VzJEarthMODEc 

VyJEarthMODEq 

VzJEarthMOEEc 

VzJEarthMOEEq 

Vz_EarthTODEc 

Vz-EarthTODEq 

Vz.EarthTOEEc 

VzJEarthTOEEq 

Vx_Mars_Centered_Mean_J2000_Earth_Ec 

Vx_Mars-Centered_Mean_J2000_Earth_Eq 

Vy_Mars_GenteredJvlean_J2000-Earth_Ee 

Vy_Mars_Centered_Mean_J20002EaxthlEq 

Vz_Mars-CenteredJMean_J2000_EfirthJEc 

Vz_Mars_CenteredJ\'lean_J2000-E;irthJEq 

Vx_MarsFixed 

Vy_MarsFixed 

Vz-MarsFixed 

VxJderairy_CenteredJVIean-J2000_Earth_Ee 

VxAlerciny .Centered Jvlean_J2flOOJEarthJEq 

Vy_Mercnry_Centered_Mean_J200QJEarth_Ec 

Vy-MercHry-Centered_Mean.-J20002BartbJBq 

Vz_Mercury_Centered_Mean_J2000_EartliJEc 

VzJvIercury_OenteredJVIean-J20002E;irthJEq 

Vx_Mcr curyF ixcd 

Vy_Mcr curyF ixcd 

VzJvlercuryFixed 


VxJNcptunc-CcntcrcdJV[can_J200Q-Earth_Ec 

Vx JN T eptune -Centered JVIean_J 2000 _E arth JEq 

Vy-Neptune-CenteredJMean_J2000_EarthJEc 

Vy .Nep tune.Cen tered_Mean_J2000 .Earth JEq 

VzJs T eptune_Centered_Mean_J2000JEarth_Ec 

Vz-Neptune_Centered_Mean_J2000_Eaxth_Eq 

Vx_N eptuneF ixed 

VyJSTeptuneFixed 

Vz JN 7 cptuncFixcd 

VxJ?hito-CenteredA'Iean_J2000JEarthJEc 

Vx_Phito_Centered-Mean.J2000-Earth-Eq 

Vy _Pluto-Centered_Mean_J 2000 -Ear th_Ec 

Vy_Pluto_Centered_Mean_J2000-Earth_Eq 

V2_Pluto_CenteredJVlean.J2000-Eart;h-Ec 

Vz-Pluto_Centered_Mea:)-J200(}_Eari;h-Eq 

VxJPiutoFixed 

VyJPlutoFixed 

Vz_PlutoFixcd 

Vx_Satum_Centered_Meaii_vT2000_EarthJEc 

Vx_Satum_CenteredJVlean_J2000-Eaxth_Eq 

Vy_Saturn_Centered_Mean_J20()0_Earth_Ee 

Vy-Saturn_Centered_Mean_,T20C)0-EarthJEq 

Vz_Satuni_CenteredJMean_J2000-Earth_Ec 

Vz-Satum_CenteredJMean-J2000Jiirrth-Eq 

Vx-SaturnFixed 

Vy.S atuniFixe d 

VzJS atxirnF ixcd 

Vx_Sun_Centered_Mean.J20Q0-Earth_Ec 

Vx-Sun..CenteredJvlean.-J2000JBartl:-Eq 

Vy-Sun_Gentered_Meaoti-J2000_Eartk_Ee 

Vy_Sun_Centered_Mean_J2000_Eartk_Eq 

Vz-Srm_CenteredJMean_J2000-Eartb_Ec 

Vz-Sim_GenteredJMe«'m_J2000JEarth-Eq 

Vx_SunFixed 

VyJSunFixcd 

Vz-SunFixed 

Vx_UranuH_Centered_Meaix_J2000_Earfch_Ec 

Vx_Uranus_Centered_Mean_J2000_Earth_Eq 

Vy_Uranus_Centered_Mean_J2000_Earth_Ek: 

Vy_Urfinus_Centered_Meim_J2000-Earth_Eq 

Vz_Urauus-CenteredJVTemi_J2000_Earth_Ec 

Vz_Uniims_CenteredJVIeaii-J2000-EarthJEq 

Vx.UranusFixed 

V y _Ur auusFixcd 

Vz-Ur anusFixed 


122 


APPENDIX D. STK SETUP 



VxJi'loou_CciitcrcclJVloaii-J2000_EarthJEc 

VxJ'.loou_Gcutcrcd_Moaii.J2000_EaxlliJEq 

Vy_Aloon_CeiLtered_M[eaxi-J2000-EarthJEc 

Vy_Moon-CenteredJVlean-J2000_Earth_Eq 

Vz-Moun_Centered_Mean_J200GJEarth_Ee 

Vz_Moon_Centered_Mean.J2000_Earth_Eq 

Vx-MooD.Fbr.ed 

VyJMoouFixed 

VzJVIoonFixed 

XJEarthGSE 

Y-EaiilxGSE 

ZJEarthGSE 

VxJSarthGSE 

VyJEarthGSE 

Vz-EarthGSE 



Vx_Vcmis-Gciitcrcd_Mcau_J2000-EartliJEc 
Vx_Vcitus_C entered JVIcan_ J 2000 -Earth. _Eq 
Vy_Veiius-Centered_Meaii_J2000_Eartli_Ec 

V y, Venus .0 entered_Mean_J2000 _E ar th_Eq 
Vz_Venus_CenfceredJMean_.T2000_Eartli_Et: 
Vz_Vemis_CenteredJ'.leaii_,] 2000 -Earth-Eci 
Vx..VenusFixed 

V y.Vem isF ixed 
Vz_VemisFixed 
X-EaxthGSM 
Y_EarthGSM 
Z_EarthGSM 
VxJSarthGSM 
VyJSarthGSM 
Vz-EarthGSM 


D.4.2 Calculation Objects: Geodetic Elements 


Altitude_Mars 

LatitudeJVlars 

Longit net e JVtars 

Altitude-Mercury 

Latitude-Mercury 

Longitude-Mercury 

AltitudcJVIoon 

Latiftick-Moon 

Longitude_Moon 

Altifcude-Neptune 

Latitude_Neptune 

Longitude-Neptune 


Altitude_Pluto 

Latitude_Pinto 

Lon gitud e. JPluto 

Altitude-Satum 

Latitude_Saturn 

Longitude-Saturn 

Altiiudc.U rams 

LatitudeJJrauus 

Longitude-LIranus 

Altitude_Vemis 

Latitude.. Venus 

Longitude-Venus 



DA. ASTROGATOR ELEMENTS 

D.4.3 Calculation Objects: Keplerian Elements 


Argiiment.ofJPeriapsis_(fixed) 

Argunient_of_Periapsis.(MJ2000Ec) 

Argument.of2Periapsis_(MODEc) 

Argumont.of_Periapsis_(MOEEc) 

Argument_of_Periapsis_(MOEEq) 

Argument_ofJPeriapsis.(TODEc) 

Argiiment_ofJPeria,pyis_(TOEEc) 

Argument_of.Periapsis_(TOEEq) 

Argument_of2Periapsis_(EarthGSE) 

Argument_of_Periapsis_(EarthGSM) 

ArgTunent_ofJPeriapsis.(M<xrsFixed) 

Argumcnt_of_Pcria.psis_(MarsMJ2000Ec) 

Argxme-ut_ofJPcriapsisJVlars_(MJ2000Eq) 

Argument_of_Periapsis.(MerciiryFixed) 

Arguns ent_ofJPeriapsis_(MercuryMJ2QOOEo) 

Argiiment_of_PeriapKisJ\lercury_(MJ2000Eq) 

Argument_of_Periapsis_(MoonFixed) 

Argiiment.of.Periapsis_(MoonMJ2000Ee) 

Argument_of-Peri apsis JVlocm_(MJ2000Eq) 

Argument_qf _Periapsis_(NeptuueFixed) 

Argiuncnt_ofJPcriapsis.(NcptuncMJ2000Ec) 

Argument_of.PeriapsisJ^eptxtne_(MJ2000Eq) 

Argument, _ofJPeriapsis-(PlutoFixed) 

Argument_ofJPeriapsis.(PhifcoMJ2000Ee) 

Argiiment_of2Periapsis_Phito_(MJ2000Eq) 

Argument.oLPeri apsis. ( SatumFixed) 

Argument_ofJ , eriapsis.(SatumMJ2000Ec) 

Argument_ofJPeriapsis.Saturii_(MJ2000Eq) 

Argument.of JPeriapsis..(UramisFixed) 

Argumcnt_of2Pcriapsis_(UranxisMJ2000Ee) 

Argument_of_Periapsis_Urarms-(MJ2000Eq) 

Argument.of J 3 eriapsis_(VenusFixed) 

Argument.oLPeri apsis. (VenusMJ2000Eo) 

Argxunent_of_Periapsis_Vemis_(MJ2000Eq) 

Eccentricity .Mars 

Eccentricity .Mercury 

Eccentricity Moon 

Eccentricity .Neptune 

Eccentricity _P luto 

Eccentricity 23 atum 

Eccentricity .Uranus 

Eccentricity .Venus 


Indi.nation_(fixed) 
Inclination_(MJ2000Ec) 
Inclination_(M ODEc) 
Inclination.(MOEEc) 

Inclinati o n_(M OEEq) 

Inclination.(TODEe) 

Inciinat.ion-(TOEEc) 

Inelination.fTOEEq) 

Inclinatiori_(EaxthGSE) 

Inclination.(EarthGSM) 

Inclination_(MarsFixed) 

.Tnclinatiou_(MarsMJ2000Ec) 

Inclination2Vlars_(MJ2000Eq) 

Inciination_(MercuryFixed) 

In.dutation_(MercuryMJ2000Ec) 

Inclination Jtfemiry.(MJ2000Eq) 

Inclination.(MoonFixed) 

Inclinati on_(MoonM J 2000Ec) 

IncLination_Moou_(MJ2000Eq) 

Inclinatioii.(NeptuneFixed) 

Inciination_(NoptuncMJ2000Ec)' 

Inclination_IS T eptime_(MJ2000Eq) 

Inciination_(PlutoFixed) 

Inclination_(PlutoMJ2000Ec) 

Inclination Jduto.(MJ20(K)Eq} 

Inclination_(Sa,turnFixed) 

Inclination_(S atuniM J2000Ec) 

Inclination_Satum_(MJ200QEq) 

Inclination.. (UranusFixed) 

Inclination.(UrarmsMJ2000Ec) 

Inclination 

Inclination. (VenusFixed) 

Inclinati c >n.( Vein xsM J 2f )00Ec) 

Inclination_Venus_(MJ2000Eq) 

Mean_Motion_Mars 

Mean_Motiou..Mercury 

Meim.Motion.Moon 

Mean JMotioiiJNcptuuc 

Mcan_Motion_Plxito 

MeanJVIotion-Satin-n 

MeanJVIotion.Uranus 

Meem_Motion_Venus 
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Orbit J?eriod_Mars 
Orbit_Pcriod_Mcrcury 
Orbit_Period Jdoou 
0 rbit -Per i o d_Nep t line 
OrbitJPeriod-Pluto 
Orbit JPeriodJSatum 
Orbit-Period-Uramis 
Orbit-Period_Venas 
RAAN.(fixed) 

RAAN_(M J 2000Ec) 

RAAN.(MODEc) 

RAAN.fMOEEc) 

RAANAMOEEq) 

RAAN.(TODEc) 

RAAN.(TOEEc) 

RAAN-(TOEEq) 

RAAN-(MarsFixed) 

RAAN_(MarsMJ2000Ec) 

RAANJMars_(MJ2000Eq) 

RAAN_(McrcuryFixcd) 

RAAN_(MercuryMJ200QEc) 

RAAN_Mercury_(MJ2000Eq) 

RAAN-(MoonFixed) 

RAAN-(MoonMJ2000Ec) 

RAANJVIoon.(MJ2000Eq) 

RAAN-(NeptuneFixed) 

RAAN_(NeptuneMJ2000Ec) 

RAAN_Ncptune_(M J2000Eq) 

RAAN-(PlutoFixed) 

RAAN_(PlutoMJ2000Ec) 

R A AN _Plu t o _ (M J2000Eq) 

RAAN.(SatiimFixed) 

RAAN.(SatumMJ2000Ec) 

RAAN-Saturn-(MJ2000Eq) 

RAAN-(UranusFixed) 

RAAN- (UrauusM J 2000Ec) 

RAAN-Uranus_(MJ2000Eq) 

RAAN _ (VenusF. bred) 

RAAN-(VenusM J2000Eei 

RAAN_Vemis_(M.T2000Eq) 


Radius_OLApoapsis_M ars 

Radius_OLApoapsisJVIercury 

Radius-OLApoapsisJYloon 

Radius-OLApoapsis-Nepfcune 

Radius-Of-Apoapsis-Pluto 

Radius_Of_Apuapsis_Satiim 

Ra< Uus_OLApo apsis_U r anus 

R.aii\is.Of-Apoapsis-Venus 

Radius_Of_Periapsis-Mars 

Radius-OfJPeriapsis-Mercury 

Radius-OLPcriapsis-Moon 

Radius-OLPeriapsis-Neptune 

Radius-Of-Periapsis -Pluto 

R.adiua_0£_Periapsis_Saturn 

Radius-OfJPeriapsis.Uranus 

Radius-OtPeriapsis-Venuo 

Semimajor_Axis_Mars 

Semimajor -Axis-Mercury 

Scmbnajor-Axis_Moon 

Sciniinajor-Axis -Neptune 

Seinimajor_Axis_Pluto 

•Sermmajor_Axis -Saturn 

Semimajor _Axxs_Uranus 

Semimajor -Axis_Venus 

True. Anomaly _M ars 

True_Anomaly_Mercurv 

True_Auomaly_Moou 

True_Anomaly_N eptuue 

True_Anomaly -Pluto 

True-Anomaly .Saturn 

True-Anomaly-Uranus 

Trae_Anomaly_ Venus 

RAAN.(EarthGSE) 

RAAN . ( E axtb GSM) 



DA. ASTROGATOR ELEMENTS 

D.4.4 Calculation Objects: Other Orbit Elements 


Beta_Angle_Mars 

Beta^AngleJVIercury 

Beta_Angle_Moon 

Bet a_Anglc_N optima 

Beta_AngleJPluto 

Beta_Angle_Satum 

Beta.Angle.Uranus 

Beta_Angle_Venus 


D.4.5 Calculation Objects: Target Vector Elements 

C3JEnergy_Mars 
G 3-Energy .Mercury 
C3_Energy_Moon 
CSJEnergyJSTeptune 
G 3 -Energy JP luto 
03_Eiiergy .Saturn 
03-Eiicrgy_Urainis 
G 3 .Energy. Venus 


D.4.6 Calculation Objects: Spherical Elements 


Declination. (EartliGSE) 

Declination_(Earth.GSM) 

Declination_(MJ2000Ec) 

Declination.(MODEe.) 

Declination.(MOEEc) 

Declination. (MOEEq) 

D eclination. (M arsF ixed) 

D eclination. (M arsM J2000Ec) 

D ccli nation. (M er cur y Fixed) 

D eclination. (McrcuryM J 2000Ec) 
Declination_(MoonFixed) 

D eclination. (MoonMJ2000Ec) 

D eclination. (N ep tnneF ixed) 

D eclination_(N eptuneM J2000Ec) 
D eclination. (P lutoFixed ) 
Declination_(PhitoMJ2000Ec) 

D eclination. (S aturnM J2000Ec) 

Dcclination_(SaturusFixcd) 

Dedination.(TODEc) 

D eclination. (T OEE-c) 


Rigid -Asc-(MercuiyMJ2000Ec) 

Right-Asc.(MoonFixed) 

Right_Asc_(MoonMJ20G0Ec) 

Right_Asc_(NeptuneFixed) 

Right„Asc_ (NeptiineM J2000Ec) 

Right-Ase_(PIutoFixed) 

Right_Asc_(PiutoM J 2000Ec) 

Right. Asc.(SaturnFixed) 

Right_Asc_(S aturnM J2000Ec) 

Right_Asc_(TODEc) 

Right_Asc_(TOEEc) 

Right_A.sc. (TOEEq) 

Righ t_Asc_(Ur ani isF ixed) 
Right_Asc_(UranusMJ2000Ec) 
Right _Asc_ (VenusFixed ) 
Right_Asc.(VennsMJ2000Ec) 
Right_Asc_Mars.(MJ2000Eq) 
Right -Asc.Mercury_(MJ2000Eqj 
Right.Asc.Moou.(MJ2000Eq) 
Riglit_AscJS T eptune_(MJ2000Eq) 
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Declination_(TOEEq) 
Declination_(UraamsFixed) 
Declination.(UrajmsMJ2000Ec) 
Declination. (Ve liusF ixed ) 
Declinatiou_(VenusMJ2000Ec) 

D cclinationJVI ars . (M J 2000Eq) 

Declination2Mercury_(MJ2000Eq) 

Declination3I o on_(M J2000Eq) 

DeclinationJNfeptnne_(MJ2000Eq) 

Declmation_PLuto_(MJ2000Eq) 

Declination_Saturn_(MJ2000Eq) 

Declination_Uramis-(MJ2000Eq) 

Declination_Venns_(MJ2000Eq) 

R31ag_(MarsFixed) - - - 

RJMag_(MercuryFfaccd) 

R31ag_(MoonFixed) 

R _Mag_(N eptuneFixed) 

R_Mag_(PlutoFixed) 

R_Mag_{SaturnFixed) 

R31ag.(UranusFixed) 

R-Mag-(VenusFixed) 

R31ag_Mars 

R.JN-lagJvlcrcury 

R31ag_Moon 

R.Mag .Neptune 

RJ\'lagJPluto 

R31ag_Saturn 

R31ag_Uranus 

RJvlag-Venus 

R_Mag_(EarthGSE) 

R31ag.{EarthGSM) 

Ri glit . Asc _ (M J2000Ec) 

Rightist, (MODEc) 
Right_Asc_(MOEEc) 

Right jLHc.(MOEEq) 

Right _Asc_(MarsFixed) 

Right 3sc.(MarsMJ200()Ec) 

Right _Asc_ (MercuryFixed) 

Rj ght -Asc: . (EarthG SM) 


Right_AscJPluto_(M.T2000Eq) 
Right Jksc_Saturn_(MJ2000Eq) 
Right_Aac_Uratms_(MJ200GEq) 
Right -Asc_Vemis.(MJ2000Eq) ' 
vjvlag_(MJ2000Ec) 
VJVlag.(MODEc) 

V31 ag_ ( M OEEc) 
VJViag.(MOEEq) 
V31ag.(MarsFixed) 

V31 ag.(M arsM J 2000Ec) 

V31 ag.(MercuryFixed) 

V 31 ag. (Mereur yM J2000Ec J 

V31ag_(MoonFixed) 

V31ag_(MoonMJ2000Ec) 

V Jdag-fNeptmieF ixed) 
V31ag.(NeptimeMJ2000Ec) 

V 31 ag. ( PlutoFixed) 
V31ag_(PIutoMJ2000Ec:) 

V 31 ag . ( S at nmM J2 OOOEc ) 
V31ag.(Sat\irnsFixed) 
V31ag.(TODEc) 
V31ag.(TOEEe) 
V31ag.(TOEEq) 
V3iag_(UranusFixed) 

V.M ag. (UranusM J 2000Ee) 
V31ag.(Venua'Fixed) 
V31ag_(VenusMJ2000Ec) 
V3iag_Mars_(MJ2000Eq) 

V 31ag31ercury_(M J 2000Eq) 
V31ag.(EarthGSE) 

V 31 ag.(EarthGSM) 
V31ag3foon.(MJ2000Eq) 
V3fagJN T eptime_(MJ2000Eq) 
V31agJ?luto.(MJ2000Eq) 
V34ag^atum_(M.T2000Eq) 
V31ag_Uranus.(MJ2000Eq) 
V31ag_Verais_(MJ2000Eq) 
Right_Asc_(EarthGSE) 
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D.4.7 Calculation Objects: Vector Elements 


(RxV)_Mag 

(RxV)JMag-(MarsFixed) 

(RxV) JVlag.(MarsM J2000) 
(RxV) _Mag_(MarsMJ2000Ec) 
(RxV) JVlag-(MercuryFixed) 
(Rx\0 -Mag.(MercuryM J2000) 
(RxV) JVlag_(MercurvM J2000Ee) 
(RxV)JVlag_(MoonFixed) 

(RxV) JVIag_(MoonM J2000) 
(RxV) _Mag_(MoonMJ2000Ec) 
(RxV) _Mag_(N eptuneFixed) 
(RxV) JVIag_(Ncpt uboM J2000) 
(RxV) _Mag_(N cptuncM J 2000Ec) 
(RxV) Jdag.(PhitoFixed) 

(RxV) _Mag_(Plu toM J2000) 
(RxV) _Mag_(PliitoM J2000Ec) 
(RxV) Jvtag_(SaturiiFixed) 

(RxV) _Mag_(SatnmM J2000) 
(RxV)Jdag.(SatumMJ2000Ec) 
(RxV) _Mag_(UramiaFixed) 

(RxV) _Mag_(UraimsMJ2000) 
(RxV) _Mag_(TJrauussM J2000Ec) 
(RxV) _Mag_(VenusFixsd) 

(RxV) _Mag_(VeminM J2000) 
(RxV) _Mag_(VermsM J2000EC) 
(RxV)JC 
(RxV) JC_(Fixed) 

(RxV) _X_(M J2000Ec) 
(RxV)JC.(MODEc) 

CRxV) GL(MODEq) 
(RxV)JL(MOEEq) 

(RxV) JC-(MaxsFixed) 

(RxV) Jt_(MaisM J2flnOEc) 

(RxV) JL(MarsM J2000Eq) 

(RxV) _X_(MercuryFixed) 

(RxV) JS-(Merc«ryMJ2000Ec) 
(RxV)_X_(MercuryMJ2000Eq) 
(RxV) -X-(MooiiFixcd) 

(RxV) _X_(MoouM J 2000Ec) 
(RxV) JC_(MoonM J2000Eq) 
(RxV) _X_(NeptuueFixed) 

(RxV) JC.(NeptuneMJ2000Ec) 
(RxV) _X_(NeptuneM J2000Eq) 
(RxV) _X_(PlutoF ixed) 

(RxV) JC.(PlutoM J2000Ec) 
(RxV) _X_(PhitoMJ2000Eq) 
(RxV) _X_(SaturnFixcd) 

(RxV) _X_(S aturnM J2000Ec) 
(Rx\0-K-(SaturnAtJ2000Eq) 
(RxV) JX.(EarthGSE) 

(RxV) _Y ..(Earth GSB) 
(RxV)JZ_(EarthGSE) 
EartliGSM-PositiondX 


(RxV) JZ.(M OEEq) 

(RxV) JZ.(MarsFixed) 

(RxV) 2ZL(MarsMJ2000Ec) 
(RxV) .Z.fMaxsM J2000Eq) 
(RxV) -Z-(MercuryFixed) 

(RxV) .Z-(MercuryMJ2000Ec) 
(RxV)-Z.(MercnryMJ2000Eq) 
(RxV)_Z_(MoonFixed) 

(RxV) -Z_(MoonMJ2000Ec) 
(RxV) -ZJMoonM J2000Eq) 
(RxV)JZ.(NeptmieFixed) 
(R>:V)-Z-(NcptuncMJ2000Ec) 
(RxV) _Z_(N cptuncM J 2000Eq) 
(RxV) -Z.(PlutoFixed) 

(RxV) _Z_(PIutoMJ2000Ec) 
(RxV) _Z_(PlutoM J2000Eq) 
(RxV) .Z.(Sat'uriiFixed) 

(RxV) _Z _ (S atiirnM J2000Ec) 

(Rx V j _Z . ( S atuniM .J 2000E< j) 

(RxV)JZ.(TODEc) 

(RxV)-Z.(TODEq) 

(RxV) -Z.(TOEEc) 
(RxV)-Z.(TOEEq) 

(RxV) .Z.(UramisFixed) 
(RxV)_Z_rUrannfiMJ2000Ec) 
(RxV) _Z_ (U rani isM J2000Eq) 
(RxV) JZ.(VenusFixed) 
(RxV)-Z_(VeiiusMJ2000Ec) 
(RxV) JZ - ( VeirnsM J2000Eq ) 
(RxV) _X_(MOEEc) 
Vector_Dec_(Fixed) 

Vector J3ec_(MJ2000Ec) 

Vector JDec_(MODEc) 

Vector JDec_(MODEq) 
Vector_Dec_(MOEEc) 

Vector JDec.(MOEEq) 

Vector JDec.(MaraFixed) 

Vector _Dcc_(MarsMJ2000Ec) 
Vector JDec_(MarsMJ2000Eq) 
Vector JDee_(Merci;ryFixed) 
Vector JDec_(MercuryMJ2000Ec) 
Vector J>io_(MercuryMJ2000Eq) 
Vector _Dec_(MoonFixed) 
Vector-Dec. (MoonMJ2000Ec) 
Vector JDec_(MooiiMJ2000Eq) 
Vector_Dec_(NeptuneFixed) 
Vector_Dec_(NcptuiioMJ2000Ec) 
VectorJ3ec_(NeptuneMJ2000Eq) 
(RxV) -Z.(MOEEc) 
(RxV)-K-(EarthGSM) 

(RxV) _Y_(EarthGSM) 
(RxV)-Z.(EarthGSM) 
EarthGSM_Velocity_X 
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EaxthGSMJPositiou_Y 
EarthGSM_Position_Z 
Vector JDec.(EarthGSE) 

Vector .RA.(EartkGSE) 

Vector JDec_(PlutoFixed) 

Vector JDec.(PlutoMJ2000Ec) 
VectorJDec_(PlutoMJ2000Eq) 
Vector _Dec_(SatumFixed) 
Vector JDcc_(SaturnMJ2000Ec) 
Vector _Dec_(SatuniMJ2000Eq) 
Vector JDec-(TODEc) 

Vector 2Dec_(TODEq) 

Vector JDec_(TOEEc) 

Vector JDec.(TOEEq) 
Vector_Dec_(UraimsFixed) 
Vector JDec_(UramisMJ2000Ec) 
Vector JDec-(UranusMJ2000Eq) 
Vector .Dcc-(VciiusFixcd) 

Vector JDee_(VeniisMJ2000Ec) 
Vector _Dec_(VeausMJ2000Eq) 
Vector JRA-(Fixed) 

Vector _RA_(MJ20Q0Ec) 

Vector -RA-(MODEc) 

Vector JRA-(MODEq) 

Vector _R A- (MO EEc) 

Vector JR.Ai(MOEEq) 

Vector _RA_(MarsFixcd) 

Vector JEtA_(MaxsM J 2000Ec) 
VectorJEtA_(Mar,sMJ2000Eq) 
Vector JRA_(MerairyFixed) 
Vector JELA_(MercirryMJ2000Ec) 
Vector _R A. (MercaryM J2000Eq) 
Vector JR A_(MoonFixed) 

Vector JEtA-(MoouMJ2000Ec) 
Vector _RA_(MoonMJ2000Eq) 
Vector _RA-(N eptuueFixed) 
Vector _RA.(NeptuneMJ2000Ec) 
Vector JRA.(NeptnmeMJ2000Eq) 
Vector JR_A_(PlutoFixed) 

Vector JR. A_(PlutoMJ2000Ec) 
Vector JRA_(PlutoM.J2000Eq) 
Vector _RA.(Sa,turnFixed) 
Vcctor_RA..(SaturnMJ2000Ec} 
Vector -RA-(SafcumMJ2000Eq) 
Vector JElA-(TODEc) 

Vector JEtA-(TODEq) 

Vector JRA.(TOEEc) 

Vector JtA_(TOEEq) 

Vector JtA.(UranusFixed) 

Vector _RA_(UraimsMJ2000Ec) 
Vec tor _R A_(Ur aousM J 200GEq) 
Vector _R A. ( VcnusFixcd) 

Vector JEtA_(VenusMJ2000Ec) 
Vector .R A. (VenusMJ2000Eq) 


Earth GSM_Vdocity_Y 

EartliGSM.Velocity_Z 

Vector_Dec_(EartliGSM) 

Vector _R. A. (Ear thGSM) 

(RxV)JC(TODEc) 

(RxV)-X-(TODEq) 

(RxV)-X-(TOEEc) 

(RxV)JC.(TOEEq) 

(RxV) -X-(UrauusFixcd) 

(RxV) _X_(UraiuisMJ2000Ec) 
(RxV)-X-(Uraau»MJ2000Eq) 
(RxV) _X_( VemtsFixed) 

(RxV) _K_( VenusM J 2000Ec) 
(RxV) JC_( VennsM J 2000Eq) 
(RxV)-Y 
(RxV) -Y-(Fixed) 

(RxV) _Y_(M J 2000Ec) 

(RxYi-Y.(MODEc) 

(RxV)-Y-(MODEq) 

(RxV)_Y_(MOEEc) 

(RxV).Y.(MOEEq) 

(RxV) _Y_(M arsFixed) 
(RxV)_Y-(MarsM,T200GEc) 
(RxV)_Y_(MarsM'J2000Eq) 
(RxV) -Y-(MercuryFixed) 
(RxV) _Y-(McrcuryM J2000Ec) 
(RxV) _Y_(McrcmyM J 2000Eq) 
(RxV) _Y_(MooiiFixed) 

(RxV) _Y_(Moom\U2000Ec) 
(RxV) _Y_(MoonM J2000Eq) 
(RxV)_Y_(NeptuneFLxed) 
(RxV) _Y.(N eptimeM J2000Ec) 
CRxV) _Y_(N eptuneM J2000Eq) 
(RxV) _Y_(PlutoFbced) 

(RxV) _Y_(PlutoM J 2000Ec) 
(RxV) _Y_(PlutoM J 2000Eq) 
(RxV) _Y_ (S atumFixed) 
(RxV)_Y_(SaturnMJ2000Ec) 
(RxV).Y.(SatiiniMJ2000Eq) 
(RxV) .Y.(TODEc) 
(RxV).Y.(TODEq) 

(RxV) .Y-_(T O EEc) 

(RxV) _Y_(T OEEq) 
(RxV).Y_(UrauuaMJ2000Ec) 
(RxV) _Y_(Ur aauaM J 2000Eq) 
(RxV).Y-(UrajaussFixed) 
(RxV)-Y-(VenufiFixed) 

(RxV j _Y_C VenusM J2000Ec) 

(RxV).Y_(VenuaMJ2000Eq) 

(RxV)2Z 

(RxV)_Z_(Fixed) 

(RxV).Z.(MJ2000Ec) 

(RxV) -Z.(MODEc) 
(RxV)-Z_(MODEq) 
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D.4.8 


D.4.9 


D.4.10 


D.4.11 


Coordinate Systems: Central Body Inertial Elements 


Mars_CenteraLMean_J2000 

Mars.Centered_MJ2000Ee 

Merc’ary.Geiitered JwTeiin_ J 2000 

Mercury -CenteredJVf J2000Ec 

Moon-Centered Jdean. J2000 

Moon_Centered_MJ2000Ec 

Neptune .Graitcrcd.Mcaoa_ J 2000 

Ncptunc_Ccntcrcd_MJ2000Ec 

Pluto_Centered_Meau.J2000 

Pluto_Centered_M J 2000Ec 

Satiim-Centered_Meaii_J2000 

Satum_Centered_MJ2000Ec 

Uranus_Centered_Mean_J2000 

Uranixs_Cenfcered_MJ2000Ec 

Venus_CenteredJV[ean_J2000 

Sune_Centered_MJ2000Ec 

Vem is_Ceiitered_M,T 2000Ec 

Propagators Elements 

Moon_LPlG5P 

MoonJLP165P-and_SRP 

CisLunar 

Moon_SRP 

EarthMoon_L2_AllPlanets 

Neptune.2-Body 

EarthMoonA2_AllPlaaets_and_SRP 

' Neptuue_AllPlanets 

EarthMoon_L2JESM 

Neptmie_SRP 

Ear fchMo ojj_L2 JBSM_J GM 

Phito_2-Body 

EarthMoon _L2 JESM API 65P 

Pluto-AllPlanets 

EarthSunA2_AllPlanets 

Pluto-SRP 

EarthSimJL2-Al]Planefcs_and_SRP 

8aturn_2-Body 

Earth. J2 

Saturn_AiIPlanets 

Earth.Point.Mass 

Satum_SRP 

Heliocentric 

Snn_AllPlauets 

Lunar 

Sun_AllPlanets.and.SRP 

Maxs_2-Body 

Uranus _2-Body 

Mars_AJlPlancts 

Uranus -AilPlancts 

Mars.MARSo0C 

Uranus -SRP 

Mars-MARS50C.and.SRP 

Venus _2-Body 

MarsJSRP 

Venus -AllPlanets 

Mercury _2-Body 

Venus JV1GNP 1 80 U 

Merc.ury_AllPlanets 

Venus_MGNP180U_iincLSRP 

Mercury _SRP 

Venus _SRP 

Moon_2-Body 

Moon-AilPIanets 


Axes Elements 

Earth_GSE Earth.GSM 


Coordinate System Elements 


Earth.GSE Earth.GSM 
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D.4.12 Vectors Elements 


Cross _Prc duct _(Rx V) 

Cross-Pro duct _ (Rx V) _ (fixed) 

Cross -Pro duct _(Rx V) _(M arsFixod ) 

Cross JProduct_(Rx' V) _(MaxsM J 2000) 

Cross _Product.(RxV)_(MaxsMJ2O0OEc) 
Cross _Produet_(RxV)_(MercuryFixed) 
Cross ..Product_(RxV).(MercuryMJ2000) 
Cross JProduct_(RxV)_(MercuryMJ2000Ec): 
Cross JProduct-(RxV) .(MoonFixed) 

Cross -Product.(RxV)_(MoouMJ2000) 

Cross JProduct-(R.xV)-(MooiiMJ2000Ec) 
Cross-Product. (RxV).(NcptuncFixcd) 
Cross_Product_(RxV)_(NeptuneMJ2000) 
Cross-Product_(RxV)-(NeptuneMJ2000Ec) 
Cross_Product.(YofGSM) 

EarthSun.J 2K_Velocity 
DiPole. Vector.GSM 


Cross_P ro duct _(Rx V) _ (Plut oFixed) 

Cross_Product_(RxV)_(PlutoMJ2000) 

Cross-Product_(RxV)-(PlutoMJ2000Ec) 

Cross JProduct_(RxV)_(SaturnFixed) 

Cross_Product_(RxV).(SaturnMJ2Q00) 

Cross_Product-(RxV)_(SaturnMJ2000Ec) 

Cross-Product ..(RxV).(UranusFixed) 

Cross_Product_(RxV)_(UramisMJ2000) 

Cross-Product_(RxV)_(UranusMJ2000Ec) 

Cross_Product.(RxV)..('VeimsFixed) 

Cross JPro duct.(RxV) .( VeuusM J 2000) 
Cross_Product_(RxV)_(VcuusMJ2000Ec) 
Cross _Product_(RxV)_(EaxthGSE) 

Cross _Product_(RxV)-(Eart.hGSM) 

Cross JProduct_(ZofGSE) 
EarthSun.Vector 


D.4.13 Vectors: Vehicle Local Elements 


Position_(MarsFixed) 

Position_(MarsMJ2000Ec) 

P osition_(MercuryFixed) 

Positiou_(MercuryMJ2000Ee) 

Positiou_(MoouFixcd) 

P ositiou_(MoonM J 2000Ec) 

Position_(NeptuueFixed) 

Position_(NeptuneMJ2OO0Ec) 

Position.(PlutoFixed) 

Position.(Plut,oMJ2000Ec) 

Position_(8atm:nFixed) 

Position_(SatumMJ2000Ec) 

Position.(UranusFixed) 

Positiou..(UrauusMJ2000Ec) 

Position- (VenusFixed) 

P osition_( VenusM J 2000Ec) 

Position._Maxs_(MJ2000) 

Position_Mercury_(MJ2000) 

Position_Moon_(MJ2000) 

Position -Neptune. (MJ2000) 

Pcsition_Sat urn_(M J 2000) 

Positioii_Urauus_(M.T2000) 

P ositiou-Vcnus. (M J2000 ) 

Position-(EarthGSE) 

Positk>n.(Eartb.GSM) 


Velocity.(MarsFixed) 

Velocity.(MarsMJ2000Ee) 

Velocit y _ (Mere! uyF ixed) 

Velocity_(Merc\iryMJ2000Ec) 

Vclocity.(MoonFixcd) 

Vdocity.(MoonMJ2000Ec) 

Velocity.(NeptutteFbced) 

Veloeity.(NeptrineMJ2Q00Ec) 

Veloci ty.(PlutoFixed) 

Velocity_(Plut.oMJ200QEc) 

Velotity.(SaturnFixed) 

Velocity. (S atumM J 2000Ec) 

Vclocity.(UranusFixcd) 

Vok)city.(Ur;xnusMJ2000Ec) 

Velocity.(VenusFixecl) 

Velocity.! VenusM J 2000Ec) 

Velocity _Mars_(MJ2000) 

Velocity J\ , Iercury_(MJ2000) 

Velocity _Moon_(MJ2000) 

Velocity _N eptune_(MJ 2000) 

Velocity_Suturii.(MJ2000) 

Velocity _Uranus_(MJ2000) 

Vclocitj’--Vciius.(MJ2000) 

Velocity.(EarthGSE) 

Veloci ty. (Ear thGSM) 
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D.5 Report Styles 





Classical_Orbit_Elements 

EarthJVIJ2000_Position_Velocity 

GMAT_ApoapsisJPeriapsis 

GM AT _C SP axameters J? bred 

GMAT_CSParameters_MJ2000Ec 

GMAT-GSP ar amctcrsJM J2000Eq 

GMAT_CSParamctcrs_MODEc 

GMAT.CSP axameters JM ODEq 

GM AT _CbP axameters 

GMAT.CbP arameters_Mars 

GMAT.CSPaxametersJVI.OEEq 

GMAT_GSParameters_MarsFixed 

G'MAT _GSP ar ametexs JMaxts Jil J2000Ec. 

GMAT-CSParametersJVtars-MJSOOOEq 

GMAT-OSPaxamctcrsJMercuryFixod 

GMAT_CSParametera_Mercury_MJ2()00Ee 

GMAT_CSParameters_Mercury_MJ2000Eq 

GMAT.CSPaxameters.JMooiiF.ixed 

GMAT_CSPaxametersJMoon_MJ2000Ec 

GMAT_CSPaxametersJMoonJMJ2000Eq 

GMATX'SParameters-N eptuneFixed 

GMAT _CSP axameters JN T eptune JM J 2000Ec. 

GMAT _CSPar ameters _N optuncJMJ 2000Eq 

Pluto -MJ2000 .FositiomVclocity 

SatumJMJ2000JPositioxi_Velocity 

Sim JvlJ2000EeJ > osition_ Velocity 

Sun .MJ2000 -Position-Velocity 

Uranus _MJ2000_Position_ Velocity 

Venus Jv02000 -Position-Velocity 

Gxeeuwicli-Hour .Angle 



GMAT_CSParameters_TODEq 
GMAT_CSParameters_TOEEc 
GMAT_CSParameters_TOEEq 
GM AT_C SP ar ame ters _U rauus Fixed 
GMAT_CSParaineters_Uraxiiis-MJ20G0Ec 
GMAT_CSParamctcrs_UxanusJVlJ2000Eq 
GM AT .0 SP ar ainctcrs .V onusFixed 
GMAT_CSParameters_Venus-MJ2000Ec 
GMAT-GSParameters-VenusJVIJ2000Eq 
GMAT-GSP arameters-MOEEc 
GMAT-CbPar ameters .Mercury 
GMAT-CbPar ameters JVIoon 
GM AT _CbP ar ameters -Neptune 
GMAT_CbParaineters -Pluto 
GMAT-CbParametcrs-Saturu 
GMAT-CbP axameters -Uranus 
G MAT-C bP axameters -Venus 
J2000-E CLPosif, ion. Velocity 
Mars-MJ2000JPosition.Velocity 
Merc iiry_M J2Q00 JPosi tion.Veloci ty 
Moon-MJ2000 -Position-Velocity 
N epti me_M J 2000 -Position.' Velocity 
GMAT-CSParamctcrsJPlutoFixcd 
GMAT.GSParamctcrsJPIuto_MJ2000Ec 
GMAT-CSPaxameters_Pluto-MJ2000Eq 
GMAT-CSParameters-SaturnFixed 
GMAT.CSPar ameters -Saturn JVI,T2000Ec 
GMAT_CSParameters_SatumJVlJ2000Eq 
GM AT-CSP arameters.T ODEc 


D.6 Scripts Used 


The main script used with STK is the STK-Repropagate.m Matlab script. STK_Repropagate.m connects 
to STK, propagates satellites in the scenario, generates reports, and outputs performance run times. The 
STK_Repropagate script was designed to reduce the time it took to generate STK report files, after modifica- 
tions to the STK scenario were made, and obtain accurate STK ran times. 

Refer to Appendix C for more details of this script and others used in the Acceptance Test Plan document. 
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Appendix E 


FF Setup 

[INSERT explanation of how the FF scenarios were setup] 
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